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Pek S PREFACE 


This book is a concise introduction to modern probability 
theory and certain of its ramifications. By deliberate 
succinctness of style and judicious selection of topics, it 
manages to be both fast-moving and self-contained. 

The present edition differs from the Russian original 
(Moscow, 1968) in several respects: 


1. It has been heavily restyled with the addition of 
some new material. Here I have drawn from my own 
background in probability theory, information theory, 
etc. 


2. Each of the eight chapters and four appendices has 
been equipped with relevant problems, many accom- 
panied by hints and answers. There are 150 of these 
problems, in large measure drawn from the excellent 
collection edited by A. A. Sveshnikov (Moscow, 1965). 


3. At the end of the book I have added a brief 
Bibliography, containing suggestions for collateral and 
supplementary reading. 

R.A. S. 
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BASIC CONCEPTS 


I. Probability and Relative Frequency 


Consider the simple experiment of tossing an unbiased coin. This 
experiment has two mutually exclusive outcomes, namely “heads’’ and 
“tails.” The various factors influencing the outcome of the experiment are 
too numerous to take into account, at least if the coin tossing is “fair.” 
Therefore the outcome of the experiment is said to be “random.” Everyone 
would certainly agree that the “probability of getting heads” and the “prob- 
ability of getting tails’ both equal $. Intuitively, this answer is based on the 
idea that the two outcomes are “equally likely” or “equiprobable,” because of 
the very nature of the experiment. But hardly anyone will bother at this 
point to clarify just what he means by “probability.” 

Continuing in this vein and taking these ideas at face value, consider an 
experiment with a finite number of mutually exclusive outcomes which are 
equiprobable, i.e., “equally likely because of the nature of the experiment.” 
Let A denote some event associated with the possible outcomes of the 
experiment. Then the probability P(A) of the event A is defined as the fraction 
of the outcomes in which A occurs. More exactly, 


— 


P(A) = (1.1) 


where N is the total number of outcomes of the experiment and N(A) is the 
number of outcomes leading to the occurrence of the event A. 


Example 1. In tossing a well-balanced coin, there are N = 2 mutually 
exclusive equiprobable outcomes (“‘heads” and “‘tails”). Let A be either of 
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these two outcomes. Then N(A) = 1, and hence 


Example 2. In throwing a single unbiased die, there are VN = 6 mutually 
exclusive equiprobable outcomes, namely getting a number of spots equal 
to each of the numbers | through 6. Let A be the event consisting of getting 
an even number of spots. Then there are N(A) = 3 outcomes leading to the 
occurrence of A (which ones ?), and hence 

ae 


Example 3. In throwing a pair of dice, there are N = 36 mutually 
exclusive equiprobable events, each represented by an ordered pair (a, 5). 
where a is the number of spots showing on the first die and 5 the number 
showing on the second die. Let A be the event that both dice show the same 
number of spots. Then A occurs whenever a = J, 1.e., n(A) = 6. Therefore 

ee 
P(A) = ae & 

Remark. Despite its seeming simplicity, formula (1.1) can lead to 
nontrivial calculations. In fact, before using (1.1) in a given problem, we 
must find all the equiprobable outcomes, and then identify all those leading 
to the occurrence of the event A in question. 


The accumulated experience of innumerable observations reveals a 
remarkable regularity of behavior, allowing us to assign a precise meaning 
to the concept of probability not only in the case of experiments with equi- 
probable outcomes, but also in the most general case. Suppose the experi- 
ment under consideration can be repeated any number of times, so that, in 
principle at least, we can produce a whole series of “‘independent trials under 
identical conditions,’’! in each of which, depending on chance, a particular 
event A of interest either occurs or does not occur. Let n be the total number 
of experiments in the whole series of trials, and let n(A) be the number of 
experiments in which A occurs. Then the ratio 


n(A) 


n 


is called the relative frequency of the event A (in the given series of trials). It 
turns out that the relative frequencies n(A)/n observed in different series of 


* Concerning the notion of independence, see Sec. 6, in particular footnote 2, p. 31. 
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trials are virtually the same for large n, clustering about some constant 


A 
fae (1.2) 
called the probability of the event A. More exactly, (1.2) means that 
A 
P(4) = lim", (1.3) 
n>o N 


Roughly speaking, the probability P(A) of the event A equals the fraction of 
experiments leading to the occurrence of A in a large series of trials.” 


_ Example 4. Table 1 shows the results of a series of 10,000 coin tosses,? 
grouped into 100 different series of nm = 100 tosses each. In every case, the 
table shows the number of tosses n(A) leading to the occurrence of a head. 
It is clear that the relative frequency of occurrence of “heads” in each set 
of 100 tosses differs only slightly from the probability P(A) = § found in 
Example 1. Note that the relative frequency of occurrence of “heads” is even 
closer to 4 if we group the tosses in series of 1000 tosses each. 


Table 1. Number of heads in a series of coin tosses 


Number of heads Number of heads 
in 100 series of in 10 series of 
100 trials each 1000 trials each* 
54. 46-33. 55 46 34 4) 48 51 53 501 
48 46 40 53 49 49 48 54 53 45 485 
43-9 er ot 50 S22 50 S53 C49 509 
Se oe > 0 «6450 AT O57 52 55 536 
48 51 51 49 44 52 50 46 53 41 485 
49-50 > 52° 52° 48 #47 = 47 «47 ~«(5i 488 
a5 47a ot Ae. 50 55: 53 50 500 
53 52 46 52 44 51 48 51 46 54 497 
45 47 46 52 47 48 59 57 45 48 494 
47° 41 734.- 468 Se 63513: 52. 55 39 41 484 


Example 5 (De Mere’s paradox). As a result of extensive observation of 
dice games, the French gambler de Méré noticed that the total number of 
spots showing on three dice thrown simultaneously turns out to be 11 (the 
event A,) more often than it turns out to be 12 (the event A,), although 
from his point of view both events should occur equally often. De Meéré 


2 For amore ee discussion of the meaning of (1.2) and (1.3), see Sec. 12 on the “law 
of large numbers.” 

3 Table 1 is taken from W. Feller, An Introduction to Probability Theory and Its Appli- 
cations, Volume I, third edition, John Wiley and Sons, Inc., New York (1968), p. 21, and 
actually stems from a table of “random numbers.” 

4 Obtained by adding the numbers on the left, row by row. 
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reasoned as follows: A, occurs in just six ways (6:4:1, 6:3:2, 5:5:1, 5:4:2, 
5:3:3, 4:4:3), and A, also occurs in just six ways (6:5:1, 6:4:2, 6:3:3, 
5:5:2, 5:4:3, 4:4:4). Therefore A, and A, have the same probability 
P(A,) = P(A)). 

The fallacy in this argument was found by Pascal, who showed that the 
outcomes listed by de Méré are not actually equiprobable. In fact, one must 
take account not only of the numbers of spots showing on the dice, but also 
of the particular dice on which the spots appear. For example, numbering 
the dice and writing the number of spots in the corresponding order, we find 
that there are six distinct outcomes leading to the combination 6:4:1, namely 
(6, 4, 1), (6, 1, 4), (4, 6, 1), (4, 1, 6), C1, 6, 4) and (1, 4, 6), whereas there is 
only one outcome leading to the combination 4:4:4, namely (4, 4, 4). The ap- 
propriate equiprobable outcomes are those described by triples of numbers 
(a, b, c), where a is the number of spots on the first die, b the number of spots 
on the second die, and c the number of spots on the third die. It is easy to 
see that there are then precisely N = 6* = 216 equiprobable outcomes. Of 
these, N(A,) = 27 are favorable to the event A, (in which the sum of all the 
spots equals 11), but only N(A,) = 25 are favorable to the event A, (in which 
the sum of all the spots equals 12).° This fact explains the tendency observed 
by de Méré for 11 spots to appear more often than 12. 


2. Rudiments of Combinatorial Analysis 


Combinatorial formulas are of great use in calculating probabilities. 
We now derive the most important of these formulas. 


THEOREM 1.1. Given n, elements a,, d.,..., a, and nz elements by, 
ba,...,6,., there are precisely nynz 

° ° 2 e * : 
distinct ordered pairs (a;, b;) contain- 
ing one element of each kind. 


Proof. Represent the elements of 
the first kind by points of the x-axis, 
and those of the second kind by points 
of the y-axis. Then the possible pairs 
(a;, b;) are points of a rectangular 
lattice in the xy-plane, as shown in 
Figure 1. The fact that there are 
just nn, such pairs is obvious from 
FiGure 1. the figure. [f° 


> To see this, note that a combination a:b:c occurs in 6 distinct ways if a, 6 and c are 
distinct, in 3 distinct ways if two (and only two) of the numbers a, b and c are distinct, and 
in only 1 way if a = 6 = c. Hence A, occurs in6 + 6+ 3+6+3 + 3 = 27 ways, while 
A, occurs in6+6+3+3-+6-+ 1 = 25 ways. . 
6 The symbol §f stands for Q.E.D. and indicates the end of a proof. 
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More generally, we have 


THEOREM 1.2. Given n, elements aj, dz,...,Qn,5 Ne elements by, 
Oar. «1, 0,,, C1C., Up ton, elements X1, X2,---»Xn,» there are precisely 
nyn,*+*n, distinct ordered r-tuples (a;,,5;,,...,%;,) containing one 


element of each kind." 


Proof. For r = 2, the theorem reduces to Theorem 1.1. Suppose 
the theorem holds for r — 1, so that in particular there are precisely 


n,* +n, (r — 1)-tuples (b;,,..., x;,) containing one element of each 
kind. Then, regarding the (r — 1)-tuples as elements of a new kind, we 
note that each r-tuple (a; ,, b;,,..., x;,) can be regarded as made up of 
a (r — 1)-tuple (5;,,..., x;,) and an element a,. Hence, by Theorem 


1.1, there are precisely 
Ny(Np***N,) = NyNg*** Ny 5 


r-tuples containing one element of each kind. The theorem now 
follows for all r by mathematical induction. J 


Example 1. What is the probability of getting three sixes in a throw 
of three dice? 


Solution. Let a be the number of spots on the first die, b the number of 
spots on the second die, and ¢ the number of spots on the third die. Then 
the result of throwing the dice is described by an ordered triple (a, 5, c), 
where each element takes values from 1 to 6. Hence, by Theorem 1.2 with 
r = 3 and ny, =n, = ng = 6, there are precisely N = 6° = 216 equiprobable 
outcomes of throwing three dice (this fact was anticipated in Example 5, 
p. 3). Three sixes can occur in only ‘one way, i.e., when a=b=c= 6. 
Therefore the probability of getting three sixes is 37. 


Example 2 (Sampling with replacement). Suppose we choose r objects 
in succession from a “population” (i.e., set) of n distinct objects a,, dg,..., 
a,,, in such a way that after choosing each object and recording the choice, we 
return the object to the population before making the next choice. This 
gives an “‘ordered sample’ of the form 


(G5 Uys e.g A}. (1.4) 
Setting 2, = m = -«:* =n; = nin: Theorem. 1:2,: we find that: there. are 
precisely 
A yt Eo 
distinct ordered samples of the form (1.4).° 
7 Two ordered r-tuples (@;;, big, ... , Xin) and (4;,, bjs, ... , Xj,) are said to be distinct 
if the elements of at least one pair a@;, and a;,, b;, and bj,,..., a;, and b;, are distinct. 
8 Two “ordered samples” (4;1, in, ... , @:,) and (@;1, Aja, ... , @;,) are said to be distinct 
if a,, A aj, for at least one kK = 1,2,...,1r. This is a special case of the definition in 


footnote 7. 
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Example 3 (Sampling without replacement). Next suppose we choose r 
objects in succession from a population of n distinct objects ay, do, ... , Ans 
in such a way that an object once chosen is removed from the population. 
Then we again get an ordered sample of the form (1.4), but now there are 
n — | objects left after the first choice, n — 2 objects left after the second 
choice, and so on. Clearly this corresponds to setting 


n= Nn, Ne =—an—1,..., 2p =H 


in Theorem 1.2. Hence, instead of n” distinct samples as in the case of sam- 
pling with replacement, there are now only 


N=n(n — 1)°+:(n—r+]1) (1.6) 
distinct samples. If r = n, then (1.6) reduces to 
N= nie = te 2 Se (1.7) 


the total number of permutations of n objects. 


Example 4. Suppose we place r distinguishable objects into n different 
“cells” (r < n), with no cell allowed to contain more than one object. Num- 
bering both the objects and the cells, let 7, be the number of the cell into which 
the first object is placed, i, the number of the cell into which the second 
object is placed, and so on. Then the arrangement of the objects in the cells 
is described by an ordered r-tuple (i,, i,,...,i,). Clearly, there are n, =n 
empty cells originally, n, = n — 1 empty cells after one cell has been occupied, 
Ng = n — 2empty cells after two cells have been occupied, and so on. Hence, 
the total number of distinct arrangements of the objects in the cells is again 
given by formula (1.6). 


Example 5. A subway train made up of n cars is boarded by r passengers 
(r <n), each entering a car completely at random. What is the probability 
of the passengers all ending up in different cars? 


Solution. By hypothesis, every car has the same probability of being 
entered by a given passenger. Numbering both the passengers and the cars, 
let i, be the number of the car entered by the first passenger, 7, the number 
of the car entered by the second passenger, and so on. Then the arrangement 
of the passengers in the cars is described by an ordered r-tuple (i, iz, . . . . i,), 
where each of the numbers i, i,,...,i, can range from 1 to n. This is 
equivalent to sampling with replacement, and hence, by Example 2, there are 


N= 


distinct equiprobable arrangements of the passengers in the cars. Let A 
be the event that “no more than one passenger enters any car.’’ Then A 
occurs if and only if all the numbers 7,, i.,...,i, are distinct. In other 
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words, if A is to occur, the first passenger can enter one of n cars, but the 
second passenger can only enter one of n — I cars, the third passenger one 
of n — 2cars, and so on. This is equivalent to sampling without replacement, 
and hence, by Example 3, there are 


NC = na 2) @ Sr) 
arrangements of passengers in the cars leadin g to the occurrence of A. There- 
fore, by (1.1), the probability of A occurring, i.e., of the passengers all ending 
up in different cars, is just 


ee ee) 


i 


n 


Any set of r elements chosen from a population of n elements, without 
regard for order, is called a subpopulation of size r of the original population. 
The number of such subpopulations is given by 


THEOREM 1.3. A population of n elements has precisely 
i (1.8) 
r'(n—r)! 
subpopulations of sizer < Nn. 


Proof. If order mattered, then the elements of each subpopulation 
could be arranged in r! distinct ways (recall Example 3). Hence there 
are r! times more “ordered samples’’ of r elements than subpopulations 
of size r. But there are precisely n(n — 1)--+ (nm — r + 1) such ordered 
samples (by Example 3 again), and hence just 

n(n —1)---(n—r+1) _ n! 
r! r!(n—r)! 
subpopulations of sizer. ff 
Remark. An expression of the form (1.8) is called a binomial coefficient, 
often denoted by 
(" 
r 


instead of C”. The number C,’ is sometimes called the number of combinations 
of n things taken r at a time (without regard for order). 


The natural generalization of Theorem 1.3 is given by 


THEOREM 1.4. Given a population of n elements, let ny, n2,..., nN, be 
positive integers such that 


Ny tng +e + +m =n. 


8 
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Then there are precisely 


n!} 
A Se (1.9) 
Ay hg) a 
ways of partitioning the population into k subpopulations, of sizes 
Ny, Ng,... Nx, respectively. 


Proof. The order of the subpopulations matters in the sense that 
nS 2, eae eS ae a oe Ny = 2,N3,..., My, (say) 
represent different partitions, but the order of elements within the 
subpopulations themselves is irrelevant. The partitioning can be effected 
in stages, as follows: First we form a group of n, elements from 
the original population. This can be done in 


Ny oe S. 


ways. Then we form a group of n, elements from the remaining n — n, 
elements. This can be done in 
Nz = Se 
ways. Proceeding in this fashion, we are left withn — ny — ++: —n,_. = 
ny,_1 + m, elements after k — 2 stages. These elements can be partitioned 
into two groups, one containing n,_, elements and the other n, elements, 
in 
ING = Shs a ie 


ways. Hence, by Theorem 1.2, there are 
N = N,No: + Nyy 


mas if Rae fl ae 
=e Ben oe Cg 


distinct ways of partitioning the given population into the indicated k 
subpopulations. But 


Se "aie | a Cer ed 
1 Me 


Ne-1 
re n! (n —n,)! = (n — n, — *-- =A! 
= (n — n,)! ng! (n — n, —n,)! n,_,!(n — ny —***— nh ge — a) 
= n! (n — n,)! aia tO 
ny! (n — ny)! ng! (n — ny — ny)! Nal BS! 
n! 
es 


in keeping with (1.9). J 
Remark. Theorem 1.4 reduces to Theorem 1.3 if 


k= 2,> = j-th oe 


SEC. 2 BASIC CONCEPTS 9 


The numbers (1.9) are called multinomial coefficients, and generalize the 
binomial coefficients (1.8). 


Example 6 (Quality control). A batch of 100 manufactured items is 
checked by an inspector, who examines 10 items selected at random. If 
none of the 10 items is defective, he accepts the whole batch. Otherwise, the 
batch is subjected to further inspection. What is the probability that a batch 
containing 10 defective items will be accepted? 


Solution. The number of ways of selecting 10 items out of a batch of 
100 items equals the number of combinations of 100 things taken 10 at a time, 
and is just 
100! 


Ws CA 
= $91 O01 - 


By hypothesis, these combinations are all equiprobable (the items being 
selected ‘‘at random’’). Let A be the event that “the batch of items is accepted 
by the inspector.” Then A occurs whenever all 10 items belong to the set of 
90 items of acceptable quality. Hence the number of combinations favorable 
to A is 

90! 
10! 80! 


It follows from (1.1) that the probability of the event A, i.e., of the batch 
being accepted, equals® 


N(A) = Cio = 


0) ee 2 


N(A) 90! 90! S182 | i) = 
N 80!100! 91-92---100 e° 


where e = 2.718... is the base of the natural logarithms. 


Example 7. What is the probability that two playing cards picked at 
random from a full deck are both aces? 


Solution. A full deck consists of 52 cards, of which 4 are aces. There are 


52 S21] 
_ = 


2! 50! 


== 1326 


ways of selecting a pair of cards from the deck. Of these 1326 pairs, there are 


® The symbol ~ means “is approximately equal to.” 
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consisting of two aces. Hence the probability of picking two aces is just 


Example 8. What is the probability that each of four bridge players 
holds an ace? 


Solution. Applying Theorem 1.4 with n = 52 andn, = n, = ng = 1g = 
13, we find that there are 
52! 
131 13013013! 
distinct deals of bridge. There are 4! = 24 ways of giving an ace to each 
player, and then the remaining 48 cards can be dealt out in 
48! 
a As Pee ay 
distinct ways. Hence there are 
48! 
(12!)* 
distinct deals of bridge such that each player receives an ace. Therefore the 
probability of each player receiving an ace is just 
F 481. .(13. 2 2403) 
Cite Ss 32 7or B49 
Remark. Most of the above formulas contain the quantity 


R= mGE— 1) 2= 1, 


~~ 0.105. 


called n factorial. For large n, it can be shown that! 


nin me awe. 


This simple asymptotic representation of n! is known as Stirling’s formula.“ 


PROBLEMS 


1. A four-volume work is placed in random order on a bookshelf. What is the 
probability of the volumes being in proper order from left to right or from 
right to left? 


‘© The symbol ~ between two variables «,, and 8, means that the ratio «,/8, > 1 as 
n— oo, 

* Proved, for example, in D. V. Widder, Advanced Calculus, second edition, Prentice- 
Hall, Inc., Englewood Cliffs, N.J. (1961), p. 386. 


PROBLEMS BASIC CONCEPTS ll 


2. A wooden cube with painted faces is sawed up into 1000 little cubes, all of the 
same size. The little cubes are then mixed up, and one is chosen at random 
What is the probability of its having just 2 painted faces? 

Ans. 0.096. 


3. A batch of m manufactured items contains k defective items. Suppose m 
items are selected at random from the batch. What is the probability that / of 
these items are defective? 


4. Ten books are placed in random order on a bookshelf. Find the probability 
of three given books being side by side: 
Ans. 7x. 


5. One marksman has an 80% probability of hitting a target, while another has 

only a 70% probability of hitting the target. What is the probability of the 

target being hit (at least once) if both marksman fire at it simultaneously ? 
Ans. 0.94. - 


6. Suppose 7 people sit down at random and independently of each other in an 
auditorium containing n + k seats. What is the probability that m seats specified 
in advance (m < n) will be occupied? 


7. Three cards are drawn at random from a full deck. What is the probability 
of getting a three, a seven and an ace? 


8. What is the probability of being able to form a triangle from three segments 
chosen at random from five line segments of lengths 1, 3, 5, 7 and 9? 

Hint. A triangle cannot be formed if one segment is longer than the sum of 
the other two. 


9. Suppose a number from 1 to 1000 is selected at random. What is the proba- 
bility that the last two digits of its cube are both 1? 

Hint There is no need to look through a table of cubes. 

Ans. 9.01. 


10. Find the probability that a randomly selected positive integer will give a 
number ending in 1 if it is 

a) Squared; 

b) Raised to the fourth power; 

c) Multiplied by an arbitrary positive integer. 

Hint. It is enough to consider one-digit numbers. 

Ans. a) 0.2; b) 0.4; c) 0.04. 


11. One of the numbers 2, 4, 6, 7, 8, 11, 12 and 13 is chosen at random as the 
numerator of a fraction, and then one of the remaining numbers is chosen at 
random as the denominator of the fraction. What is the probability of the 
fraction being in lowest terms? 

Ans. 3%. 


12. The word “‘drawer’’ is spelled with six scrabble tiles. The tiles are then 
randomly rearranged. What is the probability of the rearranged tiles spelling 
the word “reward ?”’ 

ANS. 5635- 
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13. In throwing 6n dice, what is the probability of getting each face n times? 
Use Stirling’s formula to estimate this probability for large n. 


14. A full deck of cards is divided in half at random. Use Stirling’s formula to 

estimate the probability that each half contains the same number of red and 

black cards. 

CCis 2 
Cis V26r 

15. Use Stirling’s formula to estimate the probability that all 50 states are 

represented in a committee of 50 senators chosen at random. 


Ans ew 0.22. 


16. Suppose 2n customers stand in line at a box office, m with 5-dollar bills and 
n with 10-dollar bills. Suppose each ticket costs 5 dollars, and the box office 
has no money initially. What is the probability that none of the customers has 
to wait for change??? 


17. Prove that 


3 


(Ce)? = C=". 
k=0 


Hint. Use the binomial theorem to calculate the coefficient of x” in the 
product (1 + x)"(1 + x)" =(1 + x). 


2 A detailed solution is given in B. V. Gnedenko, The Theory of Probability, fourth 
edition (translated by B. D. Seckler), Chelsea Publishing Co., New York (1967), p. 43. 


Ans oo se Citta = —— 
: co n+1- 
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COMBINATION OF EVENTS 


3. Elementary Events. The Sample Space 


The mutually exclusive outcomes of a random experiment (like throwing 
a pair of dice) will be called elementary events (or sample points), and a 
typical elementary event will be denoted by the Greek letter w. The set of 
all elementary events w associated with a given experiment will be called the 
sample space (or space of elementary events), denoted by the Greek letter 
Q. An event A is said to be “‘associated with the elementary events of ”’ if, 
given any w in Q, we can always decide whether or not leads to the occur- 
rence of A. The same symbol A will be used to denote both the event A and 
the set of elementary events leading to the occurrence of A. Clearly, an event 
A occurs if and only if one of the elementary events w in the set A occurs. 
Thus, instead of talking about the occurrence of the original event A, we can 
just as well talk about the “occurrence of an elementary event w in the set 
A.” From now on, we will not distinguish between an event associated with 
a given experiment and the corresponding set of elementary events, it being 
understood that all our events are of the type described by saying “‘one of the 
elementary events in the set A occurs.” With this interpretation, events are 
nothing more or less than subsets of some underlying sample space 2. Thus 
the certain (or sure) event, which always occurs regardless of the outcome 
of the experiment, is formally identical with the whole space {2, while the 
impossible event is just the empty set @ , containing none of the elementary 
events w. 

Given two events A, and A,, suppose A, occurs if and only if A, occurs. 
Then A, and A, are said to be identical (or equivalent), and we write A; = Ag. 


13 
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Example 1. In throwing a pair of dice, let A, be the event that “the 
total number of spots is even’’ and A, the event that “*both dice turn up even 
or both dice turn up odd.”? Then A; = Ag. 


Example 2. In throwing three dice, let A, again be the event that “the 
total number of spots is even”’ and A, the event that “‘all three dice have either 
an even number of spots or an odd number of spots.” Then A, = Ag. 


Two events A, and A, are said to be mutually exclusive or incompatible if 
the occurrence of one event precludes the occurrence of the other, i.e., if A; 
and A, cannot occur simultaneously. 

By the union of two events A, and A;, denoted by A, U Az, we mean the 
event consisting of the occurrence of at least one of the events A, and A). 
The union of several events A,, A.,... is defined in the same way, and is 
denoted by LU A,. 


k 
By the intersection of two events A, and A,, denoted by A, M A; or simply 
by A,A., we mean the event consisting of the occurrence of both events A, 
and A,. By the intersection of several events A,, Ap,... , denoted by f) 4,, 
k 


we mean the event consisting of the occurrence of all the events A,, Ag,.. .. 

Given two events A, and Ag, by the difference A, — A, we mean the event 
in which A, occurs but not A,. By the complementary event of an event A,* 
denoted by 4, we mean the event “A does not occur.”’ Clearly, 


A=Q— A. 
Example 3. In throwing a pair of dice, let A be the event that “the 
total number of spots is even,” A, the event that “both dice turn up even,” 


and A, the event that “‘both dice turn up odd.”” Then A, and A, are mutually 
exclusive, and clearly 


A= A, VU Az, A, = A — Ag, A, = A — Aj. 


Let A, A, and A, be the events complementary to A, A, and Ag, respectively. 
Then A is the event that “the total number of spots is odd,” A, the event that 
“at least one die turns up odd,” and A, the event that ‘“‘at least one die turns 
up even.”’ It is easy to see that 


Ay = ASA, TIA Sas A, — A= Ay ti Aas 


The meaning of concepts like the union of two events, the intersection 
of two events, etc., is particularly clear if we think of events as sets of ele- 
mentary events w, in the way described above. With this interpretation, 


1 To “turn up even’’ means to show an even number of spots, and similarly for to “turn 
up odd.” 
2 Synonymously, the ““complement of A” or the “event complementary to A.” 
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given events A,, A, and A, A, U A, is the union of the sets A, and Ag, 
A, ( A, is the intersection of the sets Ay and A,, A = Q — A is the comple- 
ment of the set A relative to the whole space Q, and so on. Thus the symbols 
U, , etc. have their customary set-theoretic meaning. Moreover, the 
statement that “‘the occurrence of the event A, implies that of the event A,” 
(or simply, “A, implies A,”’) means that A, © Ag, i.e., that the set A, is a 
subset of the set A,.? 


LU] 


Y 


4 
Z 
Z 


FiGuURE 2. (a) The events A, and A, are mutually exclusive; 
(b) The unshaded figure represents the union A; U Ae; (c) The 
unshaded figure represents the intersection A; A,; (d) The 
unshaded figure represents the difference A, — A,; (e) The shaded 
and unshaded events (A, and A.) are complements of each other; 
(f) Event A, implies event A,. 


To visualize relations between events, it is convenient to represent the 
sample space {2 schematically by some plane region and the elementary 
events w by points in this region. Then events, 1.e., sets of points w, become 
various plane figures. Thus Figure 2 shows various relations between two 
events A, and A,, represented by circular disks lying inside a rectangle Q, 
schematically representing the whole sample space. In turn, this way of 
representing events in terms of plane figures can be used to deduce general 
relations between events, e.g., 


avi 4, = -4.. then A, > A,; 
c) If A === A, re Ag, then A SSE Ay U A. 


* The symbol © means “‘is a subset of” or “‘is contained in,’’ while > means “contains.” 
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Quite generally, given a relation between various events, we can get an 
equivalent relation by changing events to their complements and the symbols 
A, Uand ¢ to N, U and > (the sign = is left alone). 


Example 4. The following relations are equivalent: 
=< 45 k 
NM A, = B> M Cy 
k k 
k k 


Remark. It will henceforth be assumed that all events under considera- 
tion have well-defined probabilities. Moreover, it will be assumed that all 
events obtained from a given sequence of events Aj, Ag, ... by taking unions, 
intersections, differences and complements also have well-defined probabilities. 


4. The Addition Law for Probabilities 


Consider two mutually exclusive events A, and A, associated with the 
outcomes of some random experiment, and let A = A, U A, be the union of 
the two events. Suppose we repeat the experiment a large number of times, 
thereby producing a whole series of ‘“‘independent trials under identical 
conditions.’ Let m be the total number of trials, and let n(A,), n(A,) and 
n(A) be the numbers of trials leading to the events A,, A, and A, respectively. 
If A occurs in a trial, then either A, occurs or A, occurs, but not both (since 
A, and A, are mutually exclusive). Therefore 


: n(A) = n(A,) + n(Ag), 
and hence 


m(A) _ (As) n(Aa) 
n n n 


But for sufficiently large n, the relative frequencies n(A)/n, n(A,)/n and 
n(A,)/n virtually coincide with the corresponding probabilities P(A), P(A,) 
and P(A,), as discussed on p. 3. It follows that 


P(A) = P(A,) + P(A,). (2.1) 


Similarly, if the events A,, A, and A, are mutually exclusive, then so are 
A, U A, and Az, and hence, by two applications of (2.1), 


P(A, U A, VU As) = P(A, UV Ag) = P(As) = P(A,) + P(A) + P(As3). 
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More generally, given n mutually exclusive events A,, A,,...,A,, we have 
the formula 
p(U A, = Y P(A), (2.2) 
~ k=1 
obtained by applying (2.1) n — 1 times. Equation (2.2) is called the addition 


law for probabilities. 
Next we prove some key relations involving probabilities: 


THEOREM 2.1. The formulas 


0< P(A) <1, (2.3) 
P(A, — A,) = P(A)) — P(A, 1 Ay), (2.4) 
P(A, — Aj) = P(A,) — P(A, +1 A,), (2.5) 
P(A, U Az) = P(A,) + P(Ag) — P(A, 19 Ag) (2.6) 
hold for arbitrary events A, A, and A,. Moreover, 
P(A,) < P(Ag) if A, < Ap. CAs 


Proof. Formula (2.3) follows at once from the interpretation of 
probability as the limiting value of relative frequency, since obviously 


pe 
n 


where n(A) is the number of occurrences of an event A inv trials.4 Given 
any two events A, and A,, we have 


A, = (A, — Ay) U (Ay 9 AQ), 
A, = (A, — Ay) U (A, 9 AQ), 
A, U A, = (Ay — Ag) U (Az — Ay) U (Ad 29 Ad), 
where the events A, — A,, A, — A, and A, AM A, are mutually exclusive. 
Therefore, by (2.2), 
P(A;) = P(A; — Ag) + P(A 21 Az), (2.8) 
P(A) = P(A, — A) + P(4, 9 Ag), (2.9) 
P(A, U Az) = P(A; — Ag) + P(4g — Ay) + P(41 9 AQ). (2.10) 
Formulas (2.8) and (2.9) are equivalent to (2.4) and (2.5). Then, using 
(2.4) and (2.5), we can write (2.10) in the form 
P(A, U A,) = P(A,) — P(4, 0 A,) - P(A.) 
—P(A, N Ag) + P(A, 21 AQ) 
= P(A,) + P(Ag) — P(A1 1 AQ), 


4 Note that P(@) = 0, P(Q) = 1, since n(@) = 0, n(Q) = n for all n. Thus the impos- 
sible event has probability zero, while the certain event has probability one. 
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thereby proving (2.6). Finally to prove (2.7), we note that if A, < Ag, 
then A, M A, = A, and hence (2.9) implies 


P(A,) = P(A.) — P(A, — Ay) < P(AQ), 
since P(A, — A,) > 0 by (2.3). J 


The addition law (2.2) becomes much more complicated if we drop the 
requirement that the events be mutually exclusive: 


THEOREM 2.2 Given any n events A,, Ag, ..., An; let? 
P, a > P(A,), 
71 


P, = > P(A;A;) 


1l<i<j<xn 
PP. SSS > P(A,A,A,), Ee 
1l<i<j<k<n 
Then 
p(U.4,) = Pi — Pat Pe— Peto t Pe (2.11) 
k= 


Proof. For n = 2, (2.11) reduces to formula (2.6), which we have 
already proved. Suppose (2.11) holds for any n — 1 events. Then 


p(U4,) -> P(A)— > P(A,A,) 


2<i<j<n 


+. 22 P(A,A,;A;) — = (2.12) 
I<i<j<k<n 
and 


P( U4,A,) = ¥ PAM) 7S * P44 AD 
k=2 i=2 t<t<j<n 


2<i<j<k<n 


But, by (2.6), 
p(U 4,) = P(A,) + p(U 4.) = p(U 4,4x). 


and hence, by (2.12) and (2.13), 


p(U 4.) =P(A) + 3P(4)—_ 3S _ PCAAY 


I< 95 
n 
+ >> PALA A, So = = P(A,A,) 
2<i<j<k<n i=2 
= 3 P(A, A,A,) — °° = £1 9 
2<i<j<n 


5 4,A; is shorthand for the intersection A; \ A;, A,;A;A; is shorthand for A; 0 A; Az, 
and so on. Inasum like = P(A;A;A,), each group of indices (satisfying the indicated 
ae 1<i<j<k<n 
inequalities) is encountered just once. 
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i.e., (2.11) holds for any n events. The proof for all n now follows by 
mathematical induction. J 


Example (Coincidences). Suppose n students have n identical raincoats 
which they unwittingly hang on the same coat rack while attending class. 
After class, each student selects a raincoat at random, being unable to tell it 
apart from all the others. What is the probability that at least one raincoat 
ends up with its original owner? 


Solution. We number both the students and the raincoats from 1 to n, 
with the Ath raincoat belonging to the kth student (k = 1,2,...,n). Let 
A,, be the event that the Ath student retrieves his own raincoat. Then the 
event A that “at least one raincoat ends up with its original owner’ is just 


Every outcome of the experiment consisting of “‘randomly selecting’’ the 
raincoats can be described by a permutation (4, ij,...,i,), where i, is the 
number of the raincoat selected by the Ath student. Consider the event 
A; A; °° * Ax,» Where m <n. This event occurs whenever i, = ky, i, = 
ky, ...5%,,, =k», and the other indices take the remaining n — m values in 
any order. Therefore 
tia A, A, ) (n — m)! 
P(A, Ax, oe A, ) = sea TT = ae tlt , 

where N(A,, A,,*°* A;,,) = (2 — m)! is just the total number of permuta- 
tions of n — m things, and N = n! is the total number of permutations of 
n things (m is the number of fixed indices k,, ky,...,k,,). There are pre- 
cisely 


. n! 


= m!(n — m)! 
distinct events of the type A, A;, °° A,,,, with m fixed indices, this being the 
number of combinations of n things taken m at a time (recall Theorem 1.3, 
p. 7). It follows that 
. a hm)! 
Pu= XS _PlAg Ang" ** Aig) = Ch SE = FS 


1 
1<kyj<kg<+++<km< n! m! 


Hence, by formula (2.11), 


p(U4,) | Po 
k=1 
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i.e., the desired probability P(A) is a partial sum of the power series expansion 
of the function | — e* with x = —1: 


Thus, for large n, 
P(A) = 1 — e 1 & 0.632. (2.14) 


To generalize the addition law to the case of an infinite sequence of mutu- 
ally exclusive events A,, Ao,..., we repeatedly apply (2.1). Thus 


P(A, U A, U A, U-+-)= P(A,) + P(A U Ag Oe 
= P(A,) + P(A,) + P(As U--*) 
= P(A,) + P(A,) + P(A3) +°--, 
or equivalently, 
p(U Ay) = SA). 
k=1 k=1 
We can combine this formula and (2.2) into a single formula 


P(U A, — 2 P(A,), (2.2’) 


where it will always be clear from the context whether UJ and > have finite 
or infinite limits.® : . 

The “generalized addition law” (2.2’) has a number of important con- 
sequences. We begin with two theorems expressing a kind of “continuity 
property” of probability: 


THEOREM. 2.3. If Ay, Ao, ... is an “‘increasing sequence” of events, 
i.e., a sequence such that A, < A, < +::, then 
p(U A, — lim P(A,). (2.15) 
k nN 00 


Proof. Clearly, the events 


n—1 


B,=—A, B,=A,—A,,..., 3B, =A, — ee 
k=1 


® In the last analysis, formulas (2.1), (2.2’) and (2.3) are axioms, although they are, of 
course, strongly suggested by experience, i.e., by the interpretation of probabilities as 
limiting values of relative frequencies. In this sense, they are the “‘only reasonable axioms.” 
and lead to a model of random phenomena whose consequences are fully confirmed by 
experiment. 
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are mutually exclusive and have union U A,. Moreover, 
k 


Therefore, by (2.2’), 
P(U Ay) = P(U B,) = 2 PB) = lim SPB) 


—lim P (U B,) —limP(A,). I 
N— co pa n> oo 


Similarly, we have 


THEOREM 2.3’. If A,, As,... is a “decreasing sequence’ of events, 
i.e., a sequence such that A, > A, >°--, then : 


no 


P(f) A,) =lim P(A,). 
(r ) (A,) 


Proof. Going over to complementary events, we have 4, © 4, ¢-:-, 
and hence, by (2.15), 


P(N A) = P(U 4, = 1—limP(A,) 


n> 0 


—lim [1 — P(4,)] =lim P(A,). I 


In the case of arbitrary events, we must replace = by < in (2.2’): 


THEOREM 2.4. The ineqaulity 
P(U Ay) < ¥ P(A,) 
k k 
holds for arbitrary events A,, Ag,.... 


Proof. As in the proof of Theorem 2.3, U A, is the union of the 
k 


mutually exclusive events (2.16), where obviously B, < A, and hence 
P(B,,) < P(A,), by (2.7). Therefore 


p(U4,) = P(U By) = PB) < DPA). | 

k k k k 

Finally, we prove a proposition that will be needed in Chapter 7: 
THEOREM 2.5 (First Borel-Cantellilemma).’ Given a sequence of events 


’ For the “second Borel-Cantelli lemma,’ see Theorem 3.1, p. 33. 
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Ay, Ag,..., with probabilities p, = P(A,), kK =1,2,..., suppose 


20 (249) 


i.e., suppose the series on the left converges. Then, with probability | 
only finitely many of the events A,, Ag, ..., occur. 


Proof. Let B be the event that infinitely many of the events 4,, 
As, .... occur, and let 
i Ay 
k>n 
so that B,, is the event that at least one of the events A,, Aji,. 


occurs. Clearly B occurs if and only if B, occurs for every n = 1, 
2.4.3 1 eter 


B=1B, =f) (UA). 


n k>n 
Moreover, B, > B, > :::, and hence, by Theorem 2.3’, 
P(B) = lim P(B,). 


But, by Theorem 2.4, 
PB) = > PG) => f, — 0 as n— o, 
k>n k>n 
because of (2.17). Therefore 
PS) =e) = 0, 
i.e., the probability of infinitely many of the events A, Az, ... occurring 


is 0. Equivalently, the probability of only finitely many of the events 
A,, Ae;.... OCcurTingis 1; ff 


PROBLEMS 


1. Interpret the following relations involving events A, B and C: 
a) AB=A; b) ABC=A; c) AVUBUC=A., 
2. When do the following relations involving the events A and B hold: 


a) AVB=A; b) AB=A; ¢) AUB= AB? 


3. Simplify the following expressions involving events A, B and C: 
a) (A UB\(BUC); b) (A UB)(A UB); c) (A UB)(A VU B)A UV B). 
Ans. a) ACU B; b) A; c) AB. 


4. Given two events A and B, find the event X such that 


Yu AUC A) =e 
Ans. X = B. 
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5. Let A be the event that at least one of three inspected items is defective, and B 
the event that all three items are of acceptable quality. What are the events 
A U Band AB? 


6. A whole number from 1 to 1000 is chosen at random. Let A be the event that 
the number is divisible by 5, and B the event that the number ends in a zero. 
What is the event AB? 


7. A target is made up of 10 circular disks bounded by 10 concentric circles of 
PadiT,.fo,-.-.1o wherer, <r, < -* *-< to. Let.A, be the event consisting of 
the disk of radius r;, being hit (k = 1, 2,..., 10). What are the events 


6 10 
B= U A; C= f) Ax? 
k=1 k=5 
Ans. B= Ag, C = Az. 
8. Given any event A, prove that 


P(A) = ea Ps 8 PLO. 
9. A marksman fires at a target made up of a central circular disk and two 


concentric rings. The probabilities of hitting the disk and the rings are 0.35, 
0.30 and 0.25, respectively. What is the probability of missing the target? 


10. Five items are chosen at random from a batch of 100 items and then in- 
spected. The whole batch is rejected if any of the items is found to be defective. 
What is the probability of the batch being rejected if it contains 5 defective items? 


95-94 99> 92°91 


fon ee20e oF 66 


Ans. 


11. A secretary forgets the last digit of a telephone number, and dials the last 
digit at random. What is the probability of calling no more than three wrong 
numbers? How is this probability changed if she recalls that the last digit is 
even? 


12. Given any n events Aj, Ap,..., An, prove that 


P( 14%] 5 P( U 4) 


13. A batch of 100 manufactured items contains 5 defective items. Fifty items 
are chosen at random and then inspected. Suppose the whole batch is accepted if 
no more than one of the 50 inspected items is defective. What is the probability 
of accepting the whole batch? 


14. Write an expression for the probability p(r) that among r randomly selected 
people, at least two have a common birthday. 


Comment. Rather surprisingly, it turns out that p(r) > 4 if r = pe 


®"See W. Feller, op. cit., p. 33. 
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15. Test the approximation (2.14) for n = 3, 4, 5 and 6. 


16. Use Theorem 2.2 and Stirling’s formula to find the probability that some 
player is dealt a complete suit in a game of bridge. 


16 72 72 339 397 


As. = oS ————— x» — | wx}i ae 

" Ch Cen CHCRCR” [gy = 
17. Given any n events A;, Ag,..., A,, prove that the probability of exactly 
m (m <n) of the events occurring is 


m+1 m-+2 n 
in 8 Pm + a Pg oe in Ps 


where P,,, Pmii,--. are the same as in Theorem 2.2. 


18. Let n = 10 in the example on p. 19. What is the probability that exactly 
5 raincoats end up with their original owners ? 


19. A whole number from 1 to 1000 is chosen at random. What is the proba- 
bility of its being a power (higher than the first) of another whole number? 
Hint. 31? < 1000 < 322. 
Ans. 35. 
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5. Conditional Probability 


In observing the outcomes of a random experiment, one is often inter- 
ested in how the outcome of one event A is influenced by that of another 
event B. For example, in one extreme case the relation between A and B may 
be such that A always occurs if B does, while in the other extreme case A never 
occurs if B does. To characterize the relation between A and B, we introduce 
the conditional probability of A on the hypothesis B, i.e., the “probability of 
A occurring under the condition that B is known to have occurred.” This 
quantity is defined by 


ri (3.1) 


where AB is the intersection of the events A and B, and it is assumed that 
P(B) > 0. 

To clarify the meaning of (3.1), consider an experiment with a finite 
number of equiprobable outcomes (elementary events). Let N be the total 
number of outcomes, N(B) the number of outcomes leading to the occurrence 
of the event B, and N(AB) the number of outcomes leading to the occurrence 
of both A and B. Then, as on p. 1, the probabilities of B and AB are just 


pe) =X, pap) = NAB. a) 
and hence (3.1) implies 
__ N(AB) 
P(A | B) = hay (3.3) 


25 


26 DEPENDENT EVENTS CHAP. 3 


But (3.3) is of the same form as (3.2), if we restrict the set of possible out- 
comes to those in which B is known to have occurred. In fact, the denomin- 
ator in (3.3) is the total number of such outcomes, while the numerator is the 
total number of such outcomes leading to the occurrence of A. 

It is easy to see that conditional probabilities have properties analogous 
to those of ordinary probabilities. For example, 


a) 0 < P(A| B) < 

b) If A and B are lacus so that AB = @, then P(A | B) = 0; 

c) If Biriplies A, so that B < A, then P(A | B) = 1; 

d) If A,, A,,... are mutually exclusive events wie union A = J A,, 
then k 


P(A | B) = > P(A, | B) (3.4) 
k 
(the addition law for conditional probabilities). 


Property a) is an immediate consequence of (3.1) and the formula 0 < 
P(AB) < P(B), implied by @ © ABC B. To prove b), we note that 
AB = © implies P(AB) = 0 and hence P(A | B) = 0, by (3.1). Similarly, 
c) follows from the observation that if B < A, then AB = B, P(AB) = P(B), 
and hence P(A | B) = 1, by (3.1). Finally, if 4 = U A,, where A,, As,... 
are mutually exclusive events, then 


AB = U A,B, 
k 
and hence 


P(AB) = 5 P(4,B), (3.5) 


by formula (2.2'), p. 20, the addition law for ordinary probabilities. Dividing 
(3.5) by P(B), we get (3.4), because of (3.1) and 


P(A,B) 


P(A, | B) — P(B) 


In calculating the probability of an event A, it is often convenient to use 
conditional probabilities as an intermediate step. Suppose B,, Bz,... is a 
“full set’’* of mutually exclusive events, in the sense that one (and only one) 
of the events B,, B,,... always occurs. Then we can find P(A) by using the 
“total probability formula”’ 


P(A) = > P(A | B,)P(B,). (3.6) 


* Synonymously, an “exhaustive set.” 
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To prove (3.6), we need only note that 
U B, = Q, 
k 
where {2 is the whole sample space, since one of the events B,, B,,... must 
occur. But then 
A =f AB,, 
and hence : 


P(A) = P/ U AB,) = 5 P(4B,) = > PAB) pvp), 


z P(B;) 
which is equivalent to (3.6). 


Example 1. A hiker leaves the point O shown in Figure 3, choosing one 
of the roads OB,, OB,, OB;, OB, at random. At each subsequent crossroads 
he again chooses a road at random. What is the probability of the hiker 
arriving at the point A? 


FIGURE 3 


Solution. Let the event that the hiker passes through the point B,, k = 
1,...,4, be denoted by the same symbol B, as the point itself. Then B,, By, 
B;, B, form a “full set’? of mutually exclusive events, since the hiker must 
pass through one of these points. Moreover, the events B,, Bz, Bs, By are 
equiprobable, since, by hypothesis, the hiker initially makes a completely 
random choice of one of the roads OB,, OB,, OB;,, OB,. Therefore 


] 
P(B,) = -, Fes liga. ne 
4 
Once having arrived at B,, the hiker can proceed to A only by making the 


proper choice of one of three equiprobable roads. Hence the conditional 
probability of arriving at A starting from B, is just 3. Let the event that the 
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hiker arrives at A be denoted by the same symbol A as the point itself. Then 


l 
P(A | By) = 3. 


and similarly 


2 
P(A|B) =>, P(A|B)=1, P(A|B) =; 


(consult the figure). It follows from (3.6) that the probability of arriving at 
A is 
P(A) = P(A | B,)P(B;) + P(A | B,)P(Be) 
+ P(A | B,)P(Bs) + P(A | By)P(B,) 


67 
120° 


“(; 1 , 4 
"43-4, 32 

Example 2 (The optimal choice problem). Consider a set of m objects, 
all of different quality, such that it is always possible to tell which of a given 
pair of objects is better. Suppose the objects are presented one at a time and 
at random to an observer, who at each stage either selects the object, thereby 
designating it as “the best’? and examining no more objects, or rejects the 
object once and for all and examines another one. (Of course, the observer 
may very well make the mistake of rejecting the best object in the vain hope 
of finding a better one!) For example, the observer may be a fussy young lady 
and the objects a succession of m suitors. At each stage, she can either accept 
the suitor’s proposal of marriage, thereby terminating the process of selecting 
a husband, or she may reject him (thereby losing him forever) and wait for a 
better prospect to come along. It will further be assumed that the observer 
adopts the following natural rule for selecting the best object: ““Never select 
an object inferior to those previously rejected.’’ Then the observer can select 
the first object and stop looking for a better one, or he can reject the first 
object and examine further objects one at a time until he finds one better 
than those previously examined. He can then select this object, thereby 
terminating the inspection process, or he can examine further objects in the 
hope of eventually finding a still better one, and so on. Of course, it is 
entirely possible that he will reject the very best object somewhere along the 
line, and hence never be able to make a selection at all. On the other hand, 
if the number of objects is large, almost anyone would reject the first object 
in the hope of eventually finding a better one. 

Now suppose the observer, following the above “decision rule,”’ selects 
the ith inspected object once and for all, giving up further inspection. (The 
ith object must then be better than the i — 1 previously inspected objects.) 
What is the probability that this ith object is actually the best of all m objects, 
both inspected and uninspected ? 
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Solution. Let B be the event that the last of the i inspected objects is the 
best of those inspected, and let A be the event that the ith object is the best 
of all m objects, both inspected and uninspected. Then we want the condi- 
tional probability P(A | B) of the event A given that B has already occurred. 
According to (3.1), to calculate P(A | B) we need both P(B) and P(AB). 
Obviously A < Band hence AB = A, so that P(AB) = P(A). By hypothesis, 
all possible arrangements of the objects in order of presentation are equi- 
probable (the objects are presented ‘“‘at random”). Hence P(B) is the proba- 
bility that in a random permutation of i distinguishable objects (the objects 
differ in quality) a given object (the best of all i objects) occupies the ith 
place. Since there are i! permutations of all i objects and (i — 1)! permuta- 
tions subject to the condition that a given object occupy the ith place, this 
probability is just 


pes 2. 
l 


i! 


Similarly, P(A) is the probability that in a random permutation of m dis- 
tinguishable objects, a given object (the best of all m objects) occupies the 
ith place, and hence 


meaty) 
P(A) = (m#:— 1)! ass 
m! m 
Therefore the desired conditional probability P(A | B) is just 


P(AB PA) 4 
P(B) P(B) mm 


Example 3 (The gambler’s ruin). Consider the game of “‘heads or tails,”’ 
in which a coin is tossed and a player wins 1 dollar, say, if he successfully 
calls the side of the coin which lands upward, but otherwise loses 1 dollar. 
Suppose the player’s initial capital is x dollars, and he intends to play until 
he wins m dollars but no longer. In other words, suppose the game continues 
until the player either wins the amount of m dollars, stipulated in advance, 
or else loses all his capital and is “‘ruined.’’ What is the probability that the 
player will be ruined? 


Solution. The probability of ruin clearly depends on both the initial 
capital x and the final amount m. Let p(x) be the probability of the player’s 
being ruined if he starts with a capital of x dollars. Then the probability of 
ruin given that the player wins the first call is just p(x + 1), since the player’s 
capital becomes x + 1 if he wins the first call. Similarly, the probability of 
ruin given that the player loses the first call is p(x — 1), since the player’ S 
capital becomes x — 1 if he loses the first call. In other words, if B, is the 
event that the player wins the first call and B, the event that he loses the first 
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call, while A is the event of ruin, then 
P(A | B,) = p(x + 1), P(A | B,) = p(x — 1). 


The mutually exclusive events B, and B, form a “full set,” since the player 
either wins or loses the first call. Moreover, we have 


] l 
P(B) = 5: P(B,) =>» 


assuming fair tosses of an unbiased coin (cf. Problem 1, p. 65). Hence, by 
(3.6), 

: P(A) = P(A | B,)P(B,) + P(A| B,)P(B,), 

1.€., 


po) =F [p+ D+p——-, l<x<m—1, G7) 


where obviously 
p(0o)=1, p(m)=90. = (3.8) 


The solution of (3.7) is the linear function 
p(x) = Cy + Cox, (3.9) 


where the coefficients C, and C, are determined by the boundary conditions 
(3.8), which imply 
ce eS t c + Com = 0. (3.10) 


Combining (3.9) and (3.10), we finally find that the probability of ruin given 
an initial capital of x dollars is just 


. px) =1—=, C= xe 
m 


6. Statistical Independence 


In saying that two experiments are “‘statistically independent” (or briefly, 
‘‘independent’’), we mean, roughly speaking, that the outcome of one experi- 
ment has no influence on the outcome of the other. Let A, be an event 
associated only with the first experiment, and A, an event associated only with 
the second experiment. Then the occurrence of A; has no influence on the 
probability of occurrence of Ag, and conversely. In this sense, we say that the 
events A, and A, are “‘(statistically) independent.” 

To give mathematical expression to the notion of independence, we 
calculate the probability that two independent events A, and A, both occur. 
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To this end, we again resort to the empirical fact that the relative frequency 
of an event in a large series of “‘independent trials under identical conditions’? 
virtually coincides with its probability (recall Sec. 1). Imagine a long series 
of such trials, where each trial involves carrying out both experiments. If n 
is the total number of trials and n(A,A,) the number of trials leading to 
occurrence of both A, and Ag, then 


P(A,AQ) <a 


n(A,A2) 
ree ‘ (3.11) 


Moreover, if n(A,) is the number of trials leading to occurrence of A,, then 
A 

pass) (3.12) 
n 


Suppose we confine ourselves to examining the results of the n(A,) trials 

leading to occurrence of A., and look for occurrence of A,. Then clearly A, 

will occur in precisely n(A,A,) of these trials. Moreover, if” is very large, then 

so is n(A,), and hence 

n(A;Ag) 
n(Ag) 


since A, is associated only with the second experiment, which has nothing 
whatsoever to do with the first experiment or the event A, associated with it. 
Combining (3.11)-(3.13), we find that 


P(A,) ~ (3.13) 


n(A,A,) _ n(A,Ag) n(A;) a 
~ P(A,A2) eS oa P(A,)P(A3), 


or, after going over to exact equations (in the limit as n — oo), 


P(A,A;) = P(A,)P(A,). (3.14) 


Two events A, and A, are said to be (statistically) independent if they satisfy 
(3.14) and (statistically) dependent otherwise.° 

The definition (3.14) is in keeping with the notion of conditional proba- 
bility introduced in Sec. 5. In fact, if two events A, and A, are independent, 
then, loosely speaking, the occurrence of A, should have no influence on the 
probability of occurrence of A,, and hence the conditional probability 


2 Thus there remains the problem of just what is meant by “independent trials under 
identical conditions” (a phrase already encountered on pp. 2 and 16), although the in- 
tuitive meaning of the phrase is perfectly clear, e.g., in a series of coin tosses. For a rigorous 
discussion of this whole issue, see W. Feller, op. cit., p. 128. 

3 In the last analysis, (3.14) is a definition, although one strongly suggested by experience, 
i.e., by the intuitive meaning of independence and the interpretation of probabilities as 
limiting values of relative frequencies (recall footnote 6, p. 20). 
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P(A, | A,) of A, occurring given that A, has already occurred should be the 
same as the unconditional probability of Ay, 1.e., 


P(A, | A,) = P(A) 
(and similarly with A, and A, changing places). But clearly 


P(A, | Ag) = —— — P(A) 


if and only if (3.14) holds. 


Example 1. Let A, be the event that a card picked at random from a 
full deck is a spade, and A, the event that it is a queen. Are A, and A, 
independent events ? 


Solution. The question is not easily answered on the basis of physical 
intuition alone. However, noting that a full deck (52 cards) contains 13 
spades and 4 queens, but only one queen of spades, we see at once that 


73-4 4 1 l 
P(A Se eC P(A =--—- >—- =; P(A,A = —, 
( 7) 52 4 ? ( 2) 52 13 ( 1 2) 52 
and hence P(A,A,) = P(A,)P(A,). Therefore the events A, and A, are inde- 
pendent. 


Example 2. In throwing a pair of dice, let A, be the event that “the 
first die turns up odd,” A, the event that “‘the second die turns up odd,” and 
A, the event that “‘the total number of spots is odd.”’ Clearly, the number of 
spots on one die has nothing to do with the number of spots on the other die, 
and hence the events A, and A, are independent, with probabilities 


l 1 
P(A) = 3 ? P(A;) = 5) : 


Moreover, it is clear that 


1 

P(A;) = =. 

(4s) = 5 

Given that A, has occurred, A, can occur only if the second die turns up even. 
Hence 


1 
P(A; | A;) a 4 ’ 


and similarly 
1 
P(A, | Az) a 5) « 
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It follows that 
P(A; | A,) = P(As), P(A; | A,) = P(A3). 


Therefore the events A, and A; are independent, and so are the events A, 
and As. 


Generalizing (3.14), we have the following 


DEFINITION. The events A,, Ao,..., A, are said to be (mutually) 
independent if 
P(A,A,;) = P(A,)P(A,), 


P(A,A;A,) = P(A,)P(A;)P(A,); 


P(A, A, °** An) = P(A;)P(Az) *** P(A,) 
for all combinations of indices such that! <i<j<-+:'<k<n, 


Example 3. The events 4,, A, and A, in Example 2 are not independent, 
even though they are “pairwise independent” in the sense that 


P(4,A,) ae P(A,)P(A;) 


for all 1 < i <j < 3. In fact, A, obviously cannot occur if A, and A, both 
occur, ancnence 


P(A,A,A3) — 0. 


But 
ee wees eae 
P A P A P A ee ee oes =—, 
(A)P(A)P(As) = 5°5°5 = 5 
so that 
P(A, A343) # P(A;)P(Az)P(A3). 
Given an infinite sequence of events A;, A,,..., suppose the events 
A,,...,A,, are independent for every n. Then A;, Ag,... is said to be a 


sequence of independent events. 


THEOREM 3.1 (Second Borel-Cantelli lemma). Given a sequence of 
independent events A,, Ay,..., with probabilities p, = P(A,), k = 1, 
2,..., Suppose 


yi &, (3233) 


Co 
I 
an 


i.e., suppose the series on the left diverges. Then, with probability | 
infinitely many of the events A, Ag, ... occur. 
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Proof. As in the proof of the first Borel-Cantelli lemma (Theorem 
Zo; P. 21), ret 


k2>n n nm \k>n 


so that B occurs if and only if infinitely many of the events A,, As,... 
occur. Taking complements, we have 


k>n 
Clearly, 
n+m a 
i, > M A, 
k=n 
for every m= 0,4) 2). . ; daieretore 


P(B,) < p( 1.4) Pd) Pe) 


n+m 
= (1 Py)" (1 — Prem) < exp (—¥ Pi), G6) 


where we use the inequality 1 — x < e-*, x > Oand the fact that if A,, 
Ag,... is a sequence of independent events, then so is the sequence 
of complementary events 4,, A:,...* But 

nN+m 

> Ppp> oO a m—>o, 

k=n 
because of (3.15). Therefore, passing to the limit m— oo in (3.16), 
we find that P(B,,) = 0 for every n = 1, 2,... It follows that 
P(B) < > P(B,) = 0, 


and hence 
P(B) = 1 — PGS) =, 


i.e., the probability of infinitely many of the events A,, Ag, ... occurring 
is]. 


PROBLEMS 


1. Given any events A and B, prove that the events A, ABand A UB forma 
full set of mutually exclusive events. 


4 It is intuitively clear that if the events A,,..., A, are independent, then so are their 
complements. Concerning the rigorous proof of this fact, see Problem 7 and W. Feller, 
op. cit., pp. 126, 128. 
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2. In a game of chess, let A be the event that White wins and B the event that 
Black wins. What is the event C such that A, BandC form a full set of mutually 
exclusive events ? 


3. Prove that if P(A| B) > P(A), then P(B| A) > P(B). 
4. Prove that if P(A) = P(B) = , then P(A | B) > 3. 
5, Given any three events A, B and C, prove that 
P(ABC) = P(A)P(B| A)P(C| AB). 
Generalize this formula to the case of any 7 events. 
6. Verify that 
: P(A) = P(A| B) + P(A| B) 
. a) A= 2; bp RB=o; co) B= Q;. d) B= A; e) B= A. 
7. Prove that if the events A and B are independent, then so are their comple- 


ments. 
Hint. Clearly P(B| A) + P(B| A) = 1 for arbitrary A and B. Moreover 


P(B| A) = P(B), by hypothesis. Therefore P(B| A) = 1 — P(B) = P(¥), so 
that A and B are independent. 

8. Two events A and B with positive probabilities are incompatible. Are 
they dependent? 

9. Consider n urns, each containing w white balls and b black balls. A ball is 
drawn at random from the first urn and put into the second urn, then a ball is 


drawn at random from the second urn and put into the third urn, and so on, 
until finally a ball is drawn from the last urn and examined. What is the prob- 


ability of this ball being white? 


w 
wt+b 

10. In Example 1, p. 27, find the probability of the hiker arriving at each of 
the 6 destinations other than A. Verify that the sum of the probabilities of 
arriving at all possible destinations is 1. 


Ans. 


11. Prove that the probability of ruin in Example 3, p. 29 does not change if 
the stakes are changed. 

12. Prove that the events A and B are independent if P(B| A) = P(B| A). 

13. One urn contains w, white balls and b, black balls, while another urn 
contains w, white balls and b, black balls. A ball is drawn at random from each 
urn, and then one of the two balls so obtained is chosen at random. What is the 
probability of this ball being white? 


? 1 Wy . Wo 
ea ye w, + by We + de : 
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14. Nine out of 10 urns contain 2 white balls and 2 black balls each, while the 
other urn contains 5 white balls and 1 black ball. A ball drawn from a randomly 
chosen urn turns out to be white. What is the probability that the ball came from 
the urn containing 5 white balls? 

Hint. If By,..., B, is a full set of mutually exclusive events, then 


Py] 4) < PBOPAL BY) __P@IPCA| Be) 


P(A a3 
oS PGaeaae 
k=1 
a formula known as Bayes’ rule. The events B,,..., B, are often regarded as 


“hypotheses” accounting for the occurrence of A. 
Ans. 5. 


15. One urn contains only white balls, while another urn contains 30 white 
balls and 10 black balls. An urn is selected at random, and then a ball is drawn 
(at random) from the urn. The ball turns out to be white, and is then put back 
into the urn. What is the probability that another ball drawn from the same urn 
will be black? 


3_ 
Ans. 28° 


16. Two balls are drawn from an urn containing n balls numbered from 1 to n. 
The first ball is kept if it is numbered 1, and returned to the urn otherwise. What 
is the probability of the second ball being numbered 2? 


n—n+1 
n2(n — 1) — 


17. A regular tetrahedron is made into an unbiased die, by labelling the four 
faces a, b, c and abc, respectively. Let A be the event that the die falls on either 
of the two faces bearing the letter a, B the event that it falls on either of the two 
faces bearing the letter b, and C the event that it falls on either of the two faces 
bearing the letter c. Prove that the events A, Band Care “pairwise independent”® 
but not independent. 


Ans. 


18. An urn contains w white balls, b black balls and r red balls. Find the 
probability of a white ball being drawn before a black ball if 


a) Each ball is replaced after being drawn; 
b) No balls are replaced. 


Ans. in both cases. 


w 
w+b 


® As defined in Example 3, p. 33. 
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RANDOM VARIABLES 


7. Discrete and Continuous Random Variables. 
Distribution Functions 


Given a sample space 2, by a random variable we mean a numerical 
function & = &(w) whose value depends on the elementary events w € &. 
Let P{x’ < & < x”} be the probability of the event {x’ < § < x"}, Le., the 
probability that € takes a value in the interval x’ < x < x". Then knowledge 
of P{x’ < & < x”} for all x’ and x” (x’ < x”) is said to specify the proba- 
bility distribution of the random variable ¢. 

A random variable & = &(w) is said to be discrete (or to have a discrete 
distribution) if & takes only a finite or countably infinite number of distinct 
values x, with corresponding probabilities 


P;(x) 555 P{E —_ Rts 


> P(x) = 1, 


where the summation is over all the values of x taken by €. For such random 
variables, | 


Pix’ < £<x"}=SPAOo, (4.1) 


where the summation is over the finite or countably infinite number of values 
of x which & can take in the interval x’ < x < x”. 

A random variable & = &(w) is said to be continuous (or to have a 
continuous distribution) if 


Pix’ <&<x"}=[" px) dx, (4.2) 
37 
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where p,(x) is a nonnegative integrable function, called the probability 
density of the random variable €, with unit integral 


FE pel) ax = 1, 
Clearly, if € is a continuous random variable, then 


P{E = x} =0 
for any given value x, while? 


P{E € dx} ~ p,(x) dx 


for every x with a neighborhood in which the probability density p,(x) is 
continuous. Here P{& € dx} is the probability of the event {@ € dx}, con- 
sisting of € taking any value in an infinitesimal interval dx centered at the 
point x. 

The function 


O(xy= PE < xj; > —- eer Sw 


is called the distribution function of the random variable &. If & is a discrete 
random variable, ®,(x) is the step function 


(2) = > pl 


taking a finite or countably infinite number of distinct values [the graph of 
such a function is shown in Figure 4(a)]. If & is a continuous random 
variable, ®,(x) is the continuous function 


(x) =|" x(x) dx (4.3) 
[the graph of such a function is shown in Figure 4(b)].? Clearly, 
Pix’ < §& < x"} = Ox") — Ofx’) (4.4) 


for any random variable &. 

Now consider two random variables €, and & , or equivalently, the ran- 
dom point or vector & = (&, &). First suppose €, and &, are discrete. Then 
€, and &, have a joint probability distribution, characterized by the proba- 
bilities 

De,,t(%13 X9) a Pie, = X41, Ee aay Xo}, (4.5) 


where x, and x, range over all possible values of the corresponding random 


1 The symbol € means “‘belongs to”’ or “is contained in.” 

2? By a well-known theorem on differentiation, d® e(x)/dx = P:(x) almost everywhere. 
See e.g., E. C. Titchmarsh, The Theory of Functions, second edition, Oxford University 
Press, London (1939), p. 362. 
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0 
(b) 


Ficure4. (a) A typical distribution function of a discrete random 
variable taking only the integral values ..., —2, —1,0,1,2,... 
At the points x =..., —2, a | a Oe eae ®,.(x) has jumps 
equal to the corresponding probabilities P,(x). (b) A_ typical 
distribution function of a continuous random variable. Any 


continuous monotonic function ®,(x) such that lim ®,(x) = 0, 
r—>— 0 


lim ©,(x) = 1 can serve as the distribution function of a 
t—>-+ 
continuous random variable &.* 


variables &, and &. The probability of any event of the type {(&1, &) € B}, 
ie., the “probability of the random point & = (&, &2) falling in a given set 
B,”’ is given by 
PiGi, 3) B} = = Pé,,é(%13 Xe) 
(v1,4,)EB 

where the summation is over all possible values x,, x2 of the random variables 
€., & such that the point (x, x2) lies in B. Next suppose &, and & are con- 
tinuous. Then by the joint probability density of & and &, we mean a 


3 It should be noted that there are random variables which are neither discrete nor 
continuous but a “mixture of both.” There are also continuous distribution functions more 


general than (4.3). 
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function pz, ¢,(%1, X2) of two variables x, and x, such that the probability of 
any event of the type {(&, &) € B} is given by 


P(E, £2) € B} =| | Pg, g4(Xas X2) dx, dp (4.6) 
B 


(the integral is over B). 

Given a family of random variables @,...,&,, suppose the events 
{x < &, < x,}, k=1,...,m are independent for arbiteary 2 eae 
(x;, < x;,). Then the random variables &,,..., &, are said to be (statistically) 
independent. Given an infinite sequence of random variables &,, &,..., 
suppose the random variables €,,..., &, are independent for every n, or 


equivalently that the events {x, < & < x,}, k =1,2,... are independent 
for arbitrary x;,and x,. Then &, &, ... is said to be a sequence of independent 


random variables. 
Suppose two random variables €, and €, are independent. Then clearly 
their joint probability distribution (4.5) is such that 


Ps (X1, X_) = P, (x1) P;,(%2) (4.7) 
if &, and &, are discrete, and 
Pe, ,é.(%15 X2) = De, (%1)Pe,(Xe) (4.7’) 


if €, and &, are continuous. In (4.7’), Ps, (%41) is the probability density of , 
and p;,(x2) that of €2, while p;, ¢,(%1, X2) is the joint probability density of &, 
and a figuring in (4.6). 


Example 1 (The uniform distribution). Suppose a point & is “‘tossed at 
random” onto the interval [a, b]. This means that the probability of € falling 
in a subinterval [x’, x”] < [a, b] does not depend on the location of [x’, x”]. 
Hence the probability of € falling in [x’, x”] is proportional to the length 
x” — x'.4 More exactly, we have 


x = & 
P x’ = “< x” ae = 
: b—a a eae ae = 
since then the probability of ¢ falling in [a, 5] itself is 


Pig c teem =['< 


as it must be. Clearly, & is a continuous random variable, with probability 


*Let f(s) be the probability of € falling in a subinterval of length s. Then clearly 
f(s + t) = f(s) + f(@). But it can be shown that any function f(s) satisfying this equation 
is either of the form ks (& a constant) or else unbounded in every interval (see W. Feller, 
op. cit., p. 459). 
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density 
1 
p(x) — ({b—a 
0 P= <a Gr x > 0. 
Such a random variable is said to have a uniform distribution. 


Yenc x hb. 


Example 2. Suppose two points € and & are tossed at random and 
independently onto a line segment of length L. What is the probability that 
the distance between the two points does not exceed /? 


Solution. Imagine that €, falls in an interval [0, L] of the x-axis, while 
€. falls in an interval [0, L] of the x-axis, perpendicular to the x,-axis as in 
Figure 5. Then the desired probability is just the probability that a point 
& = (&,, &) tossed at random onto the square 0 < x1, x, < L will fall in the 


FIGURE 5 


region B bounded by the lines x, =/+ x, and x, = —]+ x, (B is the 
unshaded region in Figure 5).> By hypothesis, the random variables €, and &, 
are independent and are both uniformly distributed in [0, L], i.e., both have 
probability density 


1 
<<) =, a ea. 7. 
oo ee 
Hence, by (4.6), the joint probability density of the independent random 


variables &, and &, is just 


1 
Pé,,£(X1 X2) = 12? Oxy Me f. 


> Note that |, — &| is the horizontal distance between the point (€, 2) and the line 
X_ = X,. This is the distance p shown in Figure 5, from which it is apparent that e </ if 
and only if (€1, &) lies in B. 
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Therefore the probability of the random point § = (&, &) falling in the 
region B is given by 
d : 2LI — I’ 
Pié.%)<8) || 
B L 
since L? — 2-3(L — /)*? = 2L/ — /* is the area of B (the square minus the 
two shaded triangles). 


Example 3 (Buffon’s needle problem). Suppose a needle is tossed at random 
onto a plane ruled with parallel lines a distance L apart, where by a “‘needle”’ 
we mean a line segment of length / < L. What is the probability of the needle 
intersecting one of the parallel lines? : 


Solution. Let @, be the angle between the needle and the direction of the 
rulings, and let & be the distance between the bottom point of the needle 
and the nearest line above this point [see Figure 6(a)]. Then the conditions 


FIGURE 6 


of the “‘needle tossing experiment’’ are such that the random variable &, is 
uniformly distributed in the interval [0, 7], while the random variable ¢, is 
uniformly distributed in the interval [0, L]. Hence, assuming that the random 
variables &, and & are independent, we find that their joint probability density 
is 
1 
Dé, ,64(%1 X2) = — , 0< x1< 7, O < Xp <2 
TL 

The event consisting of the needle intersecting one of the rulings occurs if 
and only if 


és < l sin Ea 


i.e., if and only if the corresponding point & = (&,, &,) falls in the region B, 
where B is the part of the rectangle 0 < x, < 7,0 < x, < L lying between 
the x,-axis and the curve x, = sin x, [B is the unshaded region in Figure 
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6(b)]. Hence, by the general formula (4.6), 


P{(Es, &4) € B} =[f 2s — = (4.8) 


where 


1| "sin 0%, = 21 


is the area of B. 

In deducing (4.8), we have assumed that @, and &» are independent 
random variables. This assumption can be tested experimentally. In fact, 
according to (4.8), if the needle is repeatedly tossed onto the ruled plane, 
then the frequency of the event A, consisting of the needle intersecting one of 
the rulings, must be approximately 2//xL. Suppose the needle is tossed n 
times, and let n(A) be the number of times A occurs, so that n(A)/n is the 
relative frequency of the event A. Then 


n(A) as 23 
n TL 
for large n, as discussed on p. 3. Hence 


should be a good approximation to 7 = 3.14... for large n. This actually 
turns out to be the case.® 


Example 4. Given two independent random variables &, and &,, with 
probability densities p, (x) and pz,(x2), find the probability density of the 
random variable 

sca ee Ee. 


Solution. By (4.7’), the joint probability distribution of €, and E, equals 
Pe,(%1)Pz,(%2), and hence, by (4.6), 


Ply <n< vy} = |) paGerdpe,%e) dx aXe 


y’ <aytegcy” 
=[F [2G — py) 4] ay. 


Therefore the probability density of the random variable 7 is given by the 
expression 


Pav) =|" Pay — Ps.) ax, 
called the composition or convolution of the functions p;, and pzg.. 


® See J. V. Uspensky, Introduction to Mathematical Probability, McGraw-Hill Book 
Co., Ine. New York (1937), p. 113. 
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For example, suppose & and &, are both uniformly distributed in the 
interval [0, 1], so that they both have the probability density 
‘te ieas ee pF 
P(x) = 
0 ix =o GAS ee 
Then 


['ax = y qt 0S ye 
0 
pO=\{' dx=2—-y if igys2 


0 ay <0 oF ee 
The graph of the density p,(y) is triangular in shape, as shown in Figure 7. 


Ply) 


0 1 £ 


FIGURE 7 


8. Mathematical Expectation 


By the mathematical expectation or mean value of a discrete random 
variable &, denoted by EZ, we mean the quantity 


Eé = xP,(x), (4.9) 


provided that the series converges absolutely.” Here the summation has the 
same meaning as on p. 37, and, as usual, P.(x) = P{€ = x}. Given a 
discrete random variable &, consider the new random variable » = 0(&), 
where ¢(x).is some function of x. Then the mathematical expectation of 7 is 
given in terms of the probability distribution of € by the formula 

00 


En = E9(8) = > 9(x) PAX). (4.10) 


bangers? 


co 
"Le., provided that > |x| P(x) < ©. 
= 60 
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In fact, 7 is a discrete random variable taking only the values y = 9(x), 
where x ranges over all possible values of the random variable . Therefore® 
Piney) = 2, Phx), 

x:e(a2)=y 


where the summation is over all x such that (x) = y, and hence 


En = Pha = =y}= > eS 2 => oP, 


—o x:9(x)=y 
as asserted. 

More generally, let ~(&, &,) be a random variable which is a function of 
two random variables &, and &, with joint probability distribution 
Pe ,5¢,(X1, Xg). Then it is easily verified that 9(&, €.) has the mathematical 
expectation 


ad 


Eo(&1, &) = bY = P(X, Xq)P sz, 2.(X1 X9). (4.11) 
It is clear from (4.9) that 


ay =i 
b) E(c&) = cEé for an arbitrary constant c; 


C) |EE| < E[€. (4.12) 


Moreover, it follows from (4.11) that 
d) E(é, + &) = E&, + E&, for arbitrary random variables €, and & 
with mathematical expectations E&, and Eé,; 
e) § > O implies EZ > 0, and more generally 


E, < & implies E&, < E&,; | (4.13) 
f) If €, and €, are independent random variables, then 
| aR Ee, Eé.. (4.14) 


For example, to prove (4.14), we write 9(&1, &:) = &1&. Then, for inde- 
pendent €, and &», (4.11) implies 


E(1%s) = y Sx X9Ps(%1)Ps,(%) 


se > XP.) > XaPs(%) = EE, Eé,. 


To define the mathematical expectation of a continuous random variable 
&, we first approximate € by a sequence of discrete random variables @,, 


8 The = should be read as “‘such that.”’ 
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n= 1,2,... Let e,, 2 = 1,2, ~.. bea sequence: of positive muimiper eae 
verging to zero. Then for eachn = 1,2,..., let 
ae, X_2\n» X_iins Xo,ns X1 ons Xo ns oi sin (4.15) 


be an infinite set of distinct points such that? 
SUP [Xin — Xx-a,nl = Ens (4.16) 
k 


and let &,, be a discrete random variable such that 
En = Xn if Xp-1.n < i= Xie n° 
It follows that 


[é ae Sa < Ens 
and hence 


len ~ Sel ee ea oe = eee 
as m,n-— > oo. Therefore, by (4.12) and (4.13), 


as m,n —> oo (provided Eé,, exists for all n). But then 
lim Eé,, 


Nn CO 
exists, by the Cauchy convergence criterion. This limit is called the mathe- 
matical expectation or mean value of the continuous random variable é, 
again denoted by Eé. Clearly, 


Ré = lim >» Mul ie ios =. . = oe 


Suppose & has the probability density p,(x). Then, choosing the points 
(4.15) to be continuity points of p(x), we have 


fe.6) = ie.e) te,n 
Tie ei = pa ie PA) dx 


CO 
> XPM OE — 
—0oo 


and hence 
Ez =|" xp) dx (4.17) 


[compare this with (4.9)]. For continuous random variables of the form 


® The symbol sup denotes the supremum or least upper bound. Therefore the left-hand 
side of (4.16) is the least upper bound of all the differences [Mere — Xe-1.0ls ee 
—1,0,1,2,... Thus (4.16) means that no two of the points (4.15) are more than <, 
apart. Note also that any closed interval of length <, contains at least two of the points 
(4.15). 
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1 = 9(&) and y = 9(&1, &), we have 


Eo(2) =| _ o(x)px) dx (4.18) 
and 


BoE, 2a) =] | elas odes gua X2) dy de (4.19) 


by analogy with (4.10) and (4.11), where pz: (x1, X2) is the joint probability 
density of the random variables £, and & . It 1s easily verified that properties 
a)-f) of the mathematical expectation continue to hold for continuous 
random variables (the details are left as an exercise). 


Remark. Other synonyms for the mathematical expectation of a random 
variable &, discrete or continuous, are the expected value of & or the average 
(value) of &. The mathematical expectation and mean value are often simply 
called the “expectation” and ‘“‘mean,”’ respectively. 


Example Let & be a random variable uniformly distributed in the interval 
[a, b], i.e., let & have the probability density 


if a<x<b, 
p(x) = {b—a 


0 eo yee of x Sb: 


Then the mathematical expectation of & is 


bx dx 8202 


a 
4b —aQ - 


A random variable of the form y = &, + i& involving two real random 
variables €, and &, (the real and imaginary parts of y) is called a complex 
random variable. The mathematical expectation of 7 = €, + i&, is defined as 


Ey = E&, + iE&. 


It is easy to see that formulas (4.10) and (4.18) remain valid for the case 
where 9(&) is a complex-valued function of a real random variable &, and 
that (4.11) and (4.19) remain valid for the case where 9(&, &) is a complex- 
valued function of two real random variables & and & . In particular, let 
0,(&) and .(&) be complex-valued functions of two independent real 
random variables &, and &. Then, choosing 9(&, &2) = 1(&1)@2(&) in 
(4.11) or (4.19), we deduce the formula 


E[91(&1) ¢2(2)] = Eo, (S)E92(&.), (4.20) 


which generalizes (4.14). 
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9. Chebyshev’s Inequality. The Variance and 
Correlation Coefficient 


By the mean square value of a (real) random variable € is meant the 
quantity E@?, equal to 


= 3 x*P (x) 
if & is discrete, or s% 
Ee =|" x p:(x) dx 
if & is continuous.!° Given any random variable and any number < > 0, 


let 
0 if. |é|<e 


ie ee Se 
Then obviously &, < &, and hence, by (4.13), 
Ee, < Ee, 


e*P{|5| > &} < Ee, 


or equivalently 


since clearly 
Ee, = ©*P{|6| > ¢}. 
It follows that 


P{lE|>e} < ug (4.21) 


a result known as Chebyshev’s inequality. According to (4.21), if EZ?/<? < 38, 
then P {|&| > <} < 6, and hence P {|&| < ¢} > 1 — 6. Therefore, if 8 is 
small, it is highly probable that |€| < e«. In particular, if E&? = 0, then 
P {|¢| > ¢} = 0 for every <« > 0, and hence & = 0 with probability 1. 

By the variance or dispersion of a random variable €, denoted by Dé, 
we mean the mean square value E(& — a)? of the difference € — a, where 
a = Eé is the mean value of &. It follows from 


E(& — a)? = E&? — 2aEé + a = EE? — 202 4+ 


that 
DE = EE? — a’. 
Obviously 
Di-= 9%, 
and 
D(c&) = c?DE 


for an arbitrary constant c. 


*° It is assumed, of course, that E@? exists. This is not always the case (see e.g., Problem 
24, p. 33). 
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If , and € are independent random variables, then 
D(é1 ar a) = Dei — Dé». 
In fact, if a, = E&, and a, = E&,, then, by (4.14), 
E(é, — @)(& — a2) = E(&, — a,)E(&, — a.) = 0, (4.22) 
and hence 
D(Ex + 2) = eG + Ee = 1 =)" 
= E(&; — a)? + 2E(é — a1)(S2 — ae) + E(&: — ay)? 
= E(E, — a,)? + E(&, — a)? = DE, + D&. 
Given two random variables €, and €, we now consider the problem of 
finding the linear expression of the form é, + é,&s, involving constants ¢, 


and é, such that ¢, + é,%» is the best “mean square approximation” to &, 
in the sense that 


Ete; —4 = Ae) = min E(E, — 4 — Cate)’, (4.23) 


where the minimum is taken with respect to all c, and c,. To solve this 
problem, we let 


a, = EG, o; = Dé,, 
a, = Eé., on = Des, 
and introduce the quantity 


_ ECGs — abo — 4s) (4.25) 


0199 


(4.24) 


called the correlation coefficient of the random variables €, and @. Going 
over for convenience to the “‘normalized’’ random variables 
ey oa ES 
SSS Ye 


O01 O» 


a ee 
otis ° 


we find that 
min E(E, — cy — ¢,&)" = o; min E(n; — ¢; — Ceo)” (4.26) 


C1,C2 C1,C2Q 

(why ?). Clearly, 

Ej = Eq, = 0, Dn, = En = Dyp = Eqn =1, En =r, 

E(q1 — Tanz = Enrne — rEn2 = 0, 

E(q, — me)” = En — 2rEnine + Ena = 1 — r° 
and hence 

E(q1 — ¢1 — Cae)? = E[(mi — M2) — C1 + (7 — C2) nel? 
= (1 — 7?) + of tr — op" 
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for arbitrary c, and c,. It follows that the minimum of E(y, — ¢, — c27%2)? 
is achieved for c, = 0, c, = r, and is equal to 
min E(y,; a Cy aa stink = 1 oa r”. (4.27) 


1,02 


But 
1 O 
oS sie ~|& aes teas r— (E> = a.) | 
0} Oo 


in terms of the original random variables &, and &. Therefore 


é, + @8 = a, + r 1 (bs — Ae) 
Og 
is the minimizing linear expression figuring in (4.23), where a, a2, oj, o5 and 
r are defined by (4.24) and (4.25). 
If €, and & are independent, then, by (4.22), 


ee E(E1 — 41)(& — 42) = E(E, — aiJE(S: — a) =" 
0199 0199 
It is clear from (4.27) that r lies in the interval —1 < r < 1. Moreover, if 
r = +1, then the random variable &, is simply a linear expression of the 
form 
G1 = C1 + C285. 

In fact, if r= +1, then, by (4.26) and (4.27), the mean square value of 
E, — 6, — 6,8, is just 

E(E, — é, — G&)" = of(1 — r*) = 0, 
and hence &, — ¢, — 6,& = 0 with probability 1 (why ?). 

The above considerations seem to suggest the use of r as a measure of the 
extent to which the random variables &, and & are dependent. However, 
although suitable in some situations (see Problem 15, p. 67), this use of r 
is not justified in general (see Problem 19, p. 53).¥ 


PROBLEMS 


1. A motorist encounters four consecutive traffic lights, each equally likely to 
be red or green. Let € be the number of green lights passed by the motorist 
before being stopped by a red light. What is the probability distribution of ¢? 


2. Give an example of two distinct random variables with the same distribution 
function. 


3. Find the distribution function of the uniformly distributed random variable 
€ considered in Example 1, p. 40. 


11 See also W. Feller, op. cit., p. 236. 
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4. A random variable & has probability density 


a 
PLO =a (— = x ok 2). 
Find 


a) The constant a; b) The distribution function of @; 
c) The probability P{—1 < & < 1}. 


P44 1 
Ans. a)—-; b)=+-arctanx; c) =. 
°19 2 ™ 2 


5. A random variable & has probability density 


ae iP Ad a e500, 
cae 0 otherwise 
(k > 0). Find 
a) The constant a; b) The distribution function of @; 
c) The probability P{O < & < 1/k}. 


6. A random variable & has distribution function 


x 
®.(x) =a + barc tan; (—o <x < o), 
Find 
a) The constants a and 6; b) The probability density of &. 


7. Two nonoverlapping circular disks of radius r are painted on a circular table 
of radius R. A point is then “tossed at random”’ onto the table. What is the 
probability of the point falling in one of the disks? 


Ans 20, RY. 


8. What is the probability that two randomly chosen numbers between 0 and 1 
will have a sum no greater than 1 and a product no greater than §? 


1 2 2/3 dx 
Ans. ~+- — 
3 ie: x 


2 in? 0.49 
ae XY UG. 


9. Given two independent random variables €, and 9, with probability densities 
se 8. tox > @, oe a 2 Se 0) 


i) = K) = 
ee) ; if By Pex\ 0 if x <0, 


find the probability density of the random variable y = & + &. 


Ans. p,(x) = F 


e@/3(1 —e-#/6) if x >0, 
tg = ®t. 
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10. Given three independent random variables €,, & and &3, each uniformly 
distributed in the interval [0, 1], find the probability density of the random 
variable &, + & + &g. 

Hint. The probability density of & + & (say) was found in Example 4, 
p. 43. 


11. A random variable & takes the values 1, 2,...,,... with probabilities 


a | 
re ee 
Find Eé. 


12. Balls are drawn from an urn containing w white balls and 6 black balls until 
a white ball appears. Find the mean value m and variance o? of the number of 
black balls drawn, assuming that each ball is replaced after being drawn. 

b = bw + b) 


= 2 ae 
Ans. m=-—,6 : 


Ww Ww 


13. Find the mean and variance of the random variable & with probability 
density 
pz (x) = dell (—«2o <x < o), 


Ans. EC =O De = 2. 
14. Find the mean and variance of the random variable € with probability 
density 
] 
se Ae a Oy 
p(x) = 47 

0 otherwise. 

Ans, ERE-= a4, BE =073. 


15. The distribution function of a random variable & is 


0 a es Ss 
®,(x) = (a + barc sin x i =) eee 
1 ae ee 


Find Eé and Dé. 
Hint. First determine a and b. 
Ans. E® =.0, Dé. = 4: 


16. Let ¢ be the number of spots obtained in throwing an unbiased die. Find 
the mean and variance of 6. 
Ans. E& =$, DE = 33. 


17. In the preceding problem, what is the probability P of € deviating from 
Eé by more than $? Show that Chebyshev’s inequality gives only a very crude 
estimate of P. 
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18. Prove that if — is a random variable such that Ee®> exists, where a is a 
positive constant, then 


Ee” 
PUES 3 <r 


Hint. Apply Chebyshev’s inequality to the random variable y = e%/?. 


19. Let € be a random variable taking each of the values —2, —1, 1 and 2 with 
probability 4, andlet 7 = &?. Provethat 6 and 7 (although obviously dependent) 
have correlation coefficient 0. 


20. Find the means and variances of two random variables €, and &» with joint 
probability density 


Tv 
sin xX, SIN X»g ft fee, as 
Pe,,€,(%1) Xo) = 
0 otherwise. 


What is the correlation coefficient of &, and &? 


21. Find the correlation coefficient r of two random variables €, and & with 
joint probability density 


1 4 7 
=~ sin (x, + x a. 022 0 Sy =; 
iG, x) - 5 (x, 2) 1<5 255 
0 otherwise. 
7 rita 
1 name 
7 2 16 1 
ns. r= == yY— A z 
ois 


22. Given a random variable &, let p(t) be a nondecreasing positive function 
such that E¢(¢) exists. Prove that 


m 
P{E>f} <TH: (4.28) 


23. Deduce Chebyshev’s inequality as a special case of (4.28). 
24. Let € be a random variable with probability density 


Pe(x) = xl ea (— © 5 ae 00), 


Show that E& and Dé fail to exist. 


> 


THREE IMPORTANT 
PROBABILITY DISTRIBUTIONS 


10. Bernoulli Trials. The Binomial and Poisson Distributions 


By Bernoulli trials we mean identical independent experiments in each of 
which an event A, say, may occur with probability 


P=F 
(p # 0) or fail to occur with probability 
q=1—p. 


“4 


Occurrence of the event A is called a “‘success,’’ and nonoccurrence of A 
(i.e., occurrence of the complementary event A) is called a “‘failure.”’ 

In the case of n consecutive Bernoulli trials, each elementary event 
can be described by a sequence like 


1011... 0001 
n times 
consisting of n digits, each a 0 or a 1, where success at the ith trial is denoted 
by a 1 in the ith place and failure at the ith trial by a 0 in the ith place. Be- 


cause of the independence of the trials, the probability of an elementary 
event w in which there are precisely k successes and n — k failures is just 


P(w) = p*qr™. 
Clearly, the various elementary events are equiprobable only if p = q. 


Now consider the random variable € equal to the total number of suc- 
cesses in n Bernoulli trials, i.e., &(@) = & if precisely k successes occur in the 


54 
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elementary event w. The number of distinct elementary events with the same 
total number of successes k is just the number of distinct sequences consisting 
of k ones and n — k zeros. But the number of such sequences is just the 


binomial coefficient 
n n! 
Cr = — —————__ 54 
: (1) k!(n — k)! ) 


equal to the number of combinations of n things taken k at a time (recall 
Theorem 1.3, p. 7). These C/’ elementary events all have the same probability 


P(w) ag aaa 
and hence the event {& = k} has probability 
Thus the probability distribution of the random variable & is given by 


Pa tp a be 0, 1,5... 8, (5.2) 


and is known as the binomial distribution. The binomial distribution is 
specified by two parameters, the probability p of a single success and the 
number of trials n. 

It should be noted that the random variable & is the sum 


CS esr (5.3) 


of n independent random variables @,,...,&,, where ¢, = 1 if “‘success” 
occurs at the kth trial and €, = 0 if “failure” occurs at the kth trial. We have 


| Dé,, = Eé* es Gey aa ae p” = p(l — p) = pq. 
Therefore 
EE=np, D&E =npga. (5.4) 


Suppose the number of trials is large while the probability of success p 
is relatively small, so that each success is a rather rare event while the average 
number of successes np is appreciable. Then it is a good approximation to 
write 


k 
P(k) ~ - pe 0 1598. 24 (5.5) 


where a = np is the average number of successes and e = 2.718... is the 
base of the natural logarithms. In fact, we know from calculus that 


fs (1 S *) — ot. 


1 mae®) 
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But p = a/n, and hence (5.2) gives 
nN a —O 
PO) = q” = (1°) ~e", 
n 
Moreover, it is easily found from (5.1) and (5.2) that 
PARK) — omp (Rk —1p ee 


Pe Say kq k 
as n— oo. Therefore 


a a 
Pe ; 


2 
a a cate 
fie ee 5 


k 
a =p 
Pek) —~T Pek — . 4, 


which proves the approximate formula (5.5). 
A random variable & taking only the integral values 0,1,2,... is said 
to have a Poisson distribution if 


k 
P,(k) = a ee (5.6) 


The distribution (5.6) is specified by a single positive parameter a, equal to the 
mean value of &: 


a = Ee = Fran 
=0 
In fact, it follows from the expansion ; 


valid for all x, that 


00 00 a® 00 a 
Be == hP (ky => ke = 6 eee ae 
. 2 ) 2 k! xk D! 

Remark. Thus the approximate formula (5.5) shows that the total 
number of successes in ” Bernoulli trials has an approximately Poisson 
distribution with parameter a = np, if n is large and the probability of 
success p is small. 


Example 1 (The lottery ticket problem). How many lottery tickets must 
be bought to make the probability of winning at least P? 
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Solution. Let N be the total number of lottery tickets and M the total 
number of winning tickets. Then M/N is the probability that a bought ticket 
is one of the winning tickets. The purchase of each ticket can be regarded as 
a separate trial with probability of “success’’ p = M/N in a series of n inde- 
pendent trials, where v is the number of tickets bought. If p is relatively small, 
as is usually the case, and the given probability P is relatively large, then it is 
clear that a rather large number of tickets must be bought to make the 
probability of buying at least one winning ticket no smaller than P. Hence 
the number of winning tickets among those purchased is a random variable 
with an approximately Poisson distribution, i.e., the probability that there 
are precisely k winning tickets among the n purchased tickets is 


gS 
Bees 3 


where 


The probability that at least one of the tickets is a winning ticket is just 
L--7fQ).—1—2-. 


Hence n must be at least as large as the smallest positive integer satisfying 
the inequality 
gra grein = 1 P. 


Example 2 (The raisin bun problem). Suppose N raisin buns of equal size 
are baked from a batch of dough into which n raisins have been carefully 
mixed. Then clearly the number of raisins will vary from bun to bun, 
although the average number of raisins per bun is just a = n/N. What is the 
probability that any given bun will contain at least one raisin? 


Solution. It is natural to assume that the volume of the raisins is much 
less than that of the dough, so that the raisins move around freely and 
virtually independently during the mixing, and hence whether or not a 
given raisin ends up in a given bun does not depend on what happens to the 
other raisins. Clearly, the raisins will be approximately uniformly distri- 
buted throughout the dough after careful mixing, i.e., every raisin has the 
same probability 


p= 


Z| 


of ending up in a given bun.! Imagine the raisins numbered from 1 to n, 


1 If v is the volume of the raisins and V that of the dough, then p = v/V. 
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and select a bun at random. Then we can interpret the problem in term as of 
series of n Bernoulli trials, where “‘success”’ at the kth trial means that the 
Ath raisin ends up in the given bun. Suppose both the number of rolls N and 
the number of raisins 7 are large, so that in particular p = 1/Nis small. Then 
the number of successes in the v trials, equal to the number of raisins in the 
given bun, has an approximately Poisson distribution, i.e., the probability 
P(k) of exactly k raisins appearing in the bun is given by 


k 
a —a 
P(k) aww) k! e 3 
where 


a=np=~— 
7 
Hence the probability P of at least one raisin appearing in the bun is 
P=1— 20 =1—<¢™. 


Example 3 (Radioactive decay). It is observed experimentally that 
radium gradually decays into radon by emitting alpha particles (helium 
nuclei). The interatomic distances are large enough to justify the assumption 
that (the nucleus of) each radium atom disintegrates independently of all the 
others. Moreover, each of the nm) radium atoms initially present clearly has 
the same small probability p(t) of disintegrating during an interval of ¢ 
seconds.* Suppose the disintegration of each radium atom is interpreted as 
a “success.” Then the random variable &(t), equal to the number of alpha 
particles emitted in ¢ seconds, equals the number of successes in a series of 
ny Bernoulli trials with probability of success p(t). The values of m, and p(t) 
are such that the distribution of &(¢) is very accurately a Poisson distribution, 
l.e., the probability of exactly k alpha particles being emitted is given by 


k 
PLE) = Kh =e, Pid, (5.7) 


where 
a = EE(t) = mpl) 


is the average number of alpha particles emitted in t seconds. 

Here we have used a model involving Bernoulli trials as a tool for showing 
that the random variable &(t) has a Poisson distribution. Another physical 
situation leading to a Poisson distribution is considered in Example 4, p. 73. 


* A gram of radium (m) ~ 10”) emits about 10!° alpha particles per second. Hence 
p(l) = 1019/10? = 10-12, 
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Il. The De Moivre-Laplace Theorem. The Normal Distribution 


Next we prove the following basic “‘limit theorem’: 


THEOREM 5.1 (De Moivre-Laplace theorem). Given n independent 
identically distributed random variables &1,..., &,, each taking the value 
1 with probability p and the value 0 with probability q = 1 — p, let 


se S, — ES 
a ee Ss. — — iG 
2 5 ./ DS, 
Then 
he xe SF th al er? dy. (5.8) 
n> co ae 


Proof. S,, is the random variable denoted by € in (5.3) and (5.4), 
i.e., S, is the number of successes in » Bernoulli trials, with mean and 
variance 


ES, = np, DS, = npq. 


Hence the ‘‘normalized sum” S* is a random variable taking the values 


ee SOT, An 
| - VnPg 
with probabilities 
ee 5 pg = py be Gs 4... tt. 
r . ki(n — k)! a 


These values divide the interval 


n n 
into n equal subintervals of length 


1 


Ax = g 
Jnpq 


Clearly, as n — ©0, 


k = np + V\npqx— o, n—k=ng—~npqx— o, 
where the convergence is uniform in every finite interval x’ < x < x’. 
Using Stirling’s formula (see p. 10), we find that 


j2un ne” . 


Sink ke" JIn(n — B) (n — byte PS 


oe orice 


P,(k) ~ 
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Moreover, 
Batt [4s nF [2s 
np np nq nq 


Therefore, using the expansion 


x 
In (1 + thy) ~ Oy = > 


(as «, 0), we have 


ae ee 
In (-) = —kln ( -+- 4. x) 
np np 


aes q fo ee 
~~ — npqx seers ao X |, 
(np + /mpax) ( = ae 


ie = = i ( = ‘Pp *) 


— 1 
~ —(ng — Jnpq x) - [E+ —— ~*) 
nq 2 nq 
Adding these expressions, we find that 
k n—k 2 
lim 1a 9, rag? ee 
N—» k n—k 2 
and hence 
.- Sapieeng ae" 
cer tae 
n—-> k i k 
uniformly in every finite interval x’ < x < x”. Since 
wee 
n(n — k) np-‘nq  ./npq : 
it follows that 
lim P {35 =x} = == el Ax, Ax = = 
n> 00 J2n npq 


Therefore 


n 
n> co n?>o << 2" 


iim Bix < So) Sli eX} 


ei a 1 ett Ay, (5.9) 


no v<xr<a” «/ 2% 
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where the sum is over all values of x in the interval x’ < x < x”. But 


clearly 
a Pe" ax (5.10) 


n~o 2 <ck2" 27 xf 7 


(why ?). Comparing (5.9) and (5.10), we finally get the desired limiting 
formula (5.8). Jj 


According to Theorem 5.1, the limiting distribution of the random 
variable S* is the distribution with probability density 


p(x) = 1 ot. (5.11) 
28. 


Such a distribution is called a normal (or Gaussian) distribution. The density 
p(x) is the “bell-shaped” curve shown in Figure 8(a). The corresponding 
distribution function is 


gph ("ule 
=e Je du, (5.12) 


p\x)= er l2 _w<x<@ 


7 


Xx 
aut 
/2 du,-O<x<@O 


FIGURE 8 
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Table 2. Values of the normal 
distribution function ®(x) given 
by formula (5.12). 
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and is shown in Figure 8(b). Since p(x) is even, it is clear that 
O(—x) = 1 — O(x). 
Representative values of (x) are given in Table 2. 
Let & be anormal (or Gaussian) random variable, i.e., a random variable 
with probability density gee Then 


EE = xe? dx = 0, 


ral 
since the integrand is odd. Moreover, 


ie.@) 


ar li 


Lae 
xe * dx 


DE = Eé? — (Eé)’ = ES’ = 


Stn is xe? dy, 
gl 27] N70 ¥—N 
Integrating by parts, we get 


T No 


—— | paeee ee aes a ax| 
wen N-> oo a=—N —N 
se 
_ ee 
Hence & has variance 1. More generally, the random variable with probability 
density 


Dé = = lim [= i ae xe") 
ie el dx = 1, 


ee ee (5.13) 
TO 


is also called a normal random variable, and has mean a and variance o? 
(show this). 


Example (Brownian motion). Suppose a tiny particle is suspended in a 
homogeneous liquid. Then the particle undergoes random collisions with the 
molecules of the liquid, and, as a result, moves about continually in a 
chaotic fashion. This is the phenomenon of Brownian motion. As a model of 
Brownian motion, we make the following simplifying assumptions, charac- 
terizing a “discrete random walk” in one dimension: 


1) The particle moves only along the x-axis. 

2) The particle moves only at the times ¢ = nAt,n =0,1,2,... 

3) Suppose the particle is at position x at time ¢. Then, regardless of its 
previous behavior, the particle moves to either of the two neigh- 
boring positions x + Ax and x — Ax (Ax > 0) with probability 3 
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In other words, at each step the particle undergoes a shift of amount 
Ax either to the right or to the left, with equal probability.® 


Now let &(t) denote the position of our “Brownian particle” at time f, 
and suppose the particle is at the point x = 0 at time t = 0, so that &(0) = 0. 
Then after t = nAt seconds, the particle undergoes n displacements of amount 
Ax, of which S,,, say, are to the right (the positive direction) and n — S,, to 
the left (the negative direction). As a result, the position of the particle at 
time ¢ = nAt is just 


E(t) = [S, Ax — (n — S,) Ax] = (2S,, — n) Ax. (5.14) 
Moreover, since &(0) = 0, we have 


&() = [&(s) — §(0)] + [E® — &G)I 


for any sin the interval 0 < s < ¢ (for the time being, sis an integral multiple 
of Ax). With our assumptions, it is clear that the increments &(s) — &(0) 
and &(t) — &(s) are independent random variables, and that the probability 
distribution of &(t) — &(s) is the same as that of &(t — s) — &(0). Therefore 
the variance o?(t) = Dé&(ft) satisfies the relation 


o*(t) = o%(s) + o%(t — 5), O-<— 2 = f. 
It follows that o?(¢) is proportional to 1, i.e.,* 
De(t) = of, (5.15) 


where o? is a constant called the diffusion coefficient. On the other hand, it is 
easy to see that after a time ¢, 1.e., after n = t/At steps, the variance of the 
displacement must be 


t 
: DEQ) === (Ay. 5.16 
a) = a) (5.16) 
Comparing (5.15) and (5.16), we obtain 


3 =n 


5 (5.17) 


The displacements of the particle are independent of one another and can 
be regarded as Bernoulli trials with probability of “success” p = 4, ‘“‘success”’ 
being interpreted as a displacement in the positive direction. In this sense, 
the number of displacements S,, in the positive direction is just the number of 


* We will eventually pass to the limit Ar + 0, Ax — 0, thereby getting the “continuous 
random walk” characteristic of the actual physical process of Brownian motion. 
* See footnote 4, p. 40. 
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successes in ” Bernoulli trials. Moreover, the relation between the particle’s 
position at time ¢ and the normalized random variable 


ae eS ae 


23 = ya 
is given by 


E(t) = St/nAx = S*/1 =e = stov/i, 


J 
because of (5.14) and (5.17). Applying Theorem 5.1, in particular formula 
(5.8), and passing to the limit At— 0 while holding o constant (so that 
Ax — 0), we find that the random variable &(t) describing the one-dimen- 
sional Brownian motion satisfies the formula 


1 
P{’ < aD "| =n Pi < Ss; < ye | 
o,/t At+0 Jon 


Therefore &(t) is a normal random variable with probability distribution 


Pie = oe = | Mee PB dy, 


a’ 


1 
om 


PROBLEMS 


1. Consider the game of “heads or tails,’ as in Example 3, p. 29. Show that 
the probability of correctly calling the side of the coin landing upward is always 
3 regardless of the call, provided the coin is unbiased. However, show that if 
the coin is biased, then “‘heads’’ should be called all the time if heads are more 
likely, while “tails” should be called all the time if tails are more likely. 


2. There are 10 children in a given family. Assuming that a boy is as likely to 
be born as a girl, what is the probability of the family having 
a) 5 boys and 5 girls; b) From 3 to 7 boys? 


3. Suppose the probability of hitting a target with a single shot is 0.001. What 
is the probability P of hitting the target 2 or more times in 5000 shots? 
Ans. P ~ 1 — 6e~° = 0.96. 


4. The page proof of a 500-page book contains 500 misprints. What is the 
probability P of 2 or more misprints appearing on the same page? 


5 
Ans. Px 1—— + 0.08. 
2e 


5. Let p be the probability of success in a series of Bernoulli trials. What is the 
probability P,, of an even number of successes in n trials? 


Ans. P, =%$[1 + (1 — 2p)"). 
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6. What is the probability of the pattern SFS appearing infinitely often in an 
infinite series of Bernoulli trials, if S denotes ‘‘success’’ and F ‘‘failure’’ ? 


Hint. Apply the second Borel-Cantelli lemma (Theorem 3.1, p. 33).° 
Ans. ~1. 


7. An electronic computer contains 1000 transistors. Suppose each transistor 
has probability 0.001 of failing in the course of a year of operation. What is the 
probability of at least 3 transistors failing in a year? 


8. Aschool has 730 students. What is the probability that exactly 4 students were 
born on January 1? 


Hint. Neglect leap years. 
9. Let € be a random variable with the Poisson distribution (5.6). Find 


— 

a) ot Dies 
1 

Ans. a) @; se ae 


10. Where is the uniform convergence used in the proof of Theorem 5.1? 


11. The probability of occurrence of an event A in one trial is 0.3. What is the 
probability P that the relative frequency of A in 100 independent trials will lie 
between 0.2 and 0.4? 


Hint. Use Theorem 5.1 and Table 2. 
Ans. P = 0.97. 


12. Suppose an event A has probability 0.4. How many trials must be performed 
to assert with probability 0.9 that the relative frequency of A differs from 0.4 by 
no more than 0.1? 


Ans. About 65. 


13. The probability of occurrence of an event A in one trial is 0.6. What is the 
probability P that A occurs in the majority of 60 trials? 


Ans. P » 0.94, 


14. Two continuous random variables &, and & are said to have a bivariate 
normal distribution if their joint probability density is 


1 
2n0,0,V 1 — r 


1 E a _ Gy. 2, ~ 8) a Ue. = allt 


Xp Ss 
P| 2(1 — r?) : om ed) 03 


Pe, E,(%15 %2) = 


fo) 


(5.18) 


* For further details, see W. Feller, op. cit., p. 202. 
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where o, > 0, o, > 0, —1 <r <1. Prove that each of the random variables 
é, and & has a univariate (i.e., one-dimensional) normal distribution of the 
form (5.13), where EE, = a, Dé, = o?, EE, = b, DE, = 95. 


Hint. Clearly, 


Pe, 1) =({". Pe, ,2,(%1> X_) dX, Pe, (2) =(". Pé1,2.(%1> X_) dx, 
(why ?). 


15. Prove that the number r in (5.18) is the correlation coefficient of the random 
variables &, and &. Prove that &, and & are independent if and only ifr = 0. 


Comment. This is a situation in which r is a satisfactory measure of the 
extent to which the random variables &, and & are dependent (the larger |r|, 
the “‘more dependent”’ €, and &»). 


16. Let &, and & be the same as in Problem 14. Find the probability distribution 
of = et + a 


Ans. The random variable 7 is normal, with probability density 


(x) 1 (x — a — b)? 
XS EEE OA SS ee ee Ls 
a V2n(o2 + 2ro,5, + 62) PL 2(63 + 2royoq + 03) 


6 


SOME LIMIT THEOREMS 


12. The Law of Large Numbers 


Consider n independent identically distributed random variables &,..., 
En Mm particular, &,..., &, have the same mean a = E&, and variance 
o” = Dé... If 


n= Ato +) 


is the arithmetic mean of the variables €,,..., &,, then 


Applying Chebyshev’s inequality (see Sec. 9) to the random variable y — a, 
we get the inequality 


1 2 
P {ly — al > ¢} < GEM — a) (6.1) 


for arbitrary « > 0. 


THEOREM 6.1 (Weak law of large numbers). Let &,,...,&, ben inde- 
pendent identically distributed random variables with mean a and variance 
o*. Then, given any 8 > 0 and < > 0, however small, there is an integer 


68 
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n such that 


a—e< “(+ cae e eae 


with probability greater than | — 6. 
Proof. The theorem is an immediate consequence of (6.1) if we choose 
n> o*/de7, fj 


Remark. Suppose 5 and «< are so small that we can practically neglect 
both the occurrence of events of probability 3 and differences between quanti- 
ties differing by no more than c. Then Theorem 6.1 asserts that for sufficiently 
large n, the arithmetic mean 


n= Gite +h) 


is an excellent approximation to the mean value a = E&, 


Now consider n consecutive Bernoulli trials, in each of which an event 
A can occur with probability p = P(A) or fail to occur with probability 
q=1-—p. Let &, be a random variable equal to | if A occurs at the Ath 
trial and 0 if A fails to occur at the kth trial. Then the random variables 
£,,..., &, are independent and identically distributed (by the very meaning 
of Bernoulli trials). Obviously 


P{&=l =p, P {&, = 0} = q. 
Moreover, each random variable &, has mean 

a= EE, = p-1+q-0=p= P(A). 
Let n(A) be the number of trials in which A occurs, so that 

n(A) 
n 
is the relative frequency of the event A. Then clearly 
Mays &, 4 +2 eas 


and hence 


oe 


It follows from Theorem 6.1 that n(A)/n virtually coincides with P(A) for 
sufficiently large n, more exactly, that given any 6 > 0 and « > 0, however 


small, there is an integer m such that 


Pi Abt << -o 
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with probability greater than 1 — 5. The justification for formula (1.2), 
p. 3 is now apparent. 


Remark. It can be shown! that with probability 1 the limit 


hit —_ 


exists and equals P(A). This result is known as the strong law of large 
numbers. 


[3. Generating Functions. Weak Convergence of Probability 
Distributions 


Let € be a discrete random variable taking the values 0,1,2,... with 
probabilities 
P(k) = P {§ = k}, ae, oe et See (6.2) 
Then the function 


Fiz) = ¥ P.(k)z*, lz eet (6.3) 


is called the generating function of the random variable € or of the corre- 
sponding probability distributions (6.2). It follows from the convergence of 
the series (6.3) for |z| = 1 and from Weierstrass’s theorem on uniformly 
convergent series of analytic functions? that F,(z) is an analytic function of z 
in |z| < 1, with (6.3) as its power series expansion. Moreover, the probability 
distribution of the random variable @ is uniquely determined by its generating 
function F,(z), and in fact 


P{k) = — FPO. bosG 87... 


where F*)(z) is Sie kth derivative of F(z). According to formula (4.10), 
p. 44, for fixed z the function F,(z) is just the mathematical expectation of 
the random variable 9(&) = 2°, ie., 


Fay = Ee, Cos 8 (6.4) 


Example 1 (The Poisson distribution). If the random variable & has a 
Poisson distribution with parameter a, so that 


qa” 
P(K)= ee, k=0,1,2,..., 


1 See e.g., W. Feller, op. cit., p. 203. 
*See e.g., R. A. Silverman, Introductory Complex Analysis, Prentice-Hall, Inc., 
Englewood Cliffs, N.J. (1967), p. 191. Also use Weierstrass’ M-test (ibid., p. 186). 
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then & has the generating function 
00 k 00 k 
F(z) ies > be etzk — oo 52) = etlz-D) (6.5) 
x=0 k! z=0 k! 


Suppose the random variable & has mean a = Eé and variance o? = DE. 
Then, differentiating (6.4) twice with respect to z behind the expectation sign 
and setting z = 1, we get 


a=EE=F (1), of = EC? — 2)? = FF’) + FQ) — [FOP 
(6.6) 


The same formulas can easily be deduced from the power series (6.3). In 
fact, differentiating (6.3) for |z| < 1, we get 


a rie", 
k=0 
and hence 
a = > kP-(k) = lim F(z) = F,(1), 
k=0 2 
and similarly for the second of the formulas (6.6). 
Next let &,,..., &, bem independent random variables taking the values 


0,1,2,... Then the random variables z,..., z’n, where zis a fixed number, 
are also independent. It follows from formula (4.20), p. 47 that 


Ezit $n) — Bizht-- + 75m — Ez1--- Ez. 


Thus we have the formula 


F(z) = F;,(2)°** F;,(2), (6.7) 
expressing the generating function F(z) = Ez ofthesum& = &,+:°:°:4+ &, 
of the n random variables ,,..., &, in terms of the generating functions 
F,(z) = Ez*,k =1,...,n of the separate summands. 


Example 2 (The binomial distribution). Suppose the random variable & 
has a binomial distribution with parameters p and n, so that 


Pik) = Cpe, gi =p, k=0,1,...,n. 


Then, as already noted in Sec. 10, € can be regarded as the sum em 
&,+.----+ &, of n independent random variables &,...,&,, where 


1 with probability p, 
"(0 with probability gq. 
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The generating function F; (z) of each summand is clearly pz + q, and hence, 
by (6.7), the generating function of & itself is 


F,(z) = (pz + q)”. (6.8) 

Now let &,, 7» =1,2,... be a sequence of discrete random variables 
taking the values 0,1, 2,..., with probability distributions P,(k) = P; (k) 
and generating functions F,(z), n = 1,2,... Then the sequence of distri- 
butions {P,,(k)} is said to converge weakly to the limiting distribution P(k) if 
fim Pak) = Ft) (6.9) 


for ait ke = 0 


Example 3 (Weak convergence of the binomial distribution to the Poisson 
distribution). Let &, &,... be a sequence of random variables such that 
€, has a binomial distribution P,,(k) with parameters p and 2, i.e., 


fin ore ee 
Suppose p depends on n in such a way that the limit 


limnp=a (6.10) 


n— oO 


exists. Then, according to formula (5.5), p. 55, the sequence of distributions 
{P,,(k)} converges weakly to the Poisson distribution 


a” 
P(ky=— e% k= 0, 12,2. 
k! 
with parameter a given by (6.10). 
In Example 3, the sequence of generating functions 


Feit gr Bee, 


of the random variables &,, &,... converges uniformly to the generating 
function F(z) = e*'*-) of the limiting Poisson distribution, i.e., 


lim F,(z) =lim [1 + p(z — 1)] =lim c 4 mp = oy 


n> co n> oO nN 0O n 
a(z — 1)]|” = 
—lim c + i225 : ee aed 
nN n 


(justify the next-to-the-last step). This is no accident, as shown by 


THEOREM 6.2. The sequence of probability distributions P,,(k), n = 1, 
2,... with generating functions F,(z),n = 1, 2,... converges weakly to 


SEC. 13 SOME LIMIT THEOREMS 73 


the limiting distribution P(k) if and only if 
baa fis) = Ftz), (6.11) 


n> 0 


where 
Hie) => Pz” 
k=0 
is the generating function of P(k) and the convergence is uniform in every 
ek iz ae <1. | 
Proof. First suppose (6.9) holds. Clearly, 
K 
|F,(z) — F(2)l < 2 1P.(k) — P(R)| + v3 Iz|" (6.12) 
k=0 


k=K+1 
for any positive integer K. Given any < > 0, we first choose K large 
enough to make 


> iz < 


k=K+1 —r 2 
and then find a positive integer NV ae that 


|PAk) — P(k)| < 


€ 
2(K + 1) 
holds fork =0,...,Kifn-> UN. It then follows from (6.12) that 
|F,(Z) io F(z)| = 


if n > N, which immediately proves (6.11). 

Conversely, suppose (6.11) holds, where the convergence is uniform 
in every disk |z| < r < 1. Then, by Weierstrass’ theorem on uniformly 
convergent sequences of analytic functions,* 


lm Fy = FG), |zf = 1 (6.13) 


nm CO 


for all & = 0.1, 2.4 Bet 


P,(k) = = 2 (0), = P(k) = <F (0), 
and hence (6.13) implies . for-an k == 6, 1, Saw 


The following example is typical of the situations where the Poisson 
distribution is encountered: 


Example 4 (Random flow of events). Suppose that events of a given kind 
occur randomly in the course of time. For example, we can think in terms 


3. R. A. Silverman, op. cit., p. 192. 
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of “service calls” (requests for service) arriving randomly at some “‘server”’ 
(service facility), like inquiries at an information desk, arrival of motorists 
at a gas station, telephone calls at an exchange, etc. Let €(A) be the number of 
events occurring during the time interval A. Then what is the distribution of 
the random variable &(A)? 

To answer this question, we will assume that our “random flow of events” 
has the following three properties: 


a) The events are independent of one another; more exactly, the random 
variables &(A,), &(A,), ... are independent if the intervals A,, A,,... 
are nonoverlapping. 

b) The flow of events is “‘stationary,”’ i.e., the distribution of the random 
variable €(A) depends only on the length of the interval A and not on 
the time of its occurrence (the initial time of A, say). 

c) The probability that at least one event occurs in a small time interval 
At is AAt + o(At), while the probability that more than one event 
occurs in At is o(At). Here o(At)is an infinitesimal of higher order than 
At, 1.e., 


and A is a positive parameter characterizing the “rate of occurrence”’ 
or “density” of the events. 


Now consider the time interval A = {0,1f], and let &(t) be the total 
number of events occurring in [0,¢]. Dividing [0,f] into n equal parts 
Ay... > 4,, we fad that 


E() =D EAs) 


where ¢(4,),... ; GG are independent random variables and &(A,) is the 
number of events occurring in the interval A,. Clearly, the generating 
function of each random variable &(A,) is 


F(Z) = (1 =S\5 ee o(*), 


where o(¢/n) is a term of order higher than Hs Hence, by (6.7), the a 
function of &(f) is 


F(z) = [F,(2)]" = c # “a” E: 0(*) |. 


But F(z) is independent of the subintervals A,,...,A,, and hence we can 
take the limit as n — oo, obtaining 


F(z) = lim E a aca — git(z-1) 


n> oO n 
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Comparing this with (6.5), we find that F(z) is the generating function of a 
Poisson distribution with parameter a = Af, so that 


P{é(1) = k} Oe Pa... 


Since 
M = EZ(1), 


the parameter A is just the average number of events occurring per unit time. 


14. Characteristic Functions. The Central Limit Theorem 


Given a real random variable &, by the characteristic function of & is 
meant the function 


ie to << t < (6.14) 


Clearly, f.(t) coincides for every fixed ¢ with the mathematical expectation of 
the complex random variable y = e***. For a discrete random variable taking 
the values 0,1,2,..., the characteristic function f.(t) coincides with the 
values of the generating function F;(z) on the boundary of the unit circle 
Pt = 1. 16; 


fut) = Fsle) = B Pale 


This formula represents f,(t) as a Fourier series, with the probabilities 
Pie) 2) 6 k;, k = 0,1, 2,... a8 its coefiicients. Thus these proba- 
bilities P.(k) are uniquely determined by the characteristic function f.(t). 

If € is a continuous random variable with probability density p(x), then, 
by formula (4.18), p. 47, the characteristic function is the Fourier transform 
of the density p,(x): 


fAt) = [Ee dx. (6.15) 


Inverting (6.15), we find that 
1 Pins aie 
p(x) = — { e*f.(t) dt, (6.16) 
27 Y—x 


at least at points where p,(x) is suitably well-behaved.* Thus p,(x) is uniquely 
determined by the characteristic function /,(f). 


* If (6.16) fails, another inversion formula can be used, giving the distribution function 
® .(x) = P{& < x} in terms of the characteristic function f z(t) (see e.g., B. V. Gnedenko, 
op. cit., Sec. 36). We can then deduce p g(x) from ® e(X) by differentiation, at least almost 
everywhere (recall footnote 2, p. 38). 
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Example 1. Let & be a normally distributed random variable, with 
probability density 


5) (6.17) 


27 


Then, by (6.15), the characteristic function of & is given by 


f(t) =|" e%**t n(x) dx = Ze se 


ja 1 [? te 2 

——— te eS 

ee eee 
The function 9(z) = e~*"/? is an analytic function of the complex variable z, 
and hence, by Cauchy’s integral theorem,’ the integral of ¢(z) along the 
rectangular contour with vertices (—N,0), (N,0), (N, —it), (—N, —it) 
equals zero. Therefore 


(6.18) 


=—€ 


I _ oe (zit)? /2 dx = lim 1 f —(x—it)?/2 dx 


af re pt No = 
1 N—it 2 2 
= |j eine e* /* dz = ji ef ax 6.19 
a Jor ae Z oe /2n ae ( ) 
--- a= el dy. 


where we use the fact that the integral of ~(z) along the vertical sides of the 
contour vanishes as N — oo (why?). But 


es: 
Jin 
as for any probability density. Hence (6.18) and (6.19) imply 
f(t) = et? (6.20) 


gee "ak =|" pOoy-ds = 1, 


Now suppose the random variable & is such that E |&|° exists. Then the 
characteristic function f;(t) has the expansion 


fO= 1+ HE 2- mE 24 RU), (6.21) 


where the remainder R(t) satisfies the estimate 


IR@| < CE |&|* - |¢/? 


* R. A. Silverman, op. cit., p. 146. 
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(C denotes a constant). In fact, we need only note that 


2 
et — 1 + i&t — x P46, (6.22) 


by Taylor’s formula, where 
JO] < CEP 2°. 


We then get (6.21) by taking the mathematical expectation of both sides of 
(6.22). In particular, it follows from (6.21) that the mean a = Eg and 
variance co? = Dé are given by the formulas 


a= —if(0), 3 = —fx0) + [fOF. (6.23) 


Example 2. According to (6.23), the normally distributed random variable 
€ with probability density (6.17) has mean : 


a = —if’(0)=0 
and variance 


ea =F (0) = 1. 


Formula (6.7) has a natural analogue for characteristic functions. In 
fact, if 2,,..., &,, are independent random variables with sum § = &, + °°: 
+ &, then, by formula (4.20), p. 47, the characteristic function of & is 


f(t) = fz (t) Sa  fz,(t). (6.24) 

Let &,,n = 1,2,... be a sequence of random variables with character- 
istic functions f,,(t), m = 1,2,... Then the sequence of probability distri- 
butions of &,, &,... is said to converge weakly to the distribution with 


density p(x) if 


Hm PS Es x" } —|* p(x) dx 
for all x’ and x” (x’ < x”). This should be compared with the definition of 
weak convergence for discrete random variables taking the values 0, 1, | Bae 
given in Sec. 13. 
Theorem 6.2 has a natural analogue for characteristic functions, whose 
proof will not be given here:® 


THEOREM 6.2.’ The sequence of probability distributions with charac- 


teristic functions f,(t), n = 1,2,... converges weakly to the limiting 
distribution with density p(x) if and only if 
lim f,(t) = f (1), (6.11’) 


6 For the proof, see e.g., B. V. Gnedenko, op. cit., Sec. 38. 
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where 
f() =| e**p(x) ax 


is the characteristic function of the limiting distribution and the con- 
vergence ts uniform in every finite interval t'’ << t < t”. 


We now prove a key proposition of probability theory, called the central 
limit theorem, which has the De Moivre-Laplace theorem (Theorem 5.1, 
p. 59) as a very special case. Roughly speaking, the central limit theorem 
asserts that the distribution of the sum of a large number of independent 
identically distributed random variables is approximately normal. 


DEFINITION. Given a sequence of random variables ©, k = 1,2,... 
with finite means a, = E&, and variances o% = Dé,, consider the 
“normalized sum’ 


e Lg TAS 
n /DS,, b] 
where 
Sy == = Ea. 
k=1 
Then the sequence &, k =1,2,... is said to satisfy the central limit 
theorem if? 
lim # {x Se eS i "et dx. (6.25) 
n> co ab fai 
THEOREM 6.3. Suppose the sequence of independent random variables 
Ex, kK =1,2,... with means a, and variances o? satisfies the Lyapunov 
condition 
lim — > Ele, — al? = 0, (6.26) 
n>o D, k=1 
where 


n 
B= BS. = >. of. 
k=1 
Then the sequence of random variables satisfies the central limit theorem. 


Proof. Equation (6.25) means that the sequence of distributions 
of the normalized sums S*, n =1,2,... converges weakly to the 


normal distribution with probability density (6.16). Hence, according 
to Theorem 6.2’, we need only show that the sequence of characteristic 


’ Cf. formula (5.8), p. 59. Note that the right-hand side of (6.25) equals B(x”) — ®(x’), 
where ®(x) is the distribution function of a normal random variable with mean 0 and vari- 
ance 1. 
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functions f,(t), 2 =1,2,... of the random variables S* converges 
uniformly in every finite interval t’<t<t" to the characteristic 
function f(t) = e~*’/? of this normal distribution (recall Example 1). 
Clearly, 


op 2 
s =. FP = DS,,. 
2 a, 


The random variable &, — a, has zero mean and variance of, and 
hence, by (6.21), has characteristic function 


2 
g() =1— FP + RO, 


where 
IR.(O)| < Cth E ae 4, |" : 


(C is some constant). Therefore the characteristic function of the 
random variable y, = (&, — @;)/B, is 


t o t 
firlt) = (5) =j1— ane a RF): 


t 
ee) 
It follows from (6.24) that the random variable S*¥ = y, +°'* + Na 
has characteristic function 


n 


where 
E ta es af 


acw 
Bi, 


oe TI fen(t. 


Hence 


in fat) =¥ foal) ~¥ | — FS + BF) | 


where, because of the hypothesis (6.26), 


2™(5) 


as n—> oo uniformly in every finite interval t’ < t < t”. Therefore 


= ares > El, — a|?—-0 
n =] 


t? n ; ¢* 
a = a 


or equivalently 


FO ees 
asn—>o. ff 
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Example 3. The Lyapunov condition is always satisfied if the random 
variables €,, 9, ... are identically distributed and if « = E |&, — a,,|°> exists. 
In fact, 


B i Dé,, sae no’, 
k=1 


where o? = Dé,, and hence 


lim — SE |é, — a,° =lim 


n— Cc 


PROBLEMS 


1. Show that the conclusion of Theorem 6.1 can be written in the form 
; oe | 


2. Let &,..., & be m independent identically distributed random variables, 
with common mean a = E€, and variance €* = D&,. Suppose a is known. Can 
the quantity 


i n 
Nk=1 


n—- © 


for arbitrary ¢ > % 


SS 

Be ee 2 

n A {5 ) 
be used to estimate o?? 


3. A random variable € has probability density® 


x™ 
= "igs TSG, 
7m 
Pa = 
5 0 otherwise, 


where m is a positive integer. Prove that 


m 


4. The probability of an event A occurring in one trial is 4. Is it true that the 
probability of A occurring between 400 and 600 times in 1000 independent 
trials exceeds 0.97? 


Ans. Yes. 
5. Let & be the number of spots obtained in throwing an unbiased die. What is 


the generating function of &? 


* It follows by repeated integration by parts that J ; ea = = 
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6. Use (6.6) and the result of the preceding problem to solve Problem 16, p. 52. 


7. Let & be a random variable with the Poisson distribution 


k 


a 
P(k) = ki e, 2g ae Se ere (6.27) 


Use (6.6) to show that EE = Dé =a. 


8. Find the generating function of the random variable & with distribution 
a” 
~ (i +ajth 


Plt =e (a > 0). 


Use (6.6) to find E& and DE. 


9. Let 7 be the sum of two independent random variables & and &, one with 
the Poisson distribution (6.27), the other with the Poisson distribution obtained 
by changing the parameter a to a’ in (6.27). Show that y also has a Poisson 
distribution, with parameter a + a’. 


10. Let S,, be the number of successes in a series of n independent trials, where 
the probability of success at the kth trial is p,. Suppose p;,..., Pn depend on x 
in such a way that 


ae 
while 


max {p1,..-5 Pn} > 0 


as n—> oo. Prove that S, has a Poisson distribution with parameter in the 
limit as n — ©. 


Hint. Use Theorem 6.2.° 


11. Find the characteristic function f; (t) of the random variable with probability 
density 


1 
pee) =ze  (-0@ <x <2), 
Ans. eee ae 


12. Use (6.23) and the result of the preceding problem to solve Problem 13, 
p. 52. 


13. Find the characteristic function of a random variable uniformly distributed 
in the interval [a, b]. 


14. A continuous random variable & has characteristic function 
fit) = ealtl (a > 0). 
Find the probability density of &. 


* For the details, see W. Feller, op. cit., p. 282. 
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Ans. Pel) = ma? + x2) . 


15. The derivatives f;(0) and fz (0) do not exist in the preceding problem. Why 
does this make sense? 


Hint. Cf. Problem 24, p. 53. 


16. Let v be the total number of spots which are obtained in 1000 independent 
throws of an unbiased die. Then Ev = 3500, because of Problem 16, p. 52. 
Estimate the probability that v is a number between 3450 and 3550. 


17. Let S, be the same as in Problem 10, and suppose > pg, = ©. Prove that 
k=1 


n 
Sn a > Pk 
t k=1 1 2 
ee ee = ee oad +— | e #12 dy 
n Vn x 
> Pr 
k=1 


assn— ©. 


Hint. Apply Theorem 6.3. 
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MARKOV CHAINS 


15. Transition Probabilities 


Consider a physical system with the following properties: 
a) The system can occupy any of a finite or countably infinite number of 


States €,, &,... 

b) Starting from some initial state at time t = 0, the system changes its 
state randomly at the times ¢ = 1, 2,... Thus, if the random variable 
E(t) is the state of the system at time ¢,’ the evolution of the system in 
time is described by the consecutive transitions (or “steps’’) 


Oo - o)> =>": 
c) At time t = 0, the system occupies the state ¢; with initial probability 
pF ee =), i= 1,2,... (7.1) 


d) Suppose the system is in the state «; at any time n. Then the proba- 
bility that the system goes into the state <, at the next step is given by 


Py =P LEG + 1 | = 28, Se eee (7.2) 


regardless of its behavior before the time n. The numbers p,,, called 
the transition probabilities, do not depend on the time n. 


1 Jn calling €(¢) a random variable, we are tacitly assuming that the states ¢), e2,... 


are numbers (random variables are numerical functions). This can always be achieved by 
the simple expedient of replacing ¢,, &2,... by the integers 1, 2,... (see W. Feller, op. cit., 
p. 419). 
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A “random process” described by this model is called a Markov chain.” 
Now let 

p(n) = P {&) = «,} (7.3) 
be the probability that the system will be in the state ¢, “after n steps.” 
To find p,(n), we argue as follows: After n — 1 steps, the system must be in 
one of the states ¢;, K = 1, 2,..., 1.¢., the events (6 — 1) = <,}, k= 
1,2,...form a full set of mutually exclusive events in the sense of p. 26. 
Hence, by formula (3.6), 


P {&(n) = ej} = == {E(n) = «,; | E(n — 1) =e} P {&(n — 1) =e}. (7.4) 


Writing (7.4) in terms of the notation (7.1)-(7.3), we get the recursion 
formulas 

p,(0) = Dis 

p(n) = 2, Pin oar 1) Dy 55 ae z; ee Jee. (7.5) 


If the system is in a definite state <, at time t = 0, the initial probability 
distribution reduces to 


Peet do S90 Ki. (7.6) 
The probability p,(n) is then the same as the probability 


Pis(n) = P {6) = ¢; | E(O) == @2 1 fa Dee ee 


that the system will go from state <,; to state <, in n steps. Hence, for the 
initial distribution (7.6), the formulas (7.5) become 


a eo 
P;(0) = ae : 
; ce eee (7.7) 


p(n) = Ze Pix — 1)Dp3; ae oe ee 
k 


The form of the sum in (7.7) suggests introducing the transition probability 
matrix 
Pui Pie 
f= 


2 More exactly, a Markov chain with stationary transition probabilities, where we allude 
to the fact that the numbers (7.2) do not depend on n. For an abstract definition of a 
Markov chain, without reference to an underlying physical system, see W. Feller, op. cit., 
p. 374. 
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and the “‘n-step transition probability matrix”’ 


Pult) Pi2(n) 
P(n) = |lpy)|| = _ a 


Then, because of the rule for matrix multiplication,° (7.7) implies 
ee 2 PZ) = PO)P =F"... 


where J is the unit matrix (with ones along the main diagonal and zeros 
everywhere else). It follows that 


PQ) =P", n= 1,2,".. (7.8) 


Example 1 (The book pile problem). Consider a pile of m books lying 
on a desk. If the books are numbered from 1 to m, the order of the books 
from the top of the pile down is described by some permutation (i,, i2,... , 
i,,) of the integers 1,2,...,m, where i, is the number of the book on top of 
the pile, 7, the number of the next book down, etc., and i,, is the number of 
the book at the bottom of the pile. Suppose each of the books is chosen with 
a definite probability, and then returned to the top of the pile. Let p, be 
the probability of choosing the kth book (k = 1,2,...,m), and suppose 
the book pile is in the state (i, i,,...,i,,). Then, at the next step, the state 
either remains unchanged, which happens with probability p;, when the top 
book (numbered /,) is chosen, or else changes to one of the m — 1 states 
of the form (i,,7%,,...), which happens with probability p;, when a book 
other than the top book is chosen. Thus we are dealing with a Markov chain, 
with states described by the permutations (i,, i.,...,i,,) and the indicated 
transition probabilities. 

For example, if m = 2, there are only two states e, = (1, 2) and «, = 
(2, 1), and the transition probabilities are 


Pu = Pu = Pi Piz = P22 = Pr 
The corresponding transition probability matrix is 
Pi P2 


Pi. Pe 
The “two-step transition probabilities’’ are 


P= 


Pu(2) = Po (2) = Pip; + Pip2 = Pil(Pr + Po) =. Dy, 
Pr2(2) = Poo(2) = Pipe + Pops = P2(Pi + Po) = Po. 


* Suitably generalized to the case of infinite matrices, if there are infinitely many states 
a 
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Hence P? = P, and more generally P” = P. Given any initial probability 
distribution p!, p§, we have 


p(n) = pipi(n) + popei(n) = pi(pi + pe) = Pr 
p(n) = prpy2(n) + pepe.(n) = pop? + Pz) = Pe- 


Example 2 (The optimal choice problem). Returning to Example 2, p. 28, 
concerning the choice of the best object among m objects all of different 
quality, let «, (k = 1,2,...,m) be the state characterized by the fact that 
the kth inspected object is the best of the first k objects inspected, and let 
Em+1 be the state characterized by the fact that the best of all m objects has 
already been examined and rejected. As the m objects are examined one by 
one at random, there are various times at which the last object examined 
turns out to be better than all previous objects examined. Denote these times 
in order of occurrence by t=0,1,...,¥v, with t= 0 corresponding to 
inspection of the first object and t = v being the time at which the best of all 
m objects is examined (v = 0 if the best object is examined first). Imagine a 
system with possible states ¢),..., &ms mii, and let &(t) be the state of the 
system at the time ¢, so that in particular ¢(0) =«,. To make the “random 
process” &(0) > &(1) ~ &(2)>--: into a Markov chain, we must define 
E(n) for n > v. This is done by the simple artifice of setting &(m) = €y4, for 
alln > v. 

The transition probabilities of this Markov chain are easily found. 
Obviously Pinit.m4i1 = 1 and p, =0 if i> j, 7 < m. To calculate p,,; for 
i<j <m, we write (7.2) in the form 


P(E,E;) 

P(E;) © 
in terms of the events E; = {&(m) = e,} and E; = {E(m + 1) = «,}. Clearly, 
P(E;) is the probability that the best object will occupy the last place in a 
randomly selected permutation of j objects, all of different quality. Since the 
total number of distinct permutations of 7 objects is j!, while the number of 
such permutations with a fixed element in the last (jth) place is (j — 1)!, 
we have 


Pi; = PE; | E;) = (7.9) 


(j= 1)! _ 
j} 


1 
PE) = a oe a ee (7.10) 
J 
Similarly, P(£;£;) is the probability that the best object occupies the jth 
place, while a definite object (namely, the second best object) occupies the 
ith place. Clearly, there are (j — 2)! permutations of j objects with fixed 
elements in two places, and hence 


(j — 2)! 1 


ee 
i! (j—-1j 


ae ae (7.11) 
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It follows from (7.7)-(7.11) that 
ae 
Goeth 


As for the transition probabilities p,; ,,,;, they have in effect already been 
calculated in Example 2, p. 28: 


Pi3 = i<j<m. 


i 
ONS eae 2s eee 
m 


Example 3 (One-dimensional random walk). Consider a particle which 
moves randomly along the x-axis, coming to rest only at the pointsx = ..., 
—2, —1,0,1,2,...with integral coordinates. Suppose the particle’s motion 
is such that once at a point i, it jumps at the next step to either the point 
i — 1 or the point i — 1, with probabilities p and gq = 1 — p, respectively.* 
Let E(n) be the particle’s position after n steps. Then the sequence &(0) — 
E(1) — £(2) — --- is a Markov chain with transition probabilities 

p fe 727, 
Pu=\q Wi g=i-kL (7.12) 
0 otherwise. 

In another kind of one-dimensional random walk, the particle comes to 
rest only at the points x = 0,1, 2,... , jumping from the point i to the point 
i + 1 with probability p, and returning to the origin with probability g, = 
| — p;. The corresponding Markov chain has transition probabilities 

Dp, oo pete 
Pu=\G fF f=9, (7.13) 


0 otherwise. 


16. Persistent and Transient States 


Consider a Markov chain with states ¢,, «,,... and transition probabi- 
lities p;;, i, 7 = 1,2,... Suppose the system is initially in the state ¢,. Let 


un = P i(M), 
and let v,, be the probability that the system returns to the initial state «, 


* Thus the particle’s motion is “‘generated’’ by an infinite sequence of Bernoulli trials 
(cf. the example on pp. 63-65, where p = q = 3). 
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for the first time after precisely n steps. Then 
Un, so Uy, ss UyVy_-1 + cote + Uy—101 = UnVos ae E, B58 ns 
(7.14) 
where we set 
Ug — i. Vo == 0 


by definition. To see this, let B, (K=1,...,m) be the event that “‘the 
system returns to ¢, for the first time after k steps,”’ B,,, the event that “the 
system does not return at all to <; during the first n steps,” and A the event 
that ‘“‘the system is in the initial state <, after n steps.’’ Then the events 
B,,..., B,, Bys1 form a full set of mutually exclusive events, and hence, by 
the “total probability formula”’ (3.6), p. 26, 


P(4) = 5 P(A | B)PCB,) (7.15) 


where clearly P(A | B,,41) = 0 and 
P(B,) — Vx P(A | B,) — Uy—ks k == 2 cee sg nN. 


Substituting these values into (7.15), we get (7.14). 
In terms of the generating functions? 


U(z) = paras V(z) = ree. jz =< f, 
k=0 k=0 


we can write (7.14) in the form 
U(z) — u = U(z)V(2), Uy = |, 
which implies 
1 


eee (7.16) 


U(z) 
The quantity 


00 
- vi 2 


n=0 
is the probability that the system sooner or later returns to the original state 
e,. The state ¢; is said to be persistent if v = 1 and transient if v < 1. 


THEOREM 7.1. The state ¢, is persistent if and only if 


co co 
> un = > pin) = ©. (7.17) 
n=0 n=0 

5 Although the numbers up, Uy, Us, ... do not correspond to a probability distribution 


cO 
as on p. 70 (in fact, we will consider the case where = U; = ©), we continue to call U(z) 
k=0 


0 
a “generating function.” The convergence of the series >} u,z* for |z| < 1 follows by com- 


k=0 
parison with the geometric series, since |ux| < 1 for every k. 


SEC. 16 MARKOV CHAINS 89 
Proof. To say that ¢, is persistent means that 


tae Sy, so lim He) =I, 
n=0 271 
or equivalently, 


lim U(z) = lim ————— =-00 
z—>1 z71 1 — V(z) 
Suppose 
SU Se (7.18) 
n=0 


Then, since the u,, are all nonnegative, 


N 12.6) 
> u, < lim U(z) < > 4, 
n=0 


271 n=0 


for every N, and hence, taking the limit as N — 00, we have 


lim U(2) = 34, 

21 n=0 
In other words, U(z) approaches a finite limit as z > 1 if and only if 
(7.18) holds. Equivalently, U(z) — oo as z—>1, ie., ¢, is persistent, if 
and only if (7.17) holds. J 


THEOREM 7.2. If the initial state <, is persistent, then with probability 
| the system returns infinitely often to ¢, as the number of steps n — o. 
If <, is transient, then with probability 1 the system returns to <, only 
finitely often, i.e., after a certain number of steps the system never again 
returns to ¢;. 


Proof. Suppose the system first returns to ¢, after v, steps, returns 
a second time to ¢; after v, steps, and so on. If there are fewer than k 
returns to ¢; as m — 00, we Set v, = oo. Then the event {v, < oo} means 
that there are at least k returns to ¢,, and the probability of the system 
returning to ¢, at least once is just 


Piy, =< © = #, 


If the event {v, < 00} occurs, the system returns to its initial state ¢«, 
after v, steps, and its subsequent behavior is the same as if it just started 
its motion in ¢,. It follows that 


Bfiy 25001 < 00} = 0. 


Clearly v; = 00 implies v, = 00, and hence vz, < oo implies v, < oo. 
Therefore 


P {v, < 0} = P {vy < c0|v, < co}P {vy < 00} = v’, 
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and similarly, 
Ph, <= Oy = OC) ae PO Se. 


If c, is transient, then v < 1 and hence 
ie.6) co 
Pi, <= oO) 353 f= 
k=1 k=1 


Therefore, by the first Borel-Cantelli lemma (Theorem 2.5, p. 21), with 
probability 1 only finitely many of the events {v, < 00} occur, 1.e., with 
probability 1 the system returns to the state ¢, only finitely often. This 
proves the second assertion in the statement of the theorem. 

On the other hand, if ¢,; is persistent, then v = 1, which implies 


P {y< ysl 
for every k. Let x be the number of times the system returns to its initial 
state ¢, asm — oo. Then obviously the events {x > k} and {v, < oo} are 
equivalent, so that if P {v, < oo} = 1 for every k, then x exceeds any 
preassigned integer k with probability 1. But then 

P {y=} seal, 
which proves the first assertion. J 

A state e, is said to be accessible from a state ¢, if the probability of the 


system going from ¢; to ¢; in some number of steps is positive, i.e., if 
Pij(M) > 0 for some M. 


THEOREM 7.3. Jf a state ¢, is accessible from a persistent state ¢,, 
then «, is in turn accessible from ¢; and «; is itself persistent. 


Proof. Suppose ¢, is not accessible from ¢;. Then the system will go 
from ¢, to ¢, with positive probability p,,() = « > 0 for some number 
of steps M, after which the system cannot return to ¢,. But then the 
probability of the system eventually returning to ¢; cannot exceed 1 — a, 
contrary to the assumption that ¢, is persistent. Hence ¢, must be acces- 
sible from ¢,, i.e., p;,(V) = B > 0 for some N. It follows from (7.8) that 

P(n + M + N)= P(M)P(n)P(N) = PW)P@)P(M), 
and hence 


Piuln + M+ N)> pij(M)pis(pi(N) = «8p;;(), 
Pin + M+ N)> py(N)pil™pis(N) = «Bp;,(n). 


These inequalities show that the series 


> Pal) Sp ji() 
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either both converge or both diverge. But 


¥ Pun) ae OS 


by Theorem 7.1, since ¢, is persistent. Therefore 


2 Pin) =e 
i.e., ¢, is also persistent (again by Theorem 7.1). Jj 


CoROLLARY. Jf a Markov chain has only a finite number of states, 
each accessible from every other state, then the states are all persistent. 


Proof. Since there are only a finite number of states, the system 
must return to at least one of them infinitely often as nm — oo. Hence at 
least one of the states, say ¢,, is persistent. But all the other states are 
accessible from ¢,. It follows from Theorem 7.3 that all the states are 
persistent. fj 


Example 1. In the book pile problem (Example 1, p. 85), if every book 
is chosen with positive probability, i.e., if p, > 0 for alli=1,...,m, then 
obviously every state is accessible from every other state. In this case, all 
m! distinct states (i,,...,i,,) are persistent. If p; = 0 for some /, then all 
states of the form (i,,...,i,,) where i, = i (the ith book lies on top of the 
pile) are transient, since at the very first step a book with a number / different 
from i will be chosen, and then the book numbered i, which can never be 
chosen from the pile, will steadily work its way downward. 


Example 2. In the optimal choice problem (Example 2, p. 86), it is 
obvious that after no more than m steps (m is the total number of objects), 
the system will arrive at the state ¢,,,,, where it will remain forever. Hence 
all the states except ¢,,,, are transient. 


Example 3. Consider the one-dimensional random walk with transition 
probabilities (7.12). Clearly, every state (i.e., every position of the particle) 
is accessible from every other state, and moreover® 


0 if k =2n + 1, 
oo Crs. af b= 2h. 
Using Stirling’s formula (see p. 10), we have 
Cay +. a Jeon One" . 1 


n a Sent = == (4 . 
Pq (nb? P4q (ian n"e-) P4q rar Pq) 


® Cf. formula (5.2), p. 55. 


2n nin 
& 
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for large n, where 
4pqg = (p +g)? — (p—Qt=1—-P-g* <I! 
(the equality holds only for p = q = 3). Therefore 


/ n 
Pii(2n) ~ Jan (4pq) 
for large n, and hence the series 


» P;(2n), 5 ~ (4pq)” 
n=0 n=0 Jan 


either both converge or both diverge. Suppose p # q, so that 4pg < 1. Then 


> Pii(2n) << 0, 
n=0 


and hence every state is transient. It is intuitively clear that if p > q (say), 
then the particle will gradually work its way out along the x axis in the posi- 
tive direction, and sooner or later permanently abandon any given state 7. 
However, if p = q = 4, we have 


> Pi(2n) ae Es . 


and the particle will return to each state infinitely often, a fact apparent from 
the symmetry of the problem in this case. 


Example 4. Next consider the one-dimensional random walk with 
transition probabilities (7.13). Obviously, if0 <p;< 1foralli=0,1,..., 
every state is accessible from every other state, and hence the states are either 
all persistent or all transient. Suppose the system is initially in the state i = 0. 
Then the probability that it does not return to the state i = 0 after n steps 
equals the product Pop, ** Pn—1, the probability of the system making the 
consecutive transitions 0 + 1—--:—n. It is easy to see that the proba- 
bility that the system never returns to its initial state i = 0 as n > oo equals 
the infinite product 

10.6) 
n=0 nN 00 
If this infinite product converges to zero, i.e., if 
lim PoP *** Pn = 9, 
n— co 
then the state i= 0 is persistent, and hence so are all the other states. 
Otherwise, the probability of return to the initial state is 
v.== 1 lim pop, ‘+: p, <1. (7.19) 
no 
Then the state i = 0 is transient, and hence so are all the other states. 
We can arrive at the same result somewhat differently by direct calculation 
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of the probability v, that the particle first returns to its initial state i = 0 in 
precisely m steps. Obviously, v,, is just the probability of the particle making 
the consecutive transitions 0 — 1 ~---—>n — 1 in the first n — 1 steps and 
then returning to the state i = 0 in the nth step. Therefore, since the transi- 
tion i — | ~ihas probability p,_,, 

—— Po» 

"a PoP: *** Pro 1 aa Pi soins 2, 3; — 
By definition, the probability of eventually returning to the initial state i = 0 
is 


Therefore 
D= 1 — po + pol — py) + Popi(l — poe) + °° = 1 —lim pop, -:: p,, 
in keeping with (7.19). 


17. Limiting Probabilities. Stationary Distributions 


As before, let p;(m) be the probability of the system occupying the state 
¢; alter m steps. Then, under certain conditions, the numbers p,(n), j = 1, 


2,... approach definite limits as n > oo: 
THEOREM 7.4. Given a Markov chain with a finite number of states 
21,---5 £m, each accessible from every other state, suppose 
min p;,(N) = 3 > 0 (7.20) 
id 


for some N.’ Then 

lim p(n) —<_ D5. 

n— CO 
where the numbers p;, j =1,...,m, called the limiting probabilities,® 
do not depend on the initial probability distribution and satisfy the in- 
equalities 

max |p,(n) — p*| < Ce", |p,(n) — pt] < Ce-P* (7.21) 

for suitable positive constants C and D. 


Proof. Let 
r(n) = min p,,(n), » Rn) = max p,,(n). 


‘In other words, suppose the probability of the system going from any state ¢,; to any 
other state ¢; in some (fixed) number of steps N is positive. 7 

* Clearly, the numbers p* are nonnegative and have the sum 1 (why ?). Hence they are 
candidates for the probabilities of a discrete probability distribution, as implicit in the 
term “limiting probabilities.” 
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Then 
r(n-+ 1) = min p,(n + 1) = min ¥ paps,(n) > min ¥ par,(n) = r4(n), 
Z 1 k= i k= 
Rjn+1)= max Pin + 1) = max >, PucPril) < max > PieR,(n) = R,(n), 
and hence z i 


ri(1) < 192) <0 <M) < ++ < Rit) < ++ '< RQ) < RD). 
Let N be the same as in (7.20). Then, for arbitrary states ¢, and é,, 


> pal) = Ppl) eB 
Therefore a 
Saal ) —¥ pul ) 
= 2"TPa(N) — Pee(N)] + 2 [Pax(N) — Dox(NY] = 0, 


where the sum 5+ ranges over all k such that p,,(N) — pg,(N) > 0 and 
>— ranges over all k such that py,(N) — pg,.(N) < 0. Clearly, (7.20) 
implies 

man e.: [Pox(N) — Pox(N)] = 4d < 1, 


for some positive number d. 
Next we estimate the differences R,(n) —r,;(n) and R,(n + N) — 
r(n + N): 
RAN) — rN) = max p, WN) — ae pe,(N) 


= [Pa N) — peN)] 
- < ee 2" Puch) = Pox(N)] = d, 


R,n+ N)—r(n+N)= es [Pan + N) — pan + N)] 
= max 5 [Pax(N) ~ Pae()IPus(™) 
< max S* [PaaN) — pse(NIR (0) 

+E [palN) — pea(NDIr(0) 


= max S* paa(N) — Poa NIRA) — rem] 


= d[R,(n) — r,(n)]. 
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It follows that 


R,(KN) = r,(KN) < a’, Ks ‘: Z. oar (725 
But, as already noted, the sequence r;(n), n=1,2,... is nonde- 
creasing while the sequence R,(n), n = 1,2,... is nonincreasing, and 


moreover r;(n) < R;(n). Hence (7.22) shows that both sequences have 
the same limit 
p; =limr,(n) = lim Rn). 


Moreover, it is clear that 
ip.An) — pil < R(n)—ri(n)<d"", i=1,...,m. (7.23) 


Therefore, given any initial distribution p®, i=1,...,m, we have 


¥ pilen) — PF 


lpAn) — p?| = >, PiPisl) — p; 


m (7.24) 
< > pi[R(n) — r,(n)] = Rn) — r{n) < . d<il. 

i=l 

But then 
lim |p(n) — p;| = 0, 
1.e., 
lim p,(n) = p? 
independently of the initial distribution, as asserted. Choosing 
C= : : D= — = Ind 
d N 


in (7.23) and (7.24), we get (7.21). Jj 


Coro.itary. The limiting probabilities p7, j=1,..., m are a 
solution of the system of linear equations 


OF = ¥ pis (> o eenee (7.25) 
Proof. According to (7.5), 
Pn) = 2 n= VPs 
But this becomes (7.25), after taking the limit asn—> oo. J 


Remark. Given an arbitrary Markov chain with states ¢,, e,,..., let 
p?, i=1,2,... be numbers such that 


pe Opis Pepe I 
a 


and 
PS= > PP J=l,2,... (7.26) 
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Choosing pf, i = 1, 2,... as the initial probability distribution, we calculate 
the probability p,;() of finding the system in the state ¢, after n steps, obtain- 
ing 


pl) = > PiPis = Pi; 


pf<2) = Z pAL)p,; = pa PsPis = Djs 


It follows that 


pin) = 9). 3 = ee (7.27) 
for all n =0, 1, 2,... [p,(O) = p? trivially], i.e., the probabilities p;(n), 
j =1,2,... remain unchanged as the system evolves in time. 


A Markov chain is said to be stationary if the probabilities p;(n), 7 = 1, 
2,...remain unchanged for alln = 0,1,2,..., and then the corresponding 
probability distribution with probabilities (7.27) is also said to be stationary. 
It follows from the corollary and the remark that a probability distribution 
P3.j =1,2,...1is stationary if and only if it satisfies the system of equations 
(7.26). Moreover, if the limiting probabilities 

pj; =lim p,(n) (7.28) 
are the same for every initial distribution, then there is a unique stationary 
distribution with probabilities 


es a 
p20, ji 3 43g 
Hence Theorem 7.4 and its corollary can be paraphrased as follows: Subject 


to the condition (7.20), the limiting probabilities (7.28) exist and are the unique 
solution of the system of linear equations (7.25) satisfying the extra conditions 


™ 
p; > 0, > Pi = 1. 
j=1 
Moreover, they form a stationary distribution for the given Markov chain. 


Example 1. In the book pile problem, it will be recalled from p. 86 that 
when m = 2, the stationary distribution 


Pil”) = Pi; P2(n) =— 5 


is established at the very first step. In the case of arbitrary m, let 
Pligsecsim).Gy.---.5) Genote the probability of the transition from the state 
(i;,...,%m) to the state (j,,...,jm), and assume that the probabilities 
Pi>+++ >» Pm are all positive. Then, as shown on p. 85, 


Pi, if Cis, es Jae a (ix, ly, os Js 


of ee aS cs 
= E 0 otherwise, 
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where the permutation (i,, i,,...) is obtained from (i,,...,i,,) by choosing 
some i, and moving it into the first position. The limiting probabilities 


Pu,..... ;,) are the solution of the system of linear equations 
PG,....i.) = Pix = = ee ee SC (7.29) 
Jyr- ceed, 


where (j;,..../j,,) ranges over the m permutations 


(ji; Jas js: 9. i8 iad (jos J1s J or ens ce eae) (ji, je; «26s ig oe a 


which give (j;,... jm) when j, is moved into the first position. 

After a sufficiently large number of steps, a stationary distribution will be 
virtually established, i.e., the book pile will occupy the states (i,,..., in) 
with virtually unchanging probabilities pj, __,_). Clearly, the probability of 
finding the ith book on top of the pile is then ; 


— x 
a z Fitstgy iss bm)? 
eae | on 


and hence, by (7.29), 


ee, ae egy 
where (ij,...,i,,) ranges over the m permutations 
ee =) AY, so te), is os fm, i) 
which give (i, i,,...,%,,) when i is moved into the first position. But then 


ieee 2 ase ae 
Pi =P; > Pic igk ee potty... Mm, 


i.e., the limiting probability p* of finding the ith book on top of the pile 
is just the probability p,; with which the ith book is chosen. Thus, the more 
often a book is chosen, the greater the probability of its ending up on top 
of the pile (which is hardly surprising). 


Example 2. Consider the one-dimensional random walk with transition 
probabilities (7.12). If p ~4q, then the particle gradually moves further and 
further away from the origin, in the positive direction if p > g and in the 
negative direction if p <q. If p = q, the particle will return infinitely often 
to each state, but for any fixed j, the probability p; () of the particle being 
at the point j approaches 0 as n + co (why?). Hence, in any case, 


lim p(n) = pj = 0 


for every j, but the numbers p;**, j = 1, 2,... cannot be interpreted as the 
limiting probabilities, since they are all zero. In particular, there is no sta- 
tionary distribution. 
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Example 3. Finally, consider the one-dimensional random walk with 
transition probabilities (7.13). Suppose 


km fo, -* *p, = 1 See (7.30) 


so that the states are all transient (see p. 92). Then as n — oo, the particle 
“moves off to infinity” in the positive direction with probability 1, and there 
is obviously no stationary distribution. If there is a stationary distribution, 
it must satisfy the system of equations (7.26), which in the present case take 
the form 


D3 == ey ee J sa i pat eee (7.31) 
It follows from (7.31) that 


PL = PoPos = P2 = PoPoPrs--+> Pn = PoPoPi*** Pn—1>- + 
Clearly a stationary distribution exists if and only if the series 


2 PoP1"** Pn = 1+ Pot PoP + °° (7.32) 
converges.® The stationary distribution is then 
1 
eee = epee are 
1. Po Pea 
pe ae Se 8 eee 
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PROBLEMS 


1. A number from 1 to m is chosen at random, at each of the timest = 1,2,... 
A system is said to be in the state ¢, if no number has yet been chosen, and in the 
state ¢,; if the largest number so far chosen is 7. Show that the random process 
described by this model is a Markov chain. Find the corresponding transition 
probabilities p;; (i, 7 =0,1,..., m). 


i 1 
Ans. | <c eeere, Pixs = Oifi > J, Pig = if <j. 


2. In the preceding problem, which states are persistent and which transient? 


3. Suppose m = 4 in Problem 1. Find the matrix P(2) = ||p;;(2)||, where 
Pi;(2) is the probability that the system will go from state ¢; to state ¢; in 2 steps. 


* Note that (7.32) automatically diverges if (7.30) holds. 
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4. An urn contains a total of N balls, some black and some white. Samples are 
drawn from the urn, m balls at a time (m < N). After drawing each sample, 
the black balls are returned to the urn, while the white balls are replaced by 
black balls and then returned to the urn. If the number of white balls in the urn 
is i, we say that the ‘“‘system”’ is in the state ¢,. Prove that the random process 
described by this model is a Markov chain (imagine that samples are drawn at 
the times t = 1, 2,.. . and that the system has some initial probability distribu- 
tion). Find the corresponding transition probabilities p;; @, 7 =0,1,..., N). 
Which states are persistent and which transient? 


Ans. p33 =O ifi<jorifi>j,j >N—m, 


OE es ae 
oe, (i>, j7 SN — m. 


m 


The state ¢ is persistent, but the others are transient. 


5. In the preceding problems, let N = 8, m = 4, and suppose there are initially 
5 white balls in the urn. What is the probability that no white balls are left after 
2 drawings (of 4 balls each)? 


6. A particle moves randomly along the interval [1, m], coming to rest only at 
the points with coordinates x = 1,...,m. The particle’s motion is described 
by a Markov chain such that 


Piz = - Pmm-1 = ., 
Pi,j+1 = P> Fiat = (= 2...., — 1), 
with all other transition probabilities equal to zero. Which states are persistent 


and which transient ? 


7. In the preceding problem, show that the limiting probabilities defined in 
Theorem 7.4 do not exist. In particular, show that the condition (7.20) does not 
hold for any N. 


Hint. py(n) = 0 if nis odd, while py.(n) = 0 if n is even. 
8. Consider the same kind of random walk as in Problem 6, but now suppose 
the nonzero transition probabilities are 
Pu ps 7; Pmm = Pp; 
Pi,i+1 =P» Pisa l= Tye. 5m); 


permitting the particle to stay at the points x = 1 and x = m. Which states 
are persistent and which transient? Show that the limiting probabilities 
px,..., p*, defined in Theorem 7.4 now exist. 


9. In the preceding problem, calculate the limiting probabilities pf, ... , fs 
Ans. Solving the system of equations 
Pr =4Pi + 4P2> 
PS =PPia +i Gj 
Pm = PPm—1 + PPm 


250s s pa AD, 
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we get 
py 
P3 - (4) Fee Se TE 
Therefore 
eer 
if p = q, while 


Pe eh (2) 

ee 

if p #q (impose the condition that S p; = 1). 
j=1 


10. Two marksmen A and B take turns shooting at a target. It is agreed that A 
will shoot after each hit, while B will shoot after each miss. Suppose A hits the 
target with probability « > 0, while B hits the target with probability 8 > 0, 
and let 7 be the number of shots fired. What is the limiting probability of hitting 
the target asn = «©? 


ae 
l1—a— 68° 


11. Suppose the condition (7.20) holds for a transition probability matrix 
whose column sums (as well as row sums) all equal unity. Find the limiting 
probabilities p?,... , p%,. 


Ans. 


ee. 

Ans... pp.= ++: eS 
12. Suppose m white balls and m black balls are mixed together and divided 
equally between two urns. A ball is then drawn at random from each urn and 
put into the other urn. Suppose this is done 7 times. If the number of white 
balls in a given urn is j, we say that the “‘system”’ is in the state ¢; (the number 
of white balls in the other urn is then m —/). Prove that the limiting prob- 
abilities pF, = eo gee defined in Theorem 7.4 exist, and calculate them. 


Hint. The only nonzero transition probabilities are 


_ 2j(m —j) 2 ite = 
Hee at ae Pin = see Seer iia = eS 
Ans. Solving the system 
P3 = P3-1Ps—a,3 + P§Pis + Phir Pi41,3 (j =0,1,...,m), 


we get p* = (Cj")®p}, and hence 


Cc”)? cm)? 
fee Ee (j =0,1,...,m) 


m 
pas = 
j=0 


(recall Problem 17, p. 12). 
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13. Find the stationary distribution p?, p°, . .. for the Markov chain whose only 
nonzero transition probabilities are 


14. Two gamblers A and B repeatedly play a game such that A’s probability 
of winning is p, while B’s probability of winning is g = 1 — p. Each bet isa 
doJlar, and the total capital of both players is m dollars. Find the probability 
of each player being ruined, given that A’s initial capital is 7 dollars. 


Hint. Let ¢; denote the state in which A has j dollars. Then the situation 
is described by a Markov chain whose only nonzero transition probabilities are 
Poo = 1, Pmm > k; 
Piitt =P» Pij31 = 4 (pt, 5 — Fp. 
Ans. Let p; = lim pjo(n) be the probability of A’s ruin, starting with an 
n— 0 
initial capital of j dollars. Then 
Pi =PP2+9, = Pm—1 = 9Pm-2> 
Ps = (Pi-1 + PPis1 (F=2, 4.8 2) 


(why ?). Solving this system of equations, we get 


ve ee (7.33) 


m 
if p = q (as in Example 3, p. 29), and 
1 — Pig’ 

pot aoe 7.34 

Pi T= (piqy™ See 
if p # g. The probability of B’s ruin is 1 — fj. 
15. In the preceding problem, prove that if p > q, then A’s probability of ruin 
increases if the stakes are doubled. 


16. Prove that a gambler playing against an adversary with unlimited capital 
is certain to be ruined unless his probability of winning in each play of the game 
exceeds $. 


Hint. Let m— © in (7.33) and (7.34). 
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CONTINUOUS MARKOV PROCESSES 


18. Definitions. The Sojourn Time 


Consider a physical system with the following properties, which are 
the exact analogues of those given on p. 83 for a Markov chain, except 
that now the time ¢ varies continuously: 


a) The system can occupy any of a finite or countably infinite number of 
states €,, €9,... 

b) Starting from some initial state at time t = 0, the system changes its 
state randomly at subsequent times. Thus, the evolution of the system 
in time is described by the “random function’’ &(t), equal to the state 
of the system at time ¢.? 

c) At time t = 0, the system occupies the state ¢,; with initial probability 


HPO = c, i= 1,2,.. 


d) Suppose the system is in the state ¢; at any time s. Then the probability 
that the system goes into the state ¢, after a time ¢ is given by 


p=? (he) =2), 9 He Se 


regardless of its behavior before the time s. The numbers p,,(t), called 
the transition probabilities, do not depend on the time s. 


A random process described by this model is called a continuous Markov 


1 Recall footnote 1, p. 83. Note that €(¢) is a random variable for any fixed 1. 
102 
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process? or simply a Markov process (as opposed to a Markov chain, which 
might be called a “discrete Markov process’’). 
Let 


p(t) = P {&(t) = ¢,}, Fe Bay Ks 
be the probability that the system will be in the state ¢, at time ¢. Then, by 
arguments which hardly differ from those given on p. 84, we have 


PAO) = Pp, 1 oes 


pis +t) = 3 PAS) Pxs(t), |e Oe eee (8.2) 
and s 
0) l ft j= i, 8.3) 
ae ee ee S 
AC 2 t) — 2 Pal) Pr), i sf si ‘, Z. cee (8.4) 


for arbitrary s and ¢ [cf. (7.5) and (7.7)]. 


THEOREM 8.1. Given a Markov process in the state < at time t = to, 
let < be the (random) time it takes the process to leave « by going to 
some other state.* Then 


Pie2 t= ot > 0, (8.5) 
where 2 is a nonnegative constant. 
Proof. Clearly P {t > t} is some function of f, say 
aye rt Se, t> 0. 


If t > s, then the process will be in the same state at time 7) + s as at 
time f), and hence its subsequent behavior will be the same as if s = 0. 
In particular, 

P{t>s+t|t>s}= (t) 


is the probability of the event {t > s+ t} given that + > s. It follows 
that 
P{r>st+t}=P{r>s4+t|t>s}P {tr > s}= oo), 


and hence 


o(s + t) = 9(s) p(t) 
or equivalently 
In o(s + t) = In o(s) + In off) 


2 More exactly, a continuous Markov process with stationary transition probabilities. 
where we allude to the fact that the numbers (8.1) do not depend on s (cf. footnote 2, 
p. 84). 

® Here we prefer to talk about states of the process rather than states of the system (as 
in Chap. 7). 
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for arbitrary s and ¢t. Therefore In 9(t) is proportional to ¢ (recall 
footnote 4, p. 40), say 


In p(t) = —At, t>0, (8.6) 


where A is some nonnegative constant (why nonnegative?). But (8.6) 
implies (8.5). Jf 


The parameter A figuring in (8.4) is called the density of the transition 
out of the state «. If A = 0, the process remains forever in «. If A > 0, the 
probability of the process undergoing a change of state in a small time 
internal At is clearly 

1 — o(At) = A At + o(Ad), (8.7) 


where o(Ar) denotes an infinitesimal of higher order than At. 
It follows from (8.5) that 


te 
P {t, <<t < te} = O(t,) — Ot, =e — e™ -| he dt (8.8) 
ty 


for arbitrary nonnegative ¢, and f, (t, < f,). Therefore the random variable 
t, called the sojourn time in state <, has the probability density 


in | ft tS, — 
pAt) = 

0 ae eae 
The distribution corresponding to (8.8) and (8.9) is called the exponential 
distribution, with parameter 4. The mean value Et, i.e., the “‘expected 
sojourn time in state <,” is given by 


Et -| ih Wpare 
0 d 

Example (Radioactive decay). In Example 3, p. 58, we gave a proba- 
bilistic model of the radioactive decay of radium (Ra) into radon (Rn). The 
behavior of each of the my radium atoms is described by a Markov process 
with two states (Ra and Rn) and one possible transition (Ra — Rn). As on 
p. 58, let p(t) be the probability that a radium atom decays into a radon 
atom in time ¢, and &(t) the number of alpha particles emitted in ¢ seconds. 
Then, according to formula (5.7), 


k 
P(E) =k} =e k =0,1,2,..., 


where 


a = E&(t) = nop(t). 
It follows from (8.5) that 


ot) = 1 = o™, iS 8 
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where 2 is the density of the transition Ra — Rn. Recalling (8.7), we see that 
2 is the constant such that the probability of the transition Ra > Rn in a 
small time interval At equals AAt + o(A?). 

The number of (undisintegrated) radium atoms left after time ¢ is clearly 
Ny — &(t), with mean value 


n(t) = E[ny — &(t)] = mo — nop(t) = ne, t > 0. (8.10) 
Let T be the half-life of radium, i.e., the amount of time required for half the 
radium to disappear (on the average). Then 
| 
a 8 5 Mo» (8.11) 
and hence, comparing (8.10) and (8.11), we find that 7 is related to the 
density of the transition Ra — Rn by the formula 


_in2 
2 


+ 


19. The Kolmogorov Equations 


Next we find differential equations satisfied by the transition probabilities 
of a Markov process: 


THEOREM 8.2. Given a Markov process with a finite number of states, 
suppose the transition probabilities p,;(t) are such that* 


1 — p,,(At) = 4, At + o(At), bee B52 oes 


P(At) = rj; At + o(At), J# i, i,jJ= i; te 
and let 


(8.12) 


Ta Mg; i= ‘; vA coe (8.13) 


Then the transition probabilities satisfy two systems of linear differential 
equations, for forward Kolmogorov equations? 


Pi;(t) oo x Pix) Anis ly J ices kr, es i (8.14) 
k 
and the backward Kolmogorov equations 
Dit) = 2 AinPrs(D, 3F — = a ee (8.15) 


subject to the initial conditions (8.3). 


4 We might call (; the “density of the transition out of the state ¢;,” and 2,,; the “density 
of the transition from the state ¢; to the state ¢;.” 
®* The prime denotes differentiation with respect to ¢. 
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Proof. It follows from (8.4) that 
P(t + At) = p3 Pix(t)P,(At) = 2 Pix( At) Dy (1). 
Hence, using (8.12) and (8.13), we have 


Di(t + At) — p;,(t) o(At) | 
wp DS pet] Des =>}, (t). 
te © pald| i a >| gr a Py jt) 
Both sums have definite limits as At — 0. In fact, 
At 
lim > pa(t) a 2 é =O: (8.16) 
At>0 x At k 
At 
lim & | du + at "]pait) = & Darr (8.17) 
At70 & At ke 
Therefore 
At70 At 
also exists, and equals (8.16) and (8.17). J 
Remark 1. It follows from (8.12) and the condition 
> P;,fAt) = 1 
“ 
that 
So (8.18) 


j#i 

Remark 2. The Kolmogorov equations hold not only in the case of a 
finite number of states, but also in the case of a countably infinite number of 
states €,, €2,... if we make certain additional assumptions. In fact, suppose 
the error terms o(Afr) in (8.12) are such that 


o(At) oe 
At 


uniformly in all i and j. Then the forward equations (8.14) hold if for any 
fixed 7, there is a constant C < oo such that 


hg eo, oe ey eer 


0 as At—O 


while the backward equations (8.15) hold if the series (8.18) converges. 


Example 1 (The Poisson process). As in Example 4, p. 73, consider 
a “random flow of events’’ with density A, and let &(t) be the number of events 
which occur in time ¢. Then &(ft) is called a Poisson process. Clearly &(t) is a 
Markov process, whose states can be described by the integers 0, 1, 2,... 
Moreover, &(t) can only leave the state i by going into the state 7+ 1. 
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Therefore the transition densities ,; are just 
OS pas 
Poe ga eee, 
hy = —A, 
where we use (8.13) and (8.18). 


The transition probabilities p;,(t) of the Poisson process &(t) clearly 
satisfy the condition 


Pii(®) = Po,s—s(t) 
Pi(t) = Pot), j=0,1,2,... 
Then the forward Kolmogorov equations take the form 
Po(t) = —Apo(t), 
DAt) = Apyalt) — Apt), =f =1,2,... 
Introducing the new functions 
thee 2d). At 2,» os 


(why ?). Let 


we find that 
Folt) = Af(t) + &*po(t) = AK(t) — rep o(t) = 0, 
F(t) = FD) + eA) 
in = If,(t) 5 re Psat) — rep (t) are Af ;-1(t), J ee e a ee 
fo(9) = I, 
fO%=0, j=1,2,..., 
because of (8.3). But the solution of the system of differential equations 
fo(t) = 9, 
Fit) +e Mf j;_s(t), j= oe oe 


subject to the initial conditions (8.19), is obviously 


where 


(8.19) 


1S ep es a 


Returning to the original functions p,(t) = cane we nina that 


ty? 
pA) = a a. 


or equivalently 
Pig) = 


just as on p. 75. 


j 
poe a ee 
J: 
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Example 2 (A service system with exponential holding times). Consider 
a random flow of service calls arriving at a server, where the incoming 
“traffic” is of the Poisson type described in Example 1, with density 4. Thus 
AAt + o(At) is the probability that at least one call arrives in a small time 
interval At. Suppose it takes a random time t to service each incoming call, 
where t has an exponential distribution with parameter wu: 


Pit>th=e* - (8.20) 


(the case of “exponential holding times’’). Then the service system has two 
states, a state c, if the server is ‘‘free’’ and a state ¢, if the server is “busy.” 
It will be assumed that a call is rejected (and is no longer a candidate for 
service) if it arrives when the server is busy. 

Suppose the system is in the state ¢) at time f). Then its subsequent 
behavior does not depend on its previous history, since the calls arrive 
independently. The probability p),(At) of the system going from the state 
Ey to the state c, during a small time interval At is just the probability AAt + 
o(At) of at least one call arriving during At. Hence the density of the 
transition from ¢, to ¢, equals 4. On the other hand, suppose the system is 
in the state c, at time ¢,;. Then the probability p,9(t) of the system going from 
the state <, to the state e) after a time ¢ is just the probability that service 
will fail to last another ¢ seconds.* Suppose that at the time f,, service has 
already been in progress for exactly s seconds. Then 


P{t>s+t} 
2 ee ey eee i+ > = 
Pio(t) west & } Pit>s} 
Using (8.20), we find that 
e uls+t) 
Pio(t) = 1 — oe 1—e™, (8.21) 
e 


regardless of the time s, i.e., regardless of the system’s behavior before the 
time ¢,.’ Hence the system can be described by a Markov process, with two 
states €) and ¢. 
The transition probabilities of this Markov process obviously satisfy the 
conditions 
Polt)=1— Polt), Pro(t)= 1 — put). (8.22) 
Moreover, 
doo = —A,; Aor = A, 


Mio = U fag es 
® For simplicity, we choose seconds as the time units. 


“It is important to note that this is true only for exponential holding times (see W. 
Feller, op. cit., p. 458). 
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where we use the fact that 


Pi(At) = 1 — eo *** = wAt + o(Ad). 


Hence in this case the forward Kolmogorov equations (8.14) become?® 


Poolt) = AooPoo(t) + AroPorl(t) = —APoolt) + uLL — Pool(t)], 
Pir(t) = Agr Prol(t) + ArwPu(t) = ALL — Poolt)] — vPir(), 


l09 


1.€., 
Poot?) + (A + 2)Poo(t) = pv, (8.23) 
Pir(t) + (A + p)pii(t) = A. 
Solving (8.23) subject to the initial conditions 
Poo(9) = pu(O) = 1, 
we get 
U. —(A+u)t U 
pat) = (1 - Je See Sora 
r r 
- x : (8.24) 
fas 11 —(A+u)t 
Pir(t) ( ae . 


20. More on Limiting Probabilities. Erlang’s Formula 


We now prove the continuous analogue of Theorem 7.4: 


THEOREM 8.3. Let &(t) be a Markov process with a finite number of 
states, €,,... 5 €m, each accessible from every other state. Then 


= pt) = pj; 


where p,(t) is the probability of &(t) being in the state <, at time t. The 
numbers Dp; , j=1,...,m, called the limiting probabilities, do not 
depend on the initial probability distribution and satisfy the inequalities 


maxipg)—pia Ce, |pft)— pil < Ce?’ (8.25) 


for suitable positive constants C and D. 


Proof. The proof is virtually the same as that of Theorem 7.4 for 
Markov chains, once we verify that the continuous analogue of the 
condition (7.20), p. 93 is automatically satisfied. In fact, we now have 


min p(t) = 6) > 0 (8.26) 


§ Because of (8.22), there is no need to write equations for po,(t) and p;i,(t). 
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for all t > 0. To show this, we first observe that p,,(t) is positive for 
sufficiently small ¢, being a continuous function (why?) satisfying the 
condition p,,(0) = 1. But, because of (8.4), 

PilS + t) > pils)pi(t) 


for arbitrary s and t, and hence p,,(t) 1s positive for all tf. 
To show that p,,(t), i A / is also positive for all ¢, thereby proving 
(8.26) and the theorem, we note that 


Pis(s) > 0 
for some 5, since ¢,; is accessible from ¢,. But 
Pit) > puWpiy(t — 4), u<tl, 
again by (8.4), where, as just shown, p,,(t — u) is always positive. 
Hence it suffices to show that p;,(u) > 0 for some u < ¢. Consider a 
Markov chain with the same states ¢,,...,¢,, and transition proba- 
bilities 
S 
Pig = Pis(7}> 


where n is an integer such that 


Since 


the state <, is accessible from ¢,. But it is easy to see that ¢, is accessible 
from ¢, not only inn steps, but also in a number of steps ny no greater 
than the total number of states m (think this through). Therefore 


where 


The limiting probabilities p*, 7 = 1,..., m forma stationary distribution 
in the same sense as on p. 96. More exactly, if we choose the initial distri- 
bution 


then 
p() = FF. ie gees. 
i.e., the probability of the system being in the state e,; remains unchanged 
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for all t > 0. In fact, taking the limit as s — 0 in (8.2), we get 
pi = Dd Pipdt), j=i,...,m. (8.27) 


But the right-hand side is just p,(t), as we see by choosing s = 0 in (8.2). 
Suppose the transition probabilities satisfy the conditions (8.12). Then 
differentiating (8.27) and setting t = 0, we find that 


> Pir = 9, ho desis ca M, (8.28) 
where 2,; is the density of the transition from the state ¢; to the state ¢;. 


Example (A service system with m servers). Consider a service system 
which can handle up to m incoming calls at once, i.e., suppose there are m 
servers and an incoming call can be handled if at least one server is free. As 
in Example 2, p. 108, we assume that the incoming traffic is of the Poisson 
type with density 4, and that the time it takes each server to service a call 
is exponentially distributed with parameter py (this is again a case of “expo- 
nential holding times”). Moreover, it will be assumed that a call is rejected 
(and is no longer a candidate for service) if it arrives when all m servers are 
busy, and that the “holding times’’ of the m servers are independent random 
variables. 

If precisely j servers are busy, we say that the service system is in the 
state ¢; (j= 0,1,...,m). In particular, ¢) means that the whole system is 
free and ¢,, that the system is completely busy. For almost the same reasons 
as on p. 108, the evolution of the system in time from state to state is described 
by a Markov process. The only nonzero transition probabilities of this 
process are 


hoo = —k, Agi =A, Amm = — MEY, 
Ayia = Je. Agg = —-AtiW), A =A G=l...,m—D). 
In fact, suppose the system is in the state <,. Then a transition from ¢; to 
¢,,, takes place if a single call arrives, which happens in a small time interval 


At with probability AAt + o(At).® Moreover, the probability that none of the 
j busy servers becomes free in time At is just 


[1 — pAt + o(At)}’, 


(8.29) 


since the holding times are independent, and hence the probability of at 
least one server becoming free in time Af equals 


1 — [1 — pAt + o(At)]}i = juAt + o(At). 


® For small Af, this is also the probability of at least one call arriving in At. 
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But for small Af, this is also the probability of a single server becoming free 
in time At, i.e., of a transition from ¢, to ¢,_,. The transitions to new states 
other than ¢;_, or ¢;,; have small probabilities of order o(At). These con- 
siderations, together with (8.12) and the formula 


> Ay = 0 
j 


implied by (8.12) and (8.13), lead at once to (8.29). 

In the case m = 1, it is clear from the formulas (8.24) that the transition 
probabilities p,,(t) approach their limiting values “exponentially fast’’ as 
t—» oo. It follows from the general formula (8.25) that the same is true in the 
case m > 1 (more than | server). To find these limiting probabilities pj, 
we use (8.28) and (8.29), obtaining the following system of linear equations: 


Po = BPs 
(A+ ju)p} =Ap3. + Gt Dep = =1,...,m—D, 
ADn—1 = MLD 


Solving this system, we get 


1 (AV 
pt == (") 26. j= 0,1, »™ 
J+\v. 
Using the “normalization condition”’ 
2.2; = 1 
j=0 


to determine p*, we finally obtain Erlang’s formula 


*— + ., Se! eres 8.30 
DP; m 4 />\i J ( ) 
j j ! ( 
j=0]-\U 
for the limiting probabilities. 
PROBLEMS 


1. Suppose each alpha particle emitted by a sample of radium has probability 

p of being recorded by a Geiger counter. What is the probability of exactly n 

particles being recorded in ¢ seconds? 

(pt)" 
n!} 


Ans e~\Pt where 2 is the same as in the example on p. 104. 
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2. A man has two telephones on his desk, one receiving calls with density A,, 
the other with density 2,.1° What is the probability of exactly n calls being 
received in ¢ seconds? 


Hint. Recall Problem 9, p. 81. Neglect the effect of the lines being found 
busy. 
[Ay + Ag)t]" 
n!} 
3. Given a Poisson process with density 2, let &(t) be the number of events 
occurring in time ¢. Find the correlation coefficient of the random variables 
E(t) and &(¢ + +t), where + > 0. 


t 
Ans. ? 
C+ 


4. Show that (8.24) leads to Erlang’s formula (8.30) for m = 1. 


5. The arrival of customers at the complaint desk of a department store is 
described by a Poisson process with density 4. Suppose each clerk takes a 
random time t to handle a complaint, where + has an exponential distribution 
with parameter »., and suppose a customer leaves whenever he finds all the clerks 
busy. How many clerks are needed to make the probability of customers 
leaving unserved less than 0.015 if A = uw? 


Hint. Use Erlang’s formula (8.30). 
Ans. Four. 


6. A single repairman services m automatic machines, which normally do not 
require his attention. Each machine has probability 4Ar + o(Ar) of breaking 
down in a small time interval At. The time required to repair each machine is 
exponentially distributed with parameter y. Find the limiting probability of 
exactly 7 machines being out of order. 


e—ArtAg) t. 


Hint. Solve the system of equations 


mpy = Upy, 
[(m —j)d + vlpf = (m —j + Ih, + UPF +1 
UPm = Pn_1: 
m! r\ 
Ans. pt = ay (=)or. eres 


™m 
where p* is determined from the condition > p* = 
j=0 
Comment. Note the similarity between this result and formula (8.30). 
7. In the preceding problem, find the average number of machines awaiting 
the repairman’s attention. 


10 Tt is assumed that the incoming calls on each line form a Poisson process. 
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8. Solve Problem 6 for the case of r repairmen, where 1 <r < m. 


9. An electric power line serves m identical machines, each operating inde- 
pendently of the others. Suppose that in a small interval of time At each machine 
has probability XAt + o(At) of being turned on and probability At + o(At) 
of being turned off. Find the limiting probability p¥ of exactly 7 machines 
being on. 


Hint. Solve the system of equations 
my = UPy> 
[(m —j)a + julpf =(m —j + Wrapf, + G+ Depa, 
mup*, = \pn_r- 


- U. sar x j 
Ans. p3 = Cj Wier Sar if = Oy he. a 9 2M. 


10. Show that the answer to the preceding problem is just what one would 
expect by an elementary argument if 4 = u. 


Appendix | 


INFORMATION THEORY 


Given a random experiment with N equiprobable outcomes A,,..., Ay, 
how much “information”’’ is conveyed on the average by a message 4% 
telling us which of the outcomes A,,..., Ay has actually occurred? As a 
reasonable measure of this information, we might take the average length of 
the message -@, provided .@ is written in an “economical way.’’ For 
example, suppose we use a “‘binary code,” representing each of the possible 
outcomes A,,..., Ay by a “code word”’ of length /, 1.e., by a sequence 


A ie8i 


where each “‘digit”’ a, is either a 0 or a 1. Obviously there are 2’ such words 
(all of the same length /), and hence to be capable of uniquely designating the 
N possible outcomes, we must choose a value of / such that 


Nm 2. (1) 
The smallest value of / satisfying (1) is just the integer such that 
0< /— log, N < 1. 
This being the case, the quantity 
I = log, N (2) 


is Clearly a reasonable definition of the average amount of information in the 
message .@ (measured in binary units or “‘bits’’). 
More generally, suppose the outcomes A,,..., Ay have different proba- 
bilities | 7 
Pi = P(Aj), ... » Pw = P(Ay). (3) 
115 
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Then it is clear that being told about a rare outcome conveys more informa- 
tion than being told about a likely outcome.’ To take this into account, we 
repeat the experiment n times, where n is very large, and send a new message 
M' conveying the result of the whole series of n trials. Each outcome is now 
a sequence 


eee (4) 


where A, is the outcome occurring at the kth trial. Of the N” possible out- 
comes a the whole series of trials, it is overwhelmingly likely that the outcome 
will belong to a much smaller set containing only 


n! 


BP cree (5) 


outcomes, where 
Nn, = Npy,...,Ny = MPN, Nt + ay =O. 


In fact, let n, = n(A,) be the number of occurrences of the event A; in N 
trials. Then 


nN; 

ee 

n 
by the law of large numbers (see Sec. 12), and hence n; ~ np;,. To get (5), 
we merely replace ~ by = and invoke Theorem 1.4, p. 7. We emphasize 
that this is a plausibility argument and not a rigorous proof,” but the basic 
idea is perfectly sound. 

Continuing in this vein, we argue that only a negligibly small amount of 

information is lost on the average if we neglect all but the set of N,, highly 
likely outcomes of the form (4), all with the same probability 


PA(,) *** P(A,,) = Pit ++ Pa 


This brings us back to the case of equiprobable outcomes, and suggests 
defining the average amount of information conveyed by the message @' as 


Hence, dividing by the number of trials, we find that the average amount of 
information in the original message -Z is just 


pa eee (6) 
: | 


1 In particular, no information at all is conveyed by being told that the sure event has 
occurred, because we already know what the message will be! 

2 Among other missing details, we note that the numbers 7,,..., my are in general not 
all integers, as assumed in (4). 
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To calculate (6), we apply Stirling’s formula (see p. 10) to the expression 
(5), obtaining 


7 J2nn ne” 
/2nn, iE ear J2nny tye 
and hence 
In N, ~nInn — np,In(np,) — - +: — npy In (npn) 
= ninn — (np, +--+ + npy) Inn — np, in p, —+:+-:—npy In py 


N 
= —n) p;In p; 
i=1 
in terms of the natural logarithm 
‘x= toe, J, 
or equivalently 


N F 
log, N,~—n = p; logs p; (7) 


in terms of the logarithm to the base 2. Changing ~ to = and substituting 
(7) into (6), we get Shannon’s formula 


N 
[= — = Dp; 10s P; (8) 


i=1 


for the average amount of information in a message -@ telling which of the 


N outcomes A,,..., Ay with probabilities (3) has occurred. Note that (8) 
reduces to (2) if the outcomes are equiprobable, since then 
a oe 
1 ee 


Example 1 (Average time of psychological reaction). One of N lamps is 
illuminated at random, where p, is the probability of the ith lamp being turned 
on, and an observer is asked to point out the lamp which is lit. In a long 
series of independent trials it turns out* that the average time required to give 
the correct answer is proportional to the quantity (8) rather than to the 
number of lamps N, as might have been expected. 


We can interpret the quantity (8) not only as the average amount of 
information conveyed by the message .@, but also the average amount of 
‘“‘uncertainty’’ residing in the given random experiment, and hence as a 
measure of the randomness of the experiment. Receiving the message 
reduces the uncertainty of the outcome of the experiment to zero, since the 


3 See A. M. Yaglom and I. M. Yaglom, Wahrscheinlichkeit und Information, second 
edition, VEB Deutscher Verlag der Wissenschaften, Berlin (1965), p. 67. 
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message tells us the result of the experiment with complete certainty. More 
generally, we might ask for the amount of information about one “full set’’ 
of mutually exclusive events A,,..., Ay conveyed by being told which of 
a related full set of mutually exclusive events B,,..., By, has occurred. 
Suppose the two sets of events have probabilities P(A,),..., P(Ay) and 
P(B,),..., P( By), where P(A,) +--+ + P(Ay) = 1, 


P(B,) + +++ + Py) = I- 


Moreover, let P(A,B;) be the probability that both events A, and B; occur, 
while P(A, | B,) is the probability of A; occurring if B; is known to have 
occurred. Then 


N 
Lap; ae ae | B;) log, P(A; | B;) 


is the amount of uncertainty about the events A,,..., Ay remaining after 
B, is known to have occurred, and hence 
N’ N’ N 
| mearmean > PB, Laz; ee 2 2, POBDE, | B;) log, P(A; | B;) 
j= j= lie 
P(A;B,) 

2, (A;B;) loge P(B,) (9) 
is the average amount of uncertainty about A,,..., Ay remaining after it is 
known which of the events B,,..., By, has occurred. Let J,, be the in- 
formation about the events A,,..., Ay conveyed by knowledge of which of 
the events B,,..., By, has occurred. Then clearly* 

lap = 1, ae (10) 


where 


N N N’ 
Iq = — > P(A) log, P(A) = —  SP(A,B,) log: P(A) — (1) 


is the quantity previously denoted by J (justify the last step). Combining (9) 
and (11), we finally get 
P(A,B; 
Ly = 3 P42) 
t,3 


P(A,)P(B;) ve 


Example 2 (Weather prediction). During a certain season it rains about 
once every five days, the weather being fair the rest of the time. Every night 
a prediction is made of the next day’s weather. Suppose a prediction of rain 
is wrong about half the time, while a prediction of fair weather is wrong 


4In words, (10) says that “the information in the message’’ equals “the uncertainty 
before the message is received” minus “the uncertainty after the message is received.” 


APP. I INFORMATION THEORY  I19 


only about one time out of ten. How much information about the weather 
is conveyed on the average by the predictions? 


Solution. Let A, denote rain, A, fair weather, B, a prediction of rain 
anc 5, a prediction of fair weather. Then, to a good approximation, 


1 4 
P(A,) == P(A,) = — 
(4)=5, P)==, 

I 1 

P(A, |B) =5, P(A, | B) = 7. 


Moreover, since 


P(A,) = P(A, | B,)P(By) + P(A, | B)P(Ba), 


we have 
oo) i ee 
a ee eee a 
and hence 
l 3 
2 a P(Bs) = 7 
1 
P(A, B,) = P(A, | B,)P(B,) = 8 ) 
3 
P(A,B,) = P(A, | B,)P(B:) = AO ) 


1 
P(A,B,) = [1 — P(A, | B,)|P(B,) = 8 : 


a 
P(4pB,) = [1 — P(4,| BYP) =F. 


It follows from (12) that 


1 a ae aX 2F 9 
Re toe. |: ioe. — + doa. — -- —__ ae 
4B m5 te OB25 + g loses + 75 logs, © 0.12 


is the average amount of weather conveyed by a prediction. In the case of 
100% accurate predictions, A, = B,, A, = B, and (12) reduces to 


1 1 4 4 
ee ie ee. 
AB 5 082 5 5 0825 672 


PROBLEMS 


1. Which conveys more information, a message telling a stranger’s birthday or 
a message telling his telephone number ? 
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2. Find the average amount of information in bits of a message telling whether 
or not the outcome of throwing a pair of unbiased dice is 
a) Anoddnumber; b) Aprimenumber; c) A number no greater than 5. 


3. An experiment has four possible outcomes, with probabilities a 


16? 16 
and ;%, respectively. What is the average amount of uncertainty about the 


outcome of the experiment ? 


4. Each of the signals 4,,..., A, has equal probability of being transmitted 
over a communication channel. In the absence of noise, the signal A; is received 
asa;(j =1,...,n), while in the presence of noise A; has probability p of being 
received as a; and equal probability of being received as any of the other symbols. 
What is the average amount of information about the symbols A,,..., An 
conveyed by receiving one of the signals a,,..., a, 

a) In the absence of noise; b) In the presence of noise? 


ee 
n=1- 


Ans. a) log,n; b) loggn + plogsp + (1 — p) log, 


Appendix 2 


GAME THEORY 


Consider the following simple model of a game played repeatedly by two 
players. Each player can choose one of two strategies determining the result 
of the game. The interests of the players are completely conflicting, e.g., 
whatever one player wins, the other loses.? Such a 
“two-person game’? can be described by the table 
shown in Figure 9, where the quantity in the ith row V4 Vi0 
and jth column is the amount gained by the first player 
if he chooses strategy i while his opponent chooses 
strategy j (i, 7 = 1, 2). For example, v,, is the amount Vo4 Vo0 
gained by the first player (the first player’s “payoff’’) if 
he chooses the first strategy and his opponent (the 
second player) chooses the second strategy, while —v., FIGURE 9 
is the second player’s payoff if he chooses strategy 1 and 
his opponent (the first player) chooses strategy 2. It is 
now natural to ask for each player’s “‘optimal strategy.” 

This question is easily answered in the case where 


MiN (041, Vy2) > MAX (Vg, Voq), (1) 


say, since then regardless of how the second player acts, the first player 


1 More generally, a game of strategy involves more than two players, each with more 
than two available strategies, but the essential features of game theory (in particular, its 
connection with probability) emerges even in this extremely simple case. 

2 Such a game, in which the algebraic sum of the players’ winnings is zero, is called a 
zero-sum game. 
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should always choose the first strategy, thereby guaranteeing himself a gain 
of at least 
Min (V1, V2). 


Assuming a “clever”? opponent, the second player should then choose the 
strategy which minimizes the first player’s maximum gain, 1.e., the strategy j 
such that 

Vy; = MIN (V44, Vy,). 


The case just described is atypical. Usually, a relation like (1) does not 
hold, and each player should adopt a “mixed strategy,’ sometimes choosing 
one of the two “pure strategies’’ available to him and sometimes choosing the 
other, with definite probabilities (found in a way to be discussed). More 
exactly, the first player should choose the ith strategy with probability p,, 
while the second player should (independently) choose the jth strategy with 
probability p,;. Then the first player’s strategy is described by a probability 
distribution P; = {p4;, Piz}, while the second player’s strategy is described 
by a probability distribution P, = {Po, Poo}. If these mixed strategies are 
adopted, the average gain to the first player is clearly just 


2 
V(P, Po) = 2, MisPusPay (2) 
41,j= 


Suppose the second player makes the optimal response to each strategy 
P, = {pi1; Piz} chosen by the first player, by adopting the strategy P} = 
{p>,> Po} Minimizing the first player’s gain. The first player then wins an 
amount 
V(Pi, PF) = min V(P,, P») = Vi(P3) 
Pe 


if he chooses the strategy P,. To maximize this gain, the first player should 
choose the strategy P® = {p°,, p®,} such that 


V,(Py) = max V,(P)), 
Py 


always, of course, under the assumption that his opponent plays in the best 
possible way. Exactly the same argument can be applied to the second player, 
and shows that his optimal strategy, guaranteeing his maximum average gain 
under the assumption of optimal play on the part of his opponent, is the 
strategy P} = {p2,. Poo} Such that 


V(P2) = Bee pee 2)» 
where 
Ve(P2) = ies {—V(P1, Pe)}. 
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To calculate the optimal strategies P? and P$, we consider the function 


V(x, y) = Vyxy + Viox(1 — y) + ¥_(l — x)y +2) — xsd = y)s 


which for x = p,, and y = pz, equals the average gain of the first player if 
the mixed strategies P, = {py, P12} and P, = {Po1, Poo} are chosen. The 
function V(x, y) is linear in each of the variables x and y, 0 < x,y < 1. 
Hence. for every fixed x, V(x, y) achieves its minimum V,(x) at one of the 
end points of the interval 0 < y < l,ie., fory=Oory=1: 


V(x) = min V(x, y) = min {v,.x + v1 — x), 0y,x + v1(1 — x)}. 
y 


As shown in Figure 10, the graph of the function V(x) is a broken line with 


VrVUWyXt+Vo4(1-x) 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


FiGurE 10. A case where min (v,;, Vi2) < max (Ve, Vee). 
vertex at the point x° such that 
VyaX° + Vo9(1 — x°) = 044x° + vq, (1 — x), 


i.e., at the point 


0 Vog — Vo (3) 
Viz + Veg — (Vig + V9q) 


The value x = x® for which the function V,(x), 0 < x < 1 takes its 
maximum is just the probability p°®, with which the first player should 
choose his first pure strategy. The corresponding optimal mixed strategy 
P) = {p?,,P{.} guarantees the maximum average gain for the first player 
under the assumption of optimal play on the part of his opponent. This 
gain is 

V(x?) = 041X° + Ve] — x) = dy_x° + Vo9(1-— x°). (4) 
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Moreover, (4) implies 
V(x®, y) = y[virx® + vy(1 — x°)] + GU — y)[Piex° + Voo(1 — x°)] = Vy(x°) 


for any y in the interval 0 < y < 1. Hence, by choosing p?, = x°, the first 
player guarantees himself an average gain V,(x°) regardless of the value of 
y, ie., regardless of how his opponent plays. However, if the first player 
deviates from this optimal strategy, by choosing p,, = x #x°, then his 
opponent need only choose p2, = y equal to 0 or 1 (as the case may be) to 
reduce the first player’s average gain to just V,(x). 

Applying the same considerations to the second player, we find that the 
second player’s optimal strategy is such that p}, = y°, where 


= Vyy + Veg — (Vy2 + Ve) 
[(5) is obtained from (3) by reversing the roles of players 1 and 2, Le., by 
interchanging the indices 1 and 2]. As in the case of the first player, this 
choice guarantees the second player an average gain V,(y°) regardless of the 
first player’s strategy, i.e., 

—V(x,y)= VA"), O<x<l. (6) 

In particular, (4) and (6) imply 

V,(x°) = VX", y"), 

Vy") = =e 


Example 1. One player repeatedly hides either a dime or a quarter, and 
the second player guesses which coin is hidden. If he guesses properly, he 
gets the coin, but otherwise he must pay the first player 15 cents. Find both 
players’ optimal strategies. 


0 Vey =~ Ue (5) 


Solution. Here 
yy = i, Vi2 = iF 
Vo, = 15, Veg = —25, 
so that, by (3), 
ee 
—35—30 13. 
Therefore the first player should hide the dime with probability #5, and hide 
the quarter with probability 7;.° Similarly, by (5), 
Pee 2 


py = = 2 = 


See 


SS 
Pi = x 


® For the first player, hiding the dime is (pure) strategy 1, and hiding the quarter is 
strategy 2. For the second player, guessing the dime is strategy 1, and guessing the quarter is 
strategy 2. 
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and hence the second player should guess that the hidden coin is a dime with 
probability <3, and that it is a quarter with probability 3°5. Then, according 
to (4). the first player’s average gain will be 


8 5 5 
Pre aA) {Q:- Mbsirs 2h sire pose 
a) ee 5; ce 


while the second player’s average gain will be 
5 

—V(x°, y°) = : 
== 


Thus this game is unfavorable to the first player, who loses an average of 
zs cents every time he plays, even if he adopts the optimal strategy. However, 
any departure from the optimal strategy will lead to an even greater loss, if 
his opponent responds properly. 


Example 2 (Aerial warfare). White repeatedly sends two-plane missions 
to attack one of Blue’s installations. One plane carries bombs, and the 
other (identical in appearance) flies cover for the plane carrying the bombs. 
Suppose the lead plane can be defended better by the guns of the plane in the 
second position than vice versa, so that the chance of the lead plane surviving 
an attack by Blue’s fighter is 80%, while the chance of the plane in the second 
position surviving such an attack is only 60%. Suppose further that Blue 
can attack just one of White’s planes and that Blue’s sole concern is the 
protection of his installation, while White’s sole concern is the destruction 
of Blue’s installation. Which of White’s planes should carry the bombs, and 
which plane should Blue attack? 


Solution. Let White’s payoff be the probability of accomplishing the 
mission. Then°® 


V1, = 0.8, Vip = I, 
Vo = 1, Voo = 0.6, 
and hence 
—0.4 2 —0.2 1 
Oo = x® = — =, Q = y? = senna 
= ee Se 


by (3) and (5). Thus always putting the bombs in the lead plane is not 
White’s best strategy, although this plane is less likely to be shot down than 


* After J. D. Williams, The Compleat Strategyst, McGraw-Hill Book Co., Inc., New 
York (1954), p. 47. 

* For White, putting the bombs in the lead plane is (pure) strategy 1, and putting the 
bombs in the other plane is strategy 2. For Blue, attacking the lead plane is strategy 1, 
and attacking the other plane is strategy 2. 
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the other. In fact, if White always puts the bombs in the lead plane, then 
Blue will always attack this plane and the resulting probability of the mission 
succeeding will be 0.8. On the other hand, if White adopts the optimal mixed 
strategy and puts the bombs in the lead plane only two times out of three, 
he will increase his probability of accomplishing the mission by ;’s, since, 
according to (4), 


a 13 
Vix8 WN ice te — _, 
ke ee +; 15 
By the same token, Blue’s best strategy is to attack the lead plane only one 
time out of three and the other plane the rest of the time. 


PROBLEMS 


1. Prove that the game considered in Example 1 becomes favorable to the first 
player if the second player’s penalty for incorrect guessing is raised to 20 cents. 


2. In Example 1, let a be the second player’s penalty for incorrect guessing. For 
what value of a does the game become “‘fair’’? 


3. Blue has two installations, only one of which he can successfully defend, 
while White can attack either but not both of Blue’s installations. Find the 
optimal strategies for White and Blue if one of the installations is three times as 
valuable as the other.® 


Ans. White should attack the less valuable installation 3 out of 4 times, 
while Blue should defend the more valuable installation 3 out of 4 times. 


® After J. D. Williams, op. cit., p. 51. 
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BRANCHING PROCESSES 


Consider a group of particles, each “randomly producing”’ more particles 
of the same type by the following process: 


a) The probability that each of the particles originally present at some 
time t = 0 produces a group of k particles after a time f is given by 
p(t), where k = 0,1, 2,... and p,(t) is the same for all the particles.’ 

b) The behavior of each particle is independent of the behavior of the 
other particles and of the events prior to the initial time t = 0. 


A random process described by this model is called a branching process. 
As concrete examples of such processes, think of nuclear chain reactions, 
survival of family names, etc.” 

Let £(t) be the total number of particles present at time ¢. Then &(f) is a 
Markov process (why?). Suppose there are exactly k particles initially 
present at time ¢ = 0, and let &,(t) be the number of particles produced by 
the ith particle after a time ¢. Then clearly 


E(t) = E(t) + -°* + &(t), (1) 


where the random variables &,(t),..., &,(t) are independent and have the 
same probability distribution 


PL) = 2) = p,(d. H= D1 2... 


1 The case k = 0 corresponds to “‘annihilation” of a particle. 
2 Concerning these examples and others, see W. Feller, op. cit., p. 294. 
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Let Pyn(t) be the probability of the k particles giving rise to a total of n 
particles after time ¢, so that the numbers p,,(t) are the transition proba- 
bilities of the Markov process &(t), and introduce the generating functions® 


F(z) => py", (2) 
Fu(t, 2) = > Peal” 3) 


Suppose the probability of a single particle giving rise to a total of n particles 
in a small time interval At is 

Pr(At) = 4, At + o(At), 
while the probability of the particle remaining unchanged is 


pri(At) = 1 — At + OAD). 
Moreover, let 


so that 
pa A, = 0. (4) 


Then the Kolmogorov equations (8.15), p. 105 for the transition probabilities 
P,At) = Pi,(t) become 


d 
— Pat) = Pel), n=0,1,2,... 
dt ke 

Next we deduce a corresponding differential equation for the generating 
function F(t, z). Clearly 
£ Ft) =£ Yaz" =z" {p=TAY Pld”) 

d t = d t n=0 n=0 d t ke n=0 

(justify the term-by-term differentiation), where F,(t,z) is the generating 
function of the random variable &(t) for the case of k original particles.’ 
But, according to (1), &(t) is the sum of k independent random variables, each 
with generating function F(t, z). Therefore, by formula (6.7), p. 71, 


F,(t, Zz) = Pt rg k=0, Po ee (6) 
(the formula is trivial for k = 0). Substituting (6) into (5), we get 


. F(t, z) = SM lF(t DF. (7) 


8 Note that F,(¢, z) = F(t, z), since clearly pin(t) = pn(t). 
* Clearly Fy(z) = 1, since new particles cannot be created in the absence of any original 
particles. 
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In what follows, we will assume that a given branching process &(f) is 
specified by giving the transition densities A,, k = 0,1,2,... Let f(x) be 
the function defined by the power series 


fx) = nx (8) 


so that in particular f(x) is analytic for 0 < x < 1. Then, according to (7), 
the generating function F(t, z) satisfies a differential equation of the form 


ee - 
eed (9) 


Moreover, since F(0,z) =z, the generating function F(t, z) coincides for 
every z in the interval 0 < z < 1 with the solution x = x(t) of (9) satisfying 
the initial condition a 

x(O) = z. (10) 
Instead of (9), it is often convenient to consider the equivalent differential 
equation 

as << 

gx Fix) 
for the inverse ¢ = f(x) of the function x = x(t). The function satisfying 
(11) and the initial condition (10) is just 


Example 1. If 


then 
F(x) = ACL — x) 


and 


= =| a=; ond —»)- In (1 — z)]. 


Hence F = F(t, z) is such that 


In (1 — F) = —At+ In (1 — 2), 
es, 


F(t, z) = 1—e(1 — 2). 


The probabilities p,,(t) are found from the expansion 


F(t, 2) => p,(2", 
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which in this case implies 


P(it)=1l—e™“, pr=e 
Pit) = 0, . Be Zaye. 


oat 


Example 2. If 
Xo = 0, Ai cal —1, 


eee ee 
(k — 1)k 
then 
(oe) xk co xk 
Ss ce pax = Ose 
== > XIn(l —x)-- teh ee — 2), 
and hence 
ao of du Se 
2 f(u) @ ==) In (1 = In(1—-z) y 


= —InIn(1 — x) + InIn(1 — z). 
It follows that F = F(t, z) is such that 


In(l =e. = 
: In (1 — 2) 
Le, 
FG, 2) = 1 -—-de. 
To find the corresponding probabilities p,,(t), we use repeated differentiation: 


Plt) = 4, Pit) = e* 
= <a = = ~(- i =n + 2; 
a oe Ee 


Turning to the analysis of the differential equation (9), where f(x) is 
given by (8), we note that 


fi’) = Dkk — Dix > Ot re 1. 


Therefore f(x) is concave upward in the interval 0 < x < 1, witha monotoni- 
cally increasing derivative. Because of (4), x = 1 is a root of the equation 
f(x) = 0. This equation can have at most one other root x = «(0 < « < 1). 
Thus f(x) must behave in one of the two ways shown in Figure 11. 

We now study the more complicated case, where f(x) = 0 has two roots 
x=a (0<a<1) and x =1, corresponding to two singular integral 
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FIGURE 11 


curves x(t)=« and x(t)=1 of the differential equations (9) and (11). 
Consider the integral curve 


* da 
z f(u) 


going through the point t=0, x =z (0<z< 4). Since the derivative 
f'(a) is finite and f(x) ~ f’(x)(x — «) for x ~ «, the value of ¢ along the 
integral curve (12) increases without limit as x — «, but the curve itself never 
intersects the other integral curve x(t) = «. The function /(x) is positive in the 
interval 0 < x < a, and hence the integral curve x = x(t) increases mono- 
tonically as t > oo, remaining bounded by the value x = «. Being a bounded 
monotonic function, x(t) has a limit 


= (12) 


So == lm x2), ec Oe es 


to 
But f(x) approaches a limit f(6) as x > 8, Le., 
f(p\) = ee [x(t] sept x'(t), 
where /(8) must vanish, since otherwise the function 


x) =2+f Fah as 


would increase without limit as t —~ oo. It follows that 6 is a root of the 
equation f(x) = 0, and hence must coincide with «. Therefore all the integral 
curves x = x(t) going through the point x =z, 0<z<a« for t=0 in- 
crease monotonically as t — oo and satisfy the condition 


lim x(t) = «. (13) 


to 
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The behavior of the integral curves going through the point x = z,a<z< 1 
for t= 0 is entirely analogous. The only difference is that x(t) now de- 
creases monotonically, since the derivative x’(t) = f[x(t)] is negative and 
f(x) <= 0 fer @ = = 1. The behavior 
x of typical integral curves in the interval 
0<z<1 is shown in Figure 12, 
: where 0 < 2, =< 4 =. = 1. 
22 The behavior of the integral curves 
at z= 1 warrants special discussion. 
First we note that in any case x(t) = 1 1s 
an integral curve corresponding to z = 
Z 1. Suppose 


FIGURE 12 xo f(x) ee: ae 


for some Xy, ® < X) < 1.5 Then an arbitrary integral curve of the form 


t=ty) + 0 <= yy =a, (15) 
: = Fea’ 
going through some point (f9, Xo), decreases without limit as x — 1, L.e., 
due 
== ft, = —> —O 
: xo f (u) 


as x > 1. This shows that given any f¢) > 0, the equation 


t(z) = to + . —- 


holds for some x = z, « < z < 1. Hence every integral curve intersects the 
axis t = 0 in a point (0, z) such that « < z < 1 (see Figure 13). It follows 
that in this case x(t) = 1 is the unique integral curve going through the point 


(0, 1). 
On the other hand, suppose 
a 
xo fi (x) 


Then for sufficiently large fy, the integral curve (15) intersects the integral 
curve x(t) = 1, and is in fact tangent to it at the point (t, 1) where 


— 0. : (16) 


® This is always the case if f’(1) < © (why ?). 
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(see Figure 13). Inthiscase, thereis x 

a whole family of integral curves x(t) 

going through the point (0,1), where 1 T 

each x_(f) is parameterized by the Xo(T) 
appropriate value of t > 0. Among 

these integral curves, the curve x9(t) ° 

shown in the figure, corresponding to 

the value + = 0, has the property of 

lying below all the other integral rt 
curves, 1.e., | Scie 94 


meat) =< x.(t), eee eee 3 


This is explained by the fact that the solution of our differential equation is 
unigue in the region 0<x<1,0<1t< ©, so that the integral curves 
do not intersect in this region. It is also easy to see that the integral curve 
x,(t) is the limit of the integral curves x(t, z) lying below it and passing 
through points (0, z) such that 0 < z < 1. In other words,° 
xi = cad ME, 2). (17) 
The above analysis of the differential equation (9) has some interesting 
implications for the corresponding branching process E(t). In general, there 
is a positive probability that no particles at all are present ata given time f¢. 
Naturally, this cannot happen if A) = 0, since then particles can only be 
“created” but not “annihilated.’’ Clearly, the probability of all particles 
having disappeared after time ¢ 1s 


Polt) ss 7 0) 
if there is only one particle originally present at time ¢ = 0, and 
Prolt) = LF(t, 0))* = [po(t)I? 


if there are k particles at time t = 0. The function p(t) is the solution of the 
differential equation (9) corresponding to the parameter z = 0: 


dpo(t) _ 

dt 
As already shown, this solution asymptotically approaches some value py = « 
as t > 00, where « is the smaller root of the equation f(x) = 0 [recall (13)]. 
Thus py) = « is the extinction probability of the branching process i e.. 
the probability that all the particles will eventually disappear. If the function 
f(x) is positive in the whole interval 0 < x < 1, the extinction probability 
equals 1. 


flpo®], po(O) = 0. 


6 Note that x(t, z) = F(t, z) fort >0,0<z<1. 
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There is also the possibility of an “explosion”? in which infinitely many 
particles are created. The probability of an explosion occurring by time f is 
just 


Po(t)=1— P {&(t) < o} = 1 — SP {&(t) = n} 


=1— p,() =1- lim F(t, 2). 


In the case where x(t) = 1 is the unique integral curve of (9) passing through 
the point (0, 1), we clearly have 


lim FG, z) = 4: 
Z2>1 


Therefore p,,(t) = 0 for arbitrary t if (14) holds, and the probability of an 
explosion ever occurring is 0. However, if (16) holds, we have (17) where 
Xo(t) is the limiting integral curve described above and shown in Figure 13. 
In this case, 


Po(t) =1— Xo(t) > 0 


and there is a positive probability of an explosion occurring. 


PROBLEMS 


1. A cosmic ray shower is initiated by a single particle entering the earth’s 
atmosphere. Find the probability p,(t) of n particles being present after time tf 
if the probability of each particle producing a new particle in a small time 
interval At is AAt + o(At). 

Hint. Ay = at Xo — A. 

Ans. py(t) = e**(1 — et)n-1, y > 1, 
2. Solve Problem 1 if each particle has probability }Ar + o(Ar) of producing a 


new particle and probability wAr + o(Ar) of being annihilated in a small time 
interval At. 


Hint. dy = 4, = —A+t B), Ag =H. 


Ans. polt) = vy; Prlt) = (1 — Ay) — py)Qy)" (n> 1), 
where 


1 — e@-vit 
uw — AeA-w)t ee 
- = 
t : 
i oh = a, 
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3. Find the extinction probability py of the branching process in the preceding 
problem. 


. 
Ans. Po =\A 
1 


Appendix 4 


PROBLEMS OF OPTIMAL CONTROL 


As in Sec. 15, consider a physical system which randomly changes its 
state at the times ¢ = 1, 2,... , starting from some initial state at time t = 0. 
Let ¢,, €,... be the possible states of the system, and &(t) the state of the 
system at time f, so that the evolution of the system in time is described by the 
consecutive transitions 


6(0) + €() > EG =, 


We will assume that &(t) is a Markov chain, whose transition probabilities 
Piss 1,7 = 1,2,... depend on a “control parameter” chosen step by step by 
an external “operator.” More exactly, if the system is in state ¢, at any time 
n and if d is the value of the control parameter chosen by the operator, then 
Pis = Pii(@) 
is the probability of the system going into the state ¢, at the next step. The 
set of all possible values of the control parameter d will be denoted by D. 
We now pose the problem of controlling this “guided random process” 
by bringing the system into a definite state, or more generally into one of a 
given set of states E, after a given number of steps n. Since the evolution of 
the process €(t) depends not only on the control exerted by the operator, but 
also on chance, there is usually only a definite probability P of bringing the 
system into one of the states of the set Z, where P depends on the “control 
program” adopted by the operator. We will assume that every such control 
program consists in specifying in advance, for all ¢; andt=0,...,n—1, 
the parameter 
ad == Hte,,f) 
136 
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to be chosen if the system is in the state ¢; at the time ¢. In other words, the 
whole control program is described by a decision rule, i.e., a function of two 
variables 

og Sg gS 


where x ranges over the states ¢,, «,,... and ¢ over the times 0,...,” — l. 
Thus the probability of the system going into the state <, at time k + 1, 
given that it is in the state <, at time k, is given by 


Pu = Pis(Q), d= d(e,, k). 


By the same token, the probability of the system being guided into one of the 
states in E depends on the choice of the control program, i.e., on the decision 
rule d = d(x, t), so that 

P == P(d). 


Control with a decision rule d° = d°(x, t) will be called optimal if 
P(d°) = max P(d), 
d 


where the maximum is taken with respect to all possible control programs, 
i.e., all possible decision rules d = d(x, t). Our problem will be to find this 
optimal decision rule d°, thereby maximizing the probability 


P(d) = P {&(n) € E} 


of the system ending up in one of the states of EF after n steps. 
We now describe a multistage procedure for finding d°. Let 


P(k, i, d) = P {&(n) € E| E(k) = ¢;} 


be the probability that after occupying the state ¢, at the kth step, the system 
will end up in one of the states of the set E after the remaining n — k steps (it 
is assumed that some original choice of the decision rule d = d(x, t) has 
been made). Then clearly 


This is a simple consequence of the total probability formula, since at the 
(kK + 1)st step the system goes into the state «,; with probability p,,(d), 
d = d(<,;, k), whence with probability P(A + 1, 7, d) it moves on (n —k — 1 
steps later) to one of the states in the set E. 

For k = n — 1, formula (1) involves the probability 


P(n,j,d) l if <¢,e8, 2) 
nN, ]; a 

: 0 otherwise, ( 
and hence 
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where the summation is over all j such that the state <, belongs to the given 
set E. Obviously, P(n — 1,i, d) does not depend on values of the control 
parameter other than the values d(<,,n — 1) chosen at the time n — 1. 
Letting d° denote the value of the control parameter at which the function 
(3) takes its maximum,! we have 


Pn — 1, i) = P(n — 1, i, d°) = max P(n — 1, i, d). (4) 
deD 


Clearly, there is a value d° = d%e,,n — 1) corresponding to every pair 
(eb), hae 7, 2 
For k = n — 2, formula (1) becomes 


P(n aia es i, d) = > p,,(d)P(n so Led; d). 
J 


Here the probabilities p,,(d) depend only on the values d = d(e,, n — 2) of 
the decision rule d = d(x, t) chosen at time n — 2, while the probabilities 
P(n — 1,j,d) depend only on the values d= d(e,,n — 1) chosen at time 
n — 1. Suppose we “‘correct’”’ the decision rule d = d(x, t) by replacing the 
original values d(e;,n — 1) by the values d%(e,, nm — 1) just found. Then the 
corresponding probabilities P(n — 1, j, d) increase to their maximum values 
Pn — 1,7), thereby increasing the probability P(n — 2,i,d) to the value 


P(n — 2,i,d) = & p,(@)P\n — 1.) (5) 


Clearly, (5) depends on the decision rule d = d(t, x) only through the de- 
pendence of the transition probabilities p,;(d) on the values d = d(e;, n — 2) 
of the control parameter at time n — 2. Again letting d° denote the value of 
the control parameter at which the function (5) takes its maximum, we have 


Pn — 2, i) = P(n — 2, i,.d°) = max Pm — 2, i, d). 
deD 


As before, there is a value d° = d%e,,n — 2) corresponding to every pair 


(¢,,n — 2),i=1,2,... Suppose we “‘correct’’ the decision rule d(x, t) by 
setting 

A(x, t) = d(x, t) (6) 
for t=n—2,n—1 and all x = &, e,... Then clearly the probabilities 


P(k, i, d) take their maximum values Pk, i)fori=1,2,...andk =n — 2, 
n — 1. Correspondingly, formula (1) becomes 


and this function of the control parameter d takes its maximum for some 
d° = d%(e,,n — 3). We can then, once again, “‘correct’”’ the decision rule 


1 It will be assumed that this maximum and the others considered below exist. 
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d = d(x, t) by requiring (6) to hold for t = n — 3 and all x = ¢,, ¢,,..., a8 
well as fort =n—2,n—landallx=¢,,6,... 

Continuing this step-by-step procedure, after n — 1 steps we eventually 
get the optimal decision rule d = d°(x, t), defined for t = 0,...,n — 1 and 
all x = ¢,, €,,..., such that the probability P(d) = P(0, i, d) satisfying the 
initial condition €(0) = <; achieves its maximum value. At the (n — k)th 
step of this procedure of “successive corrections,’’ we find the value d®° = 
d°(<;, k) maximizing the function 


where P°(k + 1,7) is the maximum value of the probability P(k + 1, j, d). 
Carrying out this maximization, we get Bellman’s equation? 


P%(k, i) = max > p;(d)P(k + 1, j), 
deD 4 
which summarizes the whole procedure just described. 


Example 1. Suppose there are just two states <, and ¢,, and suppose 
the transition probabilities are continuous functions of the control parameter 
in the intervals 


1 < py(d) < B,, Xe < Po(d) < fo. 


What is the optimal decision rule maximizing the probability of the system, 
initially in the state ¢,, going into the state ¢, two steps later? 


Solution. In this case, 


PU, 1) ae B,, PU, 2) = B., 
Piel) — _ [Pis(Z)B1 + Pro(d)Be] = Sag [P1(4)(B r= By) tr: Bal. 


If the system is initially in the state ¢,, then clearly we should maximize the 
transition probability p,, (by choosing p,,; = B,) if 8, > B., while maximizing 
the transition probability p,. = 1 — py, (by choosing p,,; = «,) if By < By.3 
There is an analogous optimal decision rule for the case where the initial 
state of the system is €. 


Example 2 (The optimal choice problem). Once again we consider the 
optimal choice problem studied on pp. 28-29 and 86-87, corresponding to 


* In keeping with (2)—-(4), we have 

P\(n,;) = 1 if ¢,E€E, 
0 otherwise. 
* Clearly, any choice of p,, in the interval «, < pis < B, is optimal if 8, = By. 
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a Markov process &(t) with transition probabilities 


0 a = = 
i 
perme seo ee (7) 
Ps =\G- Di 
ES fj mot 
m 


where, as on p. 28, choice of an object better than all those previously 
inspected causes the process &(0)— &(1) > &(2) >:-- to terminate. In 
each of the states ¢,,...,¢,, (whose meaning is explained on p. 86), the 
observer decides whether to terminate or to continue the process of inspection. 
The decision to terminate, if taken in the state ¢,, is described formally by 
the transition probabilities 
1 if “$= j, : 

Pie dt sae 29 (8) 
while the decision to continue corresponds to the transition probabilities (7): 
Hence we are dealing with a “guided Markov process,’’ whose transition 
probabilities p;; depend on the observer’s decision. Here the control param- 
eter d takes only two values, 0 and 1 say, where 0 corresponds to stopping 
the process and | to continuing it. Thus (8) gives the probabilities p,,(0) and 
(7) the probabilities p,,;(1). 

Every inspection plan is described as a decision rule d= d(x), x = 
E1,-++ 5m, Which specifies in advance for each of the states ¢,...,&, 
whether inspection should be continued or terminated by selecting the last 
inspected object. The problem consists of finding an inspection plan, or 
equivalently a decision rule d=d(x), x =&,...,€m, maximizing the 
probability of selecting the very best of all m objects. This probability is just 


where i/m is the probability that the ith inspected object is the best (recall p. 
29), p,; is the probability that the process will stop in the state <,, and the 
summation is over all the states <; in which the decision rule d = d(x) calls for 
the process to stop. 

To find the optimal decision rule d° = d°(x) maximizing (9), we consider 
the probability P(k, d) of selecting the best object, given that the number of 
previously inspected objects is no less than k, i.e., given that the process 
&(t) actually occupies the state ¢,. By the total probability formula, we have 


P(k, d) = PrP, a (10) 
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Clearly, if the process occupies the state ¢,,, then the mth object is the best of 
the first m objects inspected and hence automatically the best of all 
objects. Therefore the optimal value of the decision rule d= d(x) for 
x =¢,, is just d%e,,) = 0, and P(m, d) = 1 for this value. It follows from 
(9) and (10) that 


mls ahaa) = 0, 
m 
an — 1, d) — (11) 
ae | 
——__——— if d(e,,)=1 
(m ne 1)m 1 c= 1) 


is the probability of choosing the best object, given that the process stops in 
the state <,, and the number of previously inspected objects is no less than 
m — 1. Moreover, (11) implies that the optimal value of the decision rule 
eee tors — <., is d%(c,,_,) = 0, and that 


Pim —-)j=-———— 
m 


Now suppose the optimum values of the decision rule d = d(x) are all 
zero for x = &,..., &, corresponding to the fact that the process is termi- 
nated in any of the states ¢,,..., €,. Then what is the optimal value d%«,_,)? 
To answer this question, we note that (9) and (10) imply that 


P(k — 1, d) 
—- if d(e,s) =0, 
m 
Meek 1 ki 4 
(k—1)km k(k+1) m (m—1)m 


oe ee 


is the probability of choosing the best object, given that the process stops 
in the states ¢,,...,¢,, and the number of previously inspected objects is 
no less than k — 1. It follows that the optimal value of the decision rule 
= eae) for xX = «,_, is 


1 
0 if ——+- a +-$——— <1, 
a"(&1) = = == m—1 (12) 


1 otherwise. 


Moreover, it is easy to see that the optimal decision rule d® = d(x) has the 
structure 
0 = oe eee 
d°(x) = cg 


1 | ae ee eres 
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where my is some integer. Thus the optimal selection procedure consists 
in continuing inspection until the appearance of an object numbered k > my 
which is better than all previously inspected objects. According to (12), 
My is the largest positive integer such that 


1 1 1 
mot = “Ss * a ee 13 
Mo Mo oh 1 nh: 1 ( ) 
PROBLEMS 
1. In Example 2, prove that 
m 


if m is large, where e = 2.718 ... is the base of the natural logarithms. 
Hint. Use an integral to estimate the left-hand side of (13). 


2. Find the exact value of my for m = 50. Compare the result with (14). 


3. Consider a Markoy chain with two states ¢, and ¢, and transition proba- 
bilities p;;(d) depending on a control parameter d taking only two values 0 and 1. 
Suppose 


PuO) =§ Pu) = $3 Pull) =%, Pall) =%. 
What is the optimal decision rule maximizing the probability of the system 


initially in the state ¢,, going into the state ¢, three steps later? What is this 
maximum probability ? 
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